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The outlines of a theory of elastic loss in cross-linked polymers are developed on the basis 
of the ideas underlying the theory of dielectric loss of Fuoss and the writer and the three- 
dimensional network model of rubber structure proposed by James and Guth. The relaxation 
time distribution function for elastic loss is given in the approximation of free internal rotation 
of chain segments of the net. The influence of intermolecular hindering torques and network 
restraints on the relaxation time spectrum is discussed, but explicit calculations are postponed 


for a later report. 





HE statistical theory of the equilibrium 

elastic behavior of rubber has contributed 
greatly to our understanding of the structure of 
cross-linked polymers, although from the quanti- 
tative point of view it still leaves something to 
be desired. The basic assumption of the theory 
is that the resistance offered to deformation by 
a cross-linked polymer network is principally 
due to the decrease in configurational entropy 
attending the elongation of the chain segments 
of the net. 

If a time-dependent stress is applied to a 
polymer network, the finite time of response of 
the chain segments leads to a phase difference 
between strain and stress and a resultant energy 
dissipation in the specimen known as elastic loss. 
The elongation response of a chain segment of 
the network to the applied stress is accomplished 
by internal rotatory Brownian motion, which is 
associated with a relaxation time spectrum of the 
type described by Fuoss and the writer! in their 
theory of dielectric loss in polar polymers. It is 


aodiy Kirkwood and R. M. Fuoss, J. Chem. Phys. 9, 329 
41). 


our purpose in the present article to sketch the 
outlines of a theory of elastic loss in cross-linked 
polymers on the basis of the ideas underlying 
the theory of dielectric loss. 

As a preliminary step in the formulation of 
the theory of elastic loss, it is necessary to adopt 
a model for the polymer network spanning a 
specimen of a rubber-like material. Following 
the early qualitative theories of rubber elasticity,’ 
two types of network model have been proposed 
in efforts to refine the theory. The one is the 
three-dimensional network of James and Guth*® 
and the other is the model of Wall and Flory.* 
While it is not our purpose here to make a 
judgment on the relative merits of the two 
models in the equilibrium case, we have decided 
that the more detailed structural model of James 
and Guth offers the greater insight into the 
mechanism of the irreversible deformation of 


2 E. Guth and H. Mark, Monats. f. Chemie 65, 93 (1934); 
Kuhn, Kolloid Zeits. 76, 258 (1936). 

”3} M. James and E. Guth, J. Chem. Phys. 11, 455 
(1943). 

4F. T. Wall, J. Chem. Phys. 10, 132, 485 (1942); P. Flory, 
Chem. Rev. 35, 51 (1944); L. R. G. Treloar, Trans. Fara- 
day Soc. 40, 59 (1944). 
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the network associated with internal rotatory 
Brownian motion of the chain segments. At the 
same time, we recognize the possibility of de- 
veloping a theory of elastic loss on the basis of 
the Flory-Wall model along the lines suggested 
by Alfrey.§ 

Before proceeding to the formulation of the 
theory, we must remark that it is not designed 
to account for all of the elastic loss, but only 
that part in excess of frequency independent loss 
associated with the permanent rupture of cross- 
linking bonds and with chemical alteration of the 
specimen. There is, of course, no experimental 
difficulty in separating the two types of loss for 
independent study. 

We shall consider a rubber consisting of a 
cross-linked polymer network, possibly imbedded 
in a matrix of polymer chains not participating 
in the net structure. Following James and Guth,’ 
we replace the irregular cross-linked network by 
a regular cubic array, each lattice point of which 
is connected to its nearest neighbors by a polymer 
chain segment composed of » bonds. We shall 
refer to each skeletal atom of the chain together 
with side chain groups bonded to it as a monomer 
unit. We denote the lattice parameter of the 
network by a, the number of chains per unit 
cross-section by v, and the number, va, of chains 
per unit length by M. 

When a stress o is applied, let us say in the 
z direction, parallel to the axis of a thin cylindri- 
cal specimen of the polymeric material, we shall 
assume that it is distributed on the microscopic 
scale uniformly among the » lattice points per 
unit area in each cross-sectional plane of the 
network perpendicular to the axis of the speci- 
men.* Under this assumption, the terminal ele- 
ments of each chain segment between successive 


5 T. Alfrey, J. Chem. Phys. 12, 374 (1944). 

* This assumption is consistent with the perhaps over- 
simplified model of rubber structure of James and Guth, 
provided inertial effects are neglected. Equation (1), of 
course, represents the average terminal force on a chain 
segment arising from the applied stress, from which there 
will be fluctuations. However, we shall be interested only 
in the average force in the subsequent analysis. 

A more drastic assumption presently to be made is that 
intermolecular hindering torques, acting on internal ele- 
ments of a chain segment, have their equilibrium values 
even when the specimen is subjected to a time-dependent 
external stress. This assumption, equivalent to neglecting 
intermolecular relaxation effects, should be subjected to a 
more careful investigation, in future attempts to refine the 
preliminary theory to be described here. 


network planes are subjected to forces +F and 
— F, respectively, in the.z direction, where 


F=a/v. (1) 


Under the action of the terminal forces the 
length z of the chain segment changes by internal 
rotatory Brownian motion, and a strain €., is 
developed. 

€.,=Az/a, 

Az=z2—a. 


As the length of the chain segments oriented in 
the z direction increases, the average force 
exerted by them on the chains oriented in the x 
and y directions is altered in such a manner that 
the latter contract. The resulting Poisson effect 
in rubbers is in general adequately described by 
the assumption of incompressibility and the 
strains €,, and €,, are approximately equal to 
—e,,/2. 

If the specimen is subjected to a time-de- 
pendent stress of the form of the real part of 


o =o re"! (3) 


of frequency w/2m, and the strain response is 
e..°e', the modulus of elasticity E is defined 
by the relation, 


1/E=e.,°/o0. (4) 


It is the object of the present theory to determine 
the real and imaginary parts of the elastic 
modulus E as functions of frequency on the ba- 
sis of the schematic model of polymer structure 
which has just been described. The elastic loss 
in the specimen is of course determined by the 
imaginary part of 1/E. If we write 


1/E=1/E'-i/E", (5) 


the rate of energy dissipation per unity volume 
is given by 


9 9 
y ae) 


abdijie——, (6) 


” 
us 





Since we wish to use the concept of relaxation 
time associated with the linear stimulus-response 
theory, we shall develop the theory in the linear 
approximation and will not consider departures 
from Hooke’s law. However, it should be re- 
marked that the theory will be valid for time- 





ind 


(1) 


the 
nal 
is 


(2) 


in 
rce 
ex 
lat 
= 4 
by 
he 

to 


ne 
tic 


ire 
SS 


he 


ne 


6) 


on 
se 
ar 
es 
e- 





ELASTIC LOSS AND RELAXATION TIMES 53 


dependent stresses of small amplitude super- 
imposed upon static stresses of arbitrary magni- 
tude. In the latter case the modulus E may, of 
course, depend upon the static stress in the 
initial state. It is sufficient to develop the theory 
for an harmonic time dependent stress of a single 
frequency. When the applied stress has arbitrary 
Fourier composition, each Fourier component 
can be treated independently by the methods to 
be described, insofar as the strain response €,, is 
linear in the stress. 

Each chain segment connecting neighboring 
lattice points is supposed to consist of » bonds 
r,°-'r,, each of length b. The vector R con- 
necting the terminal elements of the chain is 


> ry. (7) 


If we assume a fixed bond angle y, the internal 
configuration of the chain is most simply de- 
scribed by coordinates , -*:@n-1, where 
denotes the angle between the plane of bonds / 
and /+1 and the plane of bonds /+1 and /+2, 
in addition to the angles @ and ¢ specifying the 
orientation of any one bond, say the central 
bond of the chain, relative to the external 
coordinate system. For brevity in notation, we 
shall denote the set of internal coordinates ¢, 
-++@n_1, g, 8 by g, when it is only necessary to 
consider them collectively. 

The ensemble representing the statistical be- 
havior of the chains comprising a specimen of 
the polymeric material is characterized by a 
probability distribution function f(g, ¢t) in the 
configuration space g of a chain. If e; is a unit 
vector in the z direction, the projection of the 
average length of a chain on the z axis is 


(n= [Resa t)dq. (8) 


ln the equilibrium unstressed state (3), is equal 
to the lattice parameter a, given by 


a= f R-esfala)da, 
fo(q) =Ao exp (— Vo/kRT), (9) 
Ar= f exp (— Vo/kT )dq, 


where Vo is the potential of the average inter- 


molecular and intramolecular forces inhibiting 
internal rotational motion of the chain. 

As has been pointed out by James and Guth,’ 
the system of forces +F and —F of Eq. (3) 
applied to the terminal elements of the chain can 
be derived from the potential, — F(R-e,). Thus, 
according to our model, the potential of average 
force in the presence of an applied stress o in 
the z direction is 


V=Vo- (R e;) F, 
F=a/y, 


(10) 


where y is the number of chains per unit cross 
section. 

By the theory of Brownian motion, the distri- 
bution function f(q, t) satisfies the generalized 
diffusion equation! 


V-D-{Vf+fVV/kT} =df/dt, 
D=kTo™, 


(11) 


where V is the gradient operator in the g-space 
of a chain and D is the diffusion tensor, related 
by the second of Eqs. (11) to the resistance 
tensor ¢. It is sufficient for our purposes to 
consider applied stresses with harmonic depend- 


ence on the time ¢, 
¢= oe'”t, 


of frequency w/2x7. When the stress has an 
arbitrary Fourier composition, each Fourier 
component of the form (3) can be treated inde- 
pendently by the methods to be described, 
insofar as the strain response e,, is linear in the 
stress. 

If we expand the distribution function f(q, ¢) 
satisfying Eq. (11) in powers of the stress 
amplitude oo in the form, 


f=Ao exp (— Vo/kT) 
do 
x 14m )+0(0¥) | (12) 
V 


we obtain the following differential equation for 
the function f;(q, ¢): 


L(fi—(R-e1)/kT) —iwf, =0, 
w0; fi-[R-e,—a]/kT, (13) 
L=exp (Vo/kT)V-Lexp (— Vo/kT)D-V], 
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which has precisely the same form as in the 
theory of dielectric loss.! We expand the solution 
of Eq. (13) in the eigenfunctions y, a complete 
orthonormal set, of the derived operator Lo. 


Lor —Ap=0, 
Lo=V-D-V+u(q), 
u(q)=(V-D-VVo)/2kT 
—(VVo)-D-(VVo)/(2kT)?, 


(14) 


and calculate the average value of the z-extension 


of a chain, 
Az=(z)y—a, 


by Eqs. (8), (9), (12), and (13). The result is 





cof lal? 
Z= : et 
vkTL X 1+iw7 
n=1/2, (15) 


a=Aet | (R-e1) exp (— Vo/2kT)y*dq, 


where the sum extends over all eigenvalues X. 
Since the details of the derivation of Eq. (15) 
are identical with those of the theory of di- 
electric polarization,! we have omitted the rou- 
tine intermediate steps.* The mean square com- 
ponent (z”)y° of the chain length along the z axis 
in the unstressed state is related to the coeffi- 
cients 2, in the following manner, 


2 | 2 | 2 = (27), 


16 
(w= Ay [ (R-e) exp (= Vo/kT )dq. ( ) 


Thus, if we define a function Q) by the relation 
Qn = | 2y|2/(2?)w?, 
2 Aa=1, 


and employ Eq. (15) to obtain the elastic 


* Properly, in order to provide for the normalization of 
fi(g, t), the 2, of Eq. (15) should be replaced by (Az), where 


(As), =2,—aAot f exp (— Vo/2kT yr*dg. 


However, for the small equilibrium extensions a in the 
unstressed state, the approximation z, is adequate and is 
introduced here largely for simplicity in the subsequent 
discussion. 


modulus £, defined as oo/e,,°, we obtain 
1 (2?) w° Qy 


—=—— 5 —. (17) 
E MkT ®* 1+27y 





If the discrete relaxation time spectrum 7) is 
closely spaced relative to 1/w, it may be approxi- 
mated by a continuous spectrum with the 
result,* 








1 1 f% Xr) 
—=— dr, 
E Eo 0 1+iwr 
dQ 
O(r)=-——; N= DA, (18) 
> A<1/r 
1 (37) py 
E) MkT’ 


where Ey is the static modulus of elasticity. 
Before analyzing the relaxation time distribution 
function ®(7r), a few remarks concerning the 
static modulus are desirable. If the average force 
potential V» is assumed constant, the free rota- 
tion approximation, the average value (z*),° has 
the well-known form, 


(2°) 0° = (R?)w°/3 =nl?/3, 


(19) 
l=b tan y/2, 

* We remark here that ®(7)dr represents the fractional 
elastic compliance in tension associated with relaxation 
times in the interval between 7 and r+dr. For an incom- 
pressible material with a Poisson ratio of 4, the shear 
modulus G is equal to E/3. Thus, if one wishes to speak of 
elastic compliance in shear, the function J(r) equal to 
36(7r)/Eo may be used. The elastic compliance in shear 
associated with relaxation times between 7 and +r+dr is 
then equal to J(r)dr. 

The theory which we have developed leads directly to a 
distribution of the elastic compliance, not the elastic 
modulus, in the relaxation time spectrum. The distribution 
is therefore not Maxwellian. In order to bridge the gap 
between the viscous fluid and the elastic solid, Maxwell 
proposed a relaxation mechanism, which for a single 
relaxation time is described by the equation, 


do .¢a de 
as a: 


Under a constant stress the material creeps. For a stress 
with harmonic dependence on time, the elastic modulus E is 


E=wwrtE./(1+iwr). 


Thus, if the frequency is large relative to 1/7, E approaches 
E., and the stress approaches proportionality with the 
strain. When the frequency is small relative to 1/7, the 
stress approaches proportionality with the rate of strain. 

Although a part of the elastic loss of cross-linked 
polymers, in which creep is observed, is no doubt due to 
mechanisms formally described by the Maxwellian scheme, 
such contributions are not included in the present theory. 
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where 7 is the number of bonds comprising the 
chain segment, b the bond length, and y the 
bond angle. The same result is obtained if each 
bond is an axis of symmetry or contains a plane 
of symmetry of the intramolecular torque hinder- 
ing rotation around that bond. The influence of 
the average intermolecular hindering torque on 
2*)w° is much more difficult to determine. If we 
make the assumption of James and Guth* that 
the average intermolecular force hindering the 
internal motion of the chain can be represented 
by an internal stress o; distributed over the 
lattice points of the net, with 


Vo= —(R-e1)0;/», 
for a chain oriented in the z direction, we obtain 


(3?)w° =nl?/9, 
(20) 
l=b tan y/2, 


a result which differs from that of Eq. (19) by 
the factor 4. The elastic modulus Ep» in the 
unstressed state then agrees with that of James 
and Guth.’ Undoubtedly, an entirely satisfactory 
quantitative theory of the elastic modulus Eo 
depends upon an exact analysis of the potential 
Vo of the intramolecular and intermolecular 
hindering torques acting on a chain segment of 
the network. 

In investigating the relaxation time distribu- 
tion, (7), we must first determine the internal 
rotatory diffusion tensor D entering into Eq. 
(11). Fuoss and the writer! have estimated D 
for the most probable configuration of a freely 
rotating polymer chain. It was assumed that 
each monomer unit, moving subject to the re- 
straints imposed upon it by the other elements 
of the chain, was subject to a resistance propor- 
tional to its velocity. The friction constant £o of 
éach such monomer unit was further supposed 
to be independent of the chain length , but of 
course dependent upon the environment of the 
chain. The existence of the friction constant £o 
implies merely the validity of the theory of 
Brownian motion. Although Fuoss and the 
writer suggested the use of the usual hydro- 
dynamic estimate of {9 by Stokes’ law, none of 
the general results of the theory of dielectric 
loss nor those of the present theory depend upon 


this estimate, as implied by Kauzmann.® If one 
desires, one may use the transition state theory 
for estimating {0, as advocated by Kauzmann 
and Eyring. The observation of Kauzmann and 
Eyring,® that the activation energy for viscous 
flow in fluids composed of long chain molecules 
becomes independent of chain length for n of the 
order of 20, lends support to our assumption 
that ¢) can be treated as independent of n. 

The estimate of the diffusion tensor D of Fuoss 
and the writer! leads to the following expression 
for the operator. V-D-V of Eq. (11), 


V-D-V= 








Dor 1 #9 rs) 
| —f{ sin )— 
2n*Lsin & A3 Ov 





sin? 3 d¢? 
4D, 
D,= , 
[n—|2s—n|]? 
Do=3kT/b7%%o, 





1 81,2 2 
|#z,2. 
atin Og," (21) 





where 6 is the bond length. We remark that D, 
is of the order Dy for bonds near the chain ends 
and of the order of D)/n? for bonds near the 
middle of the chain. This result appears reason- 
able, if we remember that D, measures the 
resisting torque experienced by the molecule 
when rotation occurs around bond s alone, the 
configuration of the other bonds remaining fixed. 

Use of Eq. (21) in Eq. (11) and the approxi- 
mation of free rotation lead to the following 
relaxation time distribution function for elastic 
loss, 


Tm 


$(1)=——, 
(r+ Tn)* 
(22) 


Tm = NTO, 
To= 2b7¢0/kT. 


The details of the determination of G(r) are 
formally identical with those of the theory of 
dielectric loss.1 It should be emphasized, how- 
ever, that the theory of elastic loss applies only 


6W. Kauzmann, Rev. Mod. Phys. 14, 12 (1942); W. 
Kauzmann and H. Eyring, J. Am. Chem. Soc. 62, 3113 
(1940). 








56 JOHN G. KIRKWOOD 


to cross-linked polymers. A polar cross-linked 
polymer should exhibit the same relaxation time 
spectrum in elastic and dielectric loss. This will 
not in general be true of a thermoplastic polymer 
such as polyvinyl chloride to which the model 
underlying the theory of elastic loss here de- 
scribed does not apply. It is possible that 
temporary cross-linkages due to intermolecular 
forces can be regarded as existing in thermo- 
plastic polymers. On this assumption, elastic 
loss might consist of a contribution of the type 
described here superposed upon creep due to 
slippage at the temporary cross links. 

The frequency of maximum loss for the 
distribution (22) is given by 


@m=1/n7O9, (23) 


and is thus inversely proportional to the length 
n of the chain segments between cross-links. In 
a rubber, the use of Stokes’ law with an effective 
viscosity coefficient for estimating 9 and 79 can 
scarcely be defended, although such an estimate 
is adequate in order of magnitude for a polymer 
chain dispersed in a low molecular weight sol- 
vent. In the absence of an adequate theoretical 


method for estimating 79, we remark that 


ro =2b?/D®, (24) 


where D° is the translational diffusion constant 
of a monomer unit of the chain. We therefore 
suggest that reasonable, although perhaps rough, 
estimates of 7) might be obtained from experi- 
mental diffusion constants in the rubber, of small 
molecules of the approximate size and shape of 
the representative monomer unit. For example 
in a Buna S rubber, one might measure the 
diffusion constant of benzene or toluene in the 
rubber itself for use in Eq. (24) to estimate 7» 
for the styrene units of the polymer, to which 
most of the loss could probably be attributed. 
In conclusion we remark that the relaxation 
time distribution, Eq. (22), is to be regarded 
only as a first rough approximation. Properly, 
the influence ‘of the intermolecular hindering 
torques should be taken into account in order to 
improve the free rotation approximation. A 
possible method of accomplishing this could be 
based upon the internal stress concept of James 
and Guth. An investigation along these lines is 
in progress and will be reported at a later time. 
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The Mercury Photosensitized Reactions of Propylene and Isoprene* 


H. E. GuNNING AND E. W. R. STEACIE 
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(Received November 21, 1945) 


An investigation has been made of the reactions of propylene and isoprene with mercury 
(3P;) atoms, at 30°C, in a static system. The reactions are similar in that the rates of con- 
sumption of both propylene and isoprene decrease with increasing pressure in the range of 
complete quenching. In addition, the quantum yields of both reactions are low. The products 
of the propylene reaction are liquid polymer and hydrogen, together with smaller amounts 
of methane, ethylene, and hexene. The principal steps in the mechanism proposed for the 
propylene reaction are C;Hg+Hg(#P1)—>C;He*+He('50); CsHe*+CsHe>2C3He; CsHs*— 
C3;H;+H; H+C;Hs>C3;H;+H:; H+C;Hs~>CH;+C:H;; followed by polymerization reac- 


tions of the type C;H;+C3;He—CeHu,, etc. 





INTRODUCTION 


HE saturated and the unsaturated hydro- 
carbons exhibit striking differences in their 
reactions with mercury (*P;) atoms. In the 
paraffins the initial step seems to be of the type! 


CrHoen+2 + Hg(?P,)-C,, He, +1 + H + Hg(}So) (1a) 


or 


— CrHongitHeH. (1b) 


The investigations of LeRoy and Steacie on 
ethylene,? and of Gunning and Steacie on buta- 
diene,’ at room temperature, suggest that the 
initial step in the reactions of these molecules 
with mercury (*P;) atoms involves the formation 
of an activated molecule, or perhaps an acti- 
vated hydrocarbon-mercury complex. LeRoy and 
Steacie,? for example, postulated the following 
sequence for ethylene: 


C.H4+Hg(@P1)-C2H,* +Hg('Sp), (2) 


C,H4*+C.Hy2C2Ha, (3) 
C,;H,*—-C2H2+Ha, (4) 
H.+He(@P;) 2H +He(!S>). (5) 


Polymerization was initiated by the H atoms 
produced in reaction (5). 
One obvious structural difference between the 


* Contribution No. 1345 from the National Research 
Laboratories, Ottawa, Canada. 

1See D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 
10, 676, 683 (1942) for a summary. 

2D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 9, 
829 (1941). 

3H. E. Gunning and E. W. R. Steacie, J. Chem. Phys. 12, 
484 (1944). , 
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paraffins, on the one hand, and ethylene and 
butadiene, on the other, is that in the latter all 
the C—H bonds are olefinic. In this connection, 
Taylor and his associates‘ found a high activation 
energy (>10 kcal.) for the reaction of methyl 
radicals with ethylene, while for the correspond- 
ing reaction with propylene they found an 
activation energy of only 3.1 kcal. They at- 
tributed the low activation energy of methane 
formation with propylene to the weakness of 
C—H bonds in a position once removed from a 
double bond. Thus, in propylene, they deduce 
that the C—H bond strength in the methyl 
group must be about 5 kcal. less than in ethane; 
or in other words about the same value as the 
secondary C—H bond in propane. 

With the above considerations in mind, an 
investigation of the reactions of mercury (*P;) 
atoms with propylene and isoprene, both pos- 
sessing C—H bonds once removed from a C=C 
bond, seemed particularly instructive. The de- 
tails of these investigations are reported below. 


EXPERIMENTAL 


For a description of the apparatus employed, 
and the experimental technique, reference should 
be made to the paper of Gunning and Steacie on 
the mercury photosensitized reactions of buta- 
diene.* No important alterations were made in 
the apparatus for the present investigations. 

The propylene was obtained from the Ohio 
Chemical and Manufacturing Company. After 
several trap-to-trap distillations, the middle 


4H.S. Taylor, Ann. N. Y. Acad. Sci. 41, 231 (1941). 
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Fic. 1. Pressure vs. time for a typical static run. 


fraction was collected and stored in a flask, 
provided with a Warrick-Fugassi valve.® 

The isoprene was an Eastman product. It was 
dried over freshly fused calcium chloride, and 
stored in a small bulb, attached to the cell 
system by a Warrick-Fugassi valve.® 

All runs were made in a static system. The 
cylindrical quartz cell employed had a volume of 
240 cc, and was immersed in a water thermostat, 
maintained at 30.00+0.01°C. 


RESULTS 


I—Propylene 
The Effect of Pressure on the Rate 


Figure 1 shows a typical one-hour run at an 
initial pressure of about 1 mm of propylene. As 
soon as the lamp is turned on, the pressure 
begins to fall linearly with time, and at the same 
time the window of the cell slowly becomes 
covered with tiny droplets of liquid polymer. 
When the amount of polymer accumulated on 
the cell window becomes sufficient to decrease 
effectively the amount of \2537 entering the cell, 
the rate of pressure decrease begins to fall off. 
Unlike butadiene or ethylene, there is neither a 
period of induction nor an initial pressure rise 
observable. The polymer can be seen on the 
windows of the cell even in the early stages of 
the run. 


§ E. Warrick and P. Fugassi, Ind. Eng. Chem. Anal. Ed. 
15, 13 (1943). 


GUNNING AND E. W. R. 


STEACIE 


In Fig. 2, the actual rate of propylene con- 
sumption, [ —d(C3H¢) ]/dt, in moles per minute, 
is plotted against the initial pressure. The rate 
of propylene consumption, as can be seen from 
the diagram, reaches a maximum at about 10 
mm. Below 10 mm, the rate is less, owing un- 
doubtedly to incomplete quenching. The decrease 
in rate above 10 mm strongly suggests collisional 
deactivation. 


Products of the Reaction 


The products of the reaction were particularly 
difficult to investigate, since over 90 percent of 
the propylene molecules decomposed went to- 
ward the formation of polymer. The procedure 
adopted was to analyze the gas, non-condensable 
in liquid air, and the residual propylene for each 
run. A complete analysis was done by collecting 
the gas from 10 to 15 runs at 10-mm initial 
pressure, and performing a careful fractionation 
of the products using the modified Ward method.* 

The products, non-condensable in liquid air, 
proved to be a mixture of hydrogen and smaller 
quantities of methane. The results of these 
analyses are tabulated in Table I. Each value 
is the average of two separate runs. 

The identity of the methane was established 
by combustion in excess oxygen. 

Analysis of the accumulated products resulted 
in the isolation of a small amount of C. fraction, 
which upon combustion proved to be ethylene. 
At 10-mm pressure, approximately 0.01 mole 
of ethylene was produced per mole of propylene 
decomposed. 

Since it was suspected that methyl acetylene 
might be present, a careful analysis of the Cs; 
fraction was made for acetylenic compounds, 
using the method of Ross and Trumbull, as 
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Fic. 2. The rate of propylene consumption in moles per 
minutes vs. the initial pressure. 
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modified by LeRoy and Steacie.? No trace of 
any such compound could be found. Allene, 
however, which is an isomer of methyl acetylene, 
might be present. 

The fractionation of the products of 10 runs at 
10 mm produced about 0.1 cc of a volatile liquid, 
whose vapor pressure corresponded roughly to 
hexene. A solution of this liquid in chloroform 
was found to discolor a dilute solution (0.05M) 
of bromide in glacial acetic acid, suggesting the 
presence of one or more double bonds in the 
compound. 

The polymer formed in the reaction was a 
clear, colorless, viscous liquid, completely soluble 
in benzene. The fact that the polymer is both 
liquid and soluble would suggest that the mole- 
cules are of the linear type, and of relatively 
low molecular weight. 


Quantum Yield 


A relative determination of the quantum yield 
was made by comparison with the initial pressure 
rise in the mercury photosensitized reaction of 
ethylene. LeRoy and Steacie® obtained the value 
0.37 at 13 mm for the quantum yield of the 
ethylene reaction. The average of four runs gave 
the value 0.15++0.05 at 14 mm for the propylene 
reaction. 


II—Isoprene 


Eight runs were made with isoprene over the 
pressure range from 40 to 225 mm. The rate of 
pressure decrease vs. initial pressure is plotted in 
Fig. 3. The rate falls off in a similar manner to 
propylene, but not quite as rapidly. 

A complete analysis of the products has not 
yet been made. However, preliminary analyses 
show that well over 90 percent of the isoprene 
decomposed forms polymer. A small amount of 











TABLE I. 
A B C 

Po —d(C3Hs) /dt d(H2)/dt d(CHa)/dt B/C 
mm Moles per min. X 107 

1.06 5.33 0.462 

1.91 7.63 0.51, 0.085 S.J 

4.23 8.63 0.535 

9.62 9.75 0.564 0.13, 4.2 
24.9 8.24 0.405 “0.102 4.0 
47.3 6.17 0.22, 0.035 6.5 
70.0 4.9; 0.127 0.036 39 








gas, non-condensable in liquid air, was formed, 
which upon analysis proved to be pure hydrogen. 
At 40 mm initial pressure, 0.0014 mole of 
hydrogen was formed per mole of isoprene 
decomposed, in a one-hour run. 
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Fic. 3. The rate of pressure decrease in mm per minute vs. 
the initial pressure. 


Since the same light intensity was used in 
both the isoprene and the propylene runs, the 
quantum yield for the isoprene reaction must 
be in the neighborhood of 0.2. 


DISCUSSION 
I—Propylene 


Any mechanism proposed to explain the reac- 
tion of propylene with mercury (*P;) atoms must 
be consistent with the following facts. 

(a) The pressure-time curve is linear even in 
the initial stages of the reaction. The rate of 
pressure decrease diminishes owing simply to 
accumulation of polymer on the window of the 
cell. 

(b) Polymer is present even in the initial 
stages of the reaction. 

(c) The products of the reaction are a liquid 
polymer, soluble in benzene and hydrogen, to- 
gether with smaller amounts of methane, ethy!l- 
ene, and possibly hexene. 

(d) The rate of propylene consumption, and 
the rates of hydrogen and methane formation 
decrease with increasing pressures above 10 mm. 

(e) The quantum yield is 0.15. 

The decrease in rate with increasing pressure 
over the range of complete quenching suggests, 
by analogy with ethylene and butadiene, a 
primary step of the type: 


C3H.+Hg(?P1)-C3He* +Heg(!Sp), (6) 


where C;H,* represents either an activated 
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molecule, or perhaps an activated propylene- 
mercury complex. Both the decreasing rate and 
the low quantum yield could then be satis- 
factorily accounted for by the collisional deacti- 
vation reaction 


C3H,*+C3H—5—2C3H.¢. (7) 


The absence of either an initial pressure rise, 
or a period of induction in the pressure-time 
curve, together with the fact that polymer is 
present even in the early stages of the reaction, 
suggests that free radical processes must occur 
early in the reaction. The reaction 


C;H.e*—C3H5;+H (8) 


would have the net effect of generating two 
radical chains for each mole of propylene decom- 
posed. Indeed in view of the presence of the 
weak C—H bonds in propylene,‘ which differ- 
entiate propylene from ethylene and butadiene, 
reaction (8) might be expected to be faster than 
the corresponding reaction for butadiene and 
ethylene. 

The H atoms generated in (8) would react 
rapidly with propylene generating propyl radicals 


C3;H,+H—-C3H7:. (9) 


The activation energy of this reaction is low, 
probably less than 5 kcal.® 

As Rabinovitch, Davis, and Winkler’ pointed 
out; reaction (9) is exothermic to the extent of 
about 40 kcal. Consequently the hydrogen, 
methane, and ethylene found in the products 
could be explained by the reactions 


C;H.+H—C;Hs+H;, (10) 
C3;Hs+H—-C.H.+CHs, (11) 
CH3;+C3Hs—-CH,4+C3Hs, (12) 


where (10) and (11) involve an energy-rich 
‘propyl radical as intermediate. Reaction (12) 
has an activation energy of 3.1 kcal.,4 and it is 
likely that reaction (10), which is the analogous 
reaction with H atoms, would have approxi- 
mately the same activation energy. In addition 
to forming methane, the methyl radicals would 


6 Annual Tables of Physical Constants, Section 602(c) 
(1942). 

7B. S. Rabinovitch, S. G. Davis, and C. A. Winkler, 
Can. J. Research B21, 251 (1943). 
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generate polymer chains by the reaction 


CH3+C3H.—CiHsg. (13) 


By analogy with the reaction of methyl radicals 
with ethylene,® this reaction should have an 
activation energy of about 1.5 kcal. 

From Table I, it can be seen that the rates of 
both the hydrogen and the methane formation 
decrease with increasing pressure over the range 
of complete quenching. Moreover, the ratio of 
these two rates is a constant within experimental 
error, and equal to about 4.8. The sequence (9) 
to (13) satisfactorily accounts for these facts 
We can express this ratio as a function of the 
individual rate constants by a _ steady-state 
calculation, from which we obtain the expression 


d(H:) /dt kif Ris 
=4.8= |! +=| 
Ruy Rie 


d(CH,) /dt 

If we assume that all these reactions have the 
same steric factor, and assign the values E1)9=3.1, 
Ey=3.1, and E,3;=1.5, to reactions (10), (12), 
and (13), respectively, reaction (11) would have 
the surprisingly low value of about 2-3 kcal. for 
its activation energy. Of course this reaction 
can occur only through the formation of an 
activated propyl radical. 

Polymerization would be initiated in the main 
by the allyl radicals formed in (8), (10), and (12), 


C3H;+CsHs-CeHu, etc... (14) 


although, of course, the propyl radicals formed 
in (9), and the methyl radicals formed in (11), 
are both capable of forming polymer chains. The 
growing polymer chains formed by a continua- 
tion of reaction (14) could be represented by the 
formula (C3H¢)n°C3Hs5. It is not likely that x 
would be very large, since as Taylor and Smith® 
pointed out, propylene is an extremely efficient 
chain breaker, and, therefore, chain-transfer 
steps of the type 


(C3He) n° C3H5+CsHs—(CsHs) n41 +C3Hs (15) 


would be expected to occur before n became 
very large. Reaction (15), it will be observed, 
would have the net effect of breaking one chain 
and initiating another, and consequently the 


8H. S. Taylor and J. O. Smith, Jr., J. Chem. Phys. 
8, 543 (1940). 
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short-chain liquid polymer isolated in the 
mercury-photosensitized reaction could be satis- 
factorily accounted for by reactions of this type. 

The small amount of liquid dimer found in 
the product might have been formed either by 
the dimerization of propylene, or by the recombi- 
nation of allyl radicals. 


II—Isoprene 


Both the low quantum yield and the decreasing 
rate with increasing pressure seem to suggest 
that isoprene follows the same general pattern of 
reaction as propylene, 


CsHs+Heg(?@P1)-CsHs*+Heg('So), (16) 
CsHs*+Cs;Hs—2C;Hs, (17) 
C;H,*->C,;H;+H, (18) 


followed by polymerization by the free radicals 
formed in (18), 


C5Hs+CsH7— CoH 15, etc. (19) 


The very small quantity of hydrogen formed 
would indicate that, if it is produced by a similar 
reaction to propylene, i.e., 


CsHs+H—-C;H;+Hz, (20) 


this reaction must occur only to a very small 
extent. 

Since a detailed analysis of the products has 
not yet been made, it is not possible at this time 
to set up a complete mechanism for the reaction. 
However, it is fairly clear from the preliminary 
investigation that isoprene reacts with mercury 
(®P,) atoms in an essentially similar manner to 
propylene. 


CONCLUSIONS 


The investigations of LeRoy and Steacie on 
ethylene,? of Gunning and Steacie on butadiene,’ 
and the present investigation on propylene and 
isoprene seem to show that the substituted and 
the unsubstituted alkenes and alkadienes all 


react with mercury (*P;) atoms at room temper- 
ature by an activated molecule mechanism. The 
essential difference between ethylene and buta- 
diene, on the one hand, and propylene and 
isoprene, on the other, lies presumably in the 
presence in the latter of weakly bound hydro- 
gens in their methyl groups. 

If we allow C,H, to represent any of the 
above-mentioned hydrocarbons, we find that the 
two initial steps, 


C,H, +Hg(?P1)-C.H,*+Hg(So) = (21) 
and 


C,H,*+C,H,—2C.H,, (22) 


are common to all. 

The activated molecules formed from ethylene 
and butadiene split mainly into hydrogen and a 
more unsaturated hydrocarbon, 


C,H,*—-C,H y_») + He. (23) 


The concentration of hydrogen builds up in the 
system until it quenches the mercury (*P1) atoms 
to an appreciable extent by reaction (5), thereby 
initiating polymerization. 

The activated propylene and isoprene mole- 
cules, on the other hand, owing to their weak 
C—H bonds, split into an H atom and a radical, 


C,H,*--C,Hy_» +H, (24) 


initiating polymerization immediately thereby. 

The low value for the quantum yield (0.15-—1.0) 
found in all these reactions suggests that reaction 
(22) must be extremely efficient. However, there 
are probably other factors involved. Under static 
conditions the concentration of mercury atoms 
in the cell is always high enough to absorb all 
the resonance radiation within a few millimeters 
from the cell window. This fact, coupled with 
the high quenching efficiencies of the unsaturated 
hydrocarbons, would seem to limit the zone of 
reaction to a thin layer close to the cell window, 
rather than throughout the body of the cell. 
Under these conditions it is likely that wall 
effects would become significant. 
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Based on classical models including the interpretation of the deBroglie wave as an electro- 
magnetic ‘‘atmosphere’’ which controls the orbit and the ‘‘action”’ of the classical electron, 
and postulating that the ‘‘atmosphere”’ of the electron suffers a permanent change only when 
the electron undergoes a change in its quantum state, an equation is derived for the He ground 
state iso-electronic sequence which agrees substantially with experiment. This view leads to 
the conclusion, moreover, that the binding energy in molecules formed without electronic 
transitions is caused entirely by a decrease in the average potential energy, the average kinetic 
energy of the molecule remaining the same as in the constituent atoms. This conclusion is 
verified in the systems H2 and H2*. Introducing the additional assumption that the resonance 
contribution to the total energy of a state, arising from an ‘‘adiabatic’’ orbital distortion, is 
proportional to the orbital frequency, equations are derived expressing the term values of a 
number of states of two-, three-, and four-electron atoms. The methods appear to apply quite 


generally. 





N a recent article’ it was shown that the 

molecules He and H,+ may be treated with 
good accuracy by classical models and old 
quantum theory employing in addition a concept 
of resonance. The purpose of this discussion is to 
apply similar methods in the calculation of 
energy states (ground and excited) of several 
many-electron atoms. In dealing with classical 
models in light of present knowledge one quite 
naturally is led to adopt an interpretation of the 
wave characteristics of matter similar to the 
views held by Thomson,? Schott,* Japolsky* 
and others, namely, that the wave associated 
with the electron is simply an electromagnetic 
“atmosphere.”’ The chief concern here is to 
justify the stationary states, i.e., the non-radia- 
tion of the field of the accelerating charged 
particle. Schott proved with Maxwell’s equations 
that non-radiating orbits of charged particles are 
indeed possible although he was not able to 
account for the Bohr orbits of hydrogen. One may 
wonder perhaps what part the nucleus plays in 
the non-radiation of the field of the electron in its 
stationary orbits.’ It seems not altogether im- 
possible, when we realize that the deBroglie 


1M. A. Cook, J. Chem. Phys. 13, 262 (1945). 

2J. J. Thomson, Phil. Mag. [7] 5, 191 (1928); 6, 254 
(1928); 8, 1073 (1929); 9, 1185 (1930); 12, 1057 (1931); 16, 
809 (1933); 27, 1 (1939). 

3G. A.Schott, Proc. Roy. Soc. (London) A139, 37 (1933) ; 
A171, 1 (1939). 

‘4N. Japolsky, Phil. Mag. 20, 417, 641 (1935); 22, 1003 
(1936). 

5 W. M. Hicks, Phil. Mag. [7] 8, 108 (1929). 
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wave-length of the proton would be the same as 
that of the electron in the classical model of 
hydrogen, that the Poynting vector may be 
nullified by an over-lapping of opposite fields 
when the boundary conditions are satisfied. At 
any rate there does not seem to be any very good 
reason in electromagnetic theory for discounting 
the possibility of standing light waves sur- 
rounding the elementary charged particle. On the 
other hand, this postulate certainly allows a 
clearer comprehension of the control of the orbit 
by conservation of energy, and even of the 
underlying meaning of quantum states. We 
visualize the ‘“‘atmosphere’”’ as not only providing 
the control of the orbit but of being the source of 
the ‘‘action,” i.e., the momentum of the particle. 
In this view there is no mysterious condensation 
of field to ponderable matter, but the light quan- 
tum, when it is absorbed by matter, merely super- 
imposes itself on the ‘atmosphere’ already 
present. However these matters, while of philo- 
sophical interest, are important only when in 
serving as basic postulates they allow us to make 
further progress in the interpretation of the laws 
of nature. 


KINETIC VS. POTENTIAL ENERGY IN SYSTEMS 
WITH RESONANCE 


According to the virial theorem the total en- 
ergy of systems operating under inverse square 
forces alone is equal to the negative of the aver- 
age kinetic energy (7), or to half the average 
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potential energy (V). However, as shown by 
Slater,® the vibrating nuclei in molecules have the 
effect of applying a (non-Coulombic) constraint 
at any nuclear separation other than the (static) 
equilibrium one. Moreover, it was shown! that 
the true (dynamical) equilibrium distance in Hy» 
and H,* differs widely from the static equi- 
librium position. Thus, in the circular model of 
He with rigid nuclei, R (the internuclear dis- 
tance) turns out to be .58A, but the true value 
is .J4A. As a consequence it is found that 
T# —E#X —1/2V. Slater pointed out in fact that 
the binding energy in molecules seemed to be 
caused largely (it appears from our analysis, 
entirely) by a decrease in V. This is of course 
contrary to the considerations of Hirschfelder 
and Kincaid’ who evaluated the various treat- 
ments of H,2 on the basis that a correct solution 
should give T= —1/2V. 

Now the wave-length of the particle, according 
to deBroglie, is a function of the classical mo- 
mentum, and a change in momentum must 
accordingly be accompanied by a change in the 
wave-length. By our interpretation of the 
deBroglie wave as a peculiar type of electro- 
magnetic ‘“‘atmosphere,” this seems to imply an 
absorption or emission of radiation whenever T 
undergoes a change. However, no such radiation 
occurs in the formation of most molecules, e.g., 
H, from two H atoms; it is necessary for the H 
atoms to lose energy to a suitable third body 
before they will form the molecule. We suspect, 
therefore, that there may be no change in T in 
the “adiabatic” formation of the molecule from 
its atoms. We recognize that the foundation of 
this postulate may seem to conflict with the fact 
that the binding energy will to a large extent 
appear as kinetic energy (translational) in other 
systems. It is true also that if resonance exists 
between the various particles of the system there 
must be a corresponding resonance interchange 
of field between the particles as required by the 
deBroglie formula. However, let us re-examine 
the Hz and H,* solutions of reference (1) to see if 
it may be verified that the average kinetic energy 
remains constant in molecular formation. 

From Eqs. (2) and (3) of reference (1), it will 


6 J. C. Slater, J. Chem. Phys. 1, 687 (1933). 
agg Hirschfelder and J. Kincaid, Phys. Rev. 52, 658 


be seen that RF=1/2RF,, since 7+1/2V 
=-—1/2RF=-—1/4RF;,. This is apparent also 
from the fact that RF,;=0 as seen by Eq. (3b). 
Thus RF =1/2ge?/# at the equilibrium distance R, 
where 7 is the root mean square average of r; and 
ro. Since T=—E-—RF, we thus have T=—E 
—1/2ge?/7. Employing the values of 7; and 72 in 
reference (1) corresponding to R=1.06A and 
.74A, we obtain for T the values 13.39 and 
26.77 ev for H2*+ and He, respectively. (In the 
case of Hz it was assumed as in reference (1) that 
6= at R=.74A.) These values differ from the 
expected ones by almost exactly the vibrational 
energy which would be entirely kinetic at the 
equilibrium position. Hence, as expected, the 
average kinetic energy of these molecules is at 
least approximately the same as in the corre- 
sponding constituent atoms. 

If as a general rule, then, we postulate that T 
remains unchanged in the ‘‘adiabatic”’ interaction 
of atoms to form molecules so far’as the internal 
energy is concerned, we conclude that the binding 
energy in molecules formed without change in the 
electronic configuration is caused entirely by a 
decrease in V. This is an interesting consideration 
if true. In fact, if it can be justified, it would 
provide a very valuable principle in the study of 
molecules, the interaction of light and matter, 
and in other problems. It thus becomes of interest 
to extend the postulates to other resonance 
systems. In a restricted sense the same sort of 
thing is involved in the structure of many 
electron atoms as in molecules. That is, as a 
second electron enters an atomic system, the 
orbital of the first one must suffer a distortion 
caused by the electron-electron interaction. This 
would be “adiabatic” in the same sense as in 
molecules since it involves no change in the 
quantum numbers of the first electron. 


TWO-ELECTRON ATOMS (GROUND STATES) 


We shall adopt a modified circular model for 
this system similar to that employed in Hg. It is 
evident, however, that the model must be some- 
what different in this case as the resonance 
phenomenon must be different in character from 
that found in molecules, the latter being asso- 
ciated to a large extent with nuclear vibrations. 
If there were no resonance or exchange energy 
involved, the ground state energy of two-electron 
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TABLE I. Ground states of helium and helium-like 











atoms (ev). 

Atom —E (k =.725) —E (obs.) 
H- 14.29 

He (78.61) 78.61 
Lit 197.07 197.10 
Be?* 369.70 369.75 
Bs+ 596.46 596.59 
C+ 877.38 877.67 








atoms according to the circular model would be 
given by the equation 


E' = —(4Z—1)@/4r = —(4Z—1)*Rhc/8, (1) 


with r=4a9/(4Z—1) for the coincident orbital 
radius, the other quantities having their usual 
significance. We consider the circular model to be 
merely transitory in the dynamic structure of the 
atom, representing an equilibrium structure as in 
H». The energy in (1) is, therefore, dashed to 
show that it is not the true value. According to 
the above postulate at the configuration where 
the circular model applies there should be an 
additional energy (the resonance energy) to 
account for the ‘‘adiabatic”’ change or distortion 
of the original one-electron orbit. That is, we 
consider that radiation takes place only from the 
field of the second electron entering the system, 
and that there should thus be an additional 
kinetic energy given by 


T=} f ‘(Z—be/rdr, (2) 


where ¢ and 79 are the radii of the coincident two- 
electron and the original one-electron orbits, re- 
spectively, and & is a factor accounting for the 
electron-electron interaction. The factor } ac- 
counts for the fact that we are dealing with 
kinetic energy, the integral being of the form of a 
potential energy term. We shall leave k unde- 
termined, noting that it is of the same form as Z, 
and therefore is the same for all members of the 
sequence. Hence the arbitrary assignment of any 
value for k should give values for E (= E’+7") 
for the normal He iso-electronic sequence which 
differ from the true values by a constant amount 
plus the spin and relativity corrections, the 
integral (2) turning out to be a constant (Rhc/4) 
times Z—k. Hence, 


E=E'+T’=[2(Z—k) —(4Z—1)"]Rhc/8. (3) 


Table I shows values of E obtained from (3) for 
k=.725, the value needed to give agreement with 
experiment in normal He. It is seen that the 
agreement is about what should be expected with 
spin and relativity contributions neglected. This 
treatment does not of course give the detailed 
mechanism of the resonance effect although one 
might guess that it involves orbital pulsations 
and/or oscillations through the coincident circu- 
lar model orbit as an equilibrium position. 


EXCITED STATES IN MANY-ELECTRON ATOMS 


The nature of the solution for normal helium 
suggests that the energy of other configurations 
in many-electron atoms may also be expressed in 
terms of (a) an ideal H-like or He-like term or 
terms, the He-like term of the form shown in 
(1). applying with coupled .electrons, plus (b) a 
term or terms due to ‘‘adiabatic”’ orbital distor- 
tion for electrons not undergoing radiation in 
reaching the state. In considering excited 1S 
states of two-electron atoms the ideal term should 
be an H-like term with effective atomic number 
Z'=Z-—1, i.e., with perfect screening. We shall 
designate it H, for the state of quantum number 
n. Now let us postulate on the basis of classical 
models that the average distortion of the orbit in 
excited states is proportional to the orbital fre- 
quency which seems to be a logical assumption. 
This suggests also that since the ground state is 
really a! S state, the orbital distortion or resonance 
term should be proportional to (Z—k)Rhc/4, 
although it is apparent that there is an abrupt 
change in the characteristics of the orbit due to 
breaking of the electron coupling. Let us suppose 
for the moment, however, that this merely intro- 
duces a constant factor. Since the frequency 
factor should be given to a close approximation at 
least by E,/nE,, where E, is the term value 
taken as positive by convention and £, is the 
ionization potential, the following equation 
should hold according to our postulates for the 1S 


TABLE II. Term values of 4S series in He and Li* (ev). 








He Li* 





State calc. obs. calc. obs. 
2’s 3.94 3.95 14.74 14.80 
3's 1.66 1.66 6.35 6.33 
4's 91 91 3.53 3.51 
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TABLE III. Term values of various He-like series from Eq. (4) (ev). 











He Li* Be?* B3+ c+ 
Term calc. obs. calc. obs. cale. obs. calc. obs. calc. obs. 
2'p (3.35)  <3.d0 13.34 13.36 30.06 30.07 53.54 53.53 (83.79) 83.79 
3'p 1.49 1.49 5.96 5.96 13.42 13.43 23.89 23.95 37.37 37.39 
4’p .84 .84 3.36 3.36 1.59 7.57 13.46 13.47 21.05 21.06 
(ki = —.52; k=1.705) 

23s (4.74) 4.74 (16.53) 16.53 

38s 1.86 1.86 6.84 6.82 

43s .99 .99 Re 3.41 

(ki =3.92; k=.94) 

2p (3.61) 3.61 (14.29) 14.29 

3%p 1.57 157 6.24 6.24 

43p .88 .88 3.48 3.48 


(ki =1.44; k=1.36) 








TABLE IV. Ground and excited states in Li and Li-like atoms (ev). 














Li Be* B+ C3+ Nét O5+ F6+ 
calc. obs. calc. obs. cale. obs. calc. obs. calc. obs. calc’ obs. calc. obs. 
22s (5.36) 5.36 18.03 18.12 37.63 37.74 64.13 64.17 97.33 97.40 137.42 137.43 (184.26) 184.26 
3?s 2.00 2.01 7.21 7.24 15.50 15.51 26.84 26.81 41.18 41.13 58.56 58.64 78.95 78.84 
425 1.04 1.05 3.87 3.88 8.41 8.42 14.67 14.66 22.63 22.60 32.29 32.29 
(ki =6.94; k =1.724) 

2°p (3.53) 3.53 14.21 14.18 31.80 31.78 56.23 56.22 87.46 87.47 125.50 125.54 (170.30) 170.30 
3p 1.55 1.55 6.22 6.22 13.93 13.94 24.68 24.69 38.46 38.46 55.22 55.25 75.02 75.07 
42s 87 87 3.47 3.47 7.78 7.78 13.79 13.80 21.51 21.50 30.91 30.92 42.01 42.02 








states of two-electron atoms: 
E,=H,t+ki(Z—k)Rhc/4-E,/nE,. (4) 


If we employ ‘values of E; obtained from (3) 
by subtracting one-electron atom ionization po- 
tentials, k=.725 as in Table I, and evaluate k, 
from the 2’s term in helium, values of E, are 
obtained which agree with observation within 
+.01 ev in He and +.04 ev in Lit. Other ex- 
perimental data are not available. 

In order to investigate the nature of k; let us 
develop (3) in the same form as (4) by rearrange- 
ment and subtraction of one-electron ionization 
potentials. This gives 


E.=H; 
+[(4Z—3.5)/(Z—k) -1](Z—k)Rhc/4. (5) 


Here ((4Z —3.5)/(Z—k) —1) takes the place of hk, 
in (4), but as noted above it differs from k; by the 
effect of electron coupling. It is interesting that 
this factor while not a constant shows large 
deviation only in H~ and approaches the value 
3.0 for large values of Z. Now in normal He the 
energy assigned to a given electron should be 
39.3 ev rather than 24.47 and 54.14 ev, since both 
the electrons are assumed to be identical. Hence, 
to be strictly comparable to (4), Eq. (5) for He 
would apparently contain the frequency factor 


39.3/24.47 which if factored from the coefficient 
of the resonance term would leave 1.585 which 
agrees closely with the value of k; obtained from 
the 2’s state of He. Similarly in Lit the factor 
98.55/75.3-Z2'Zo/Z, where Zo is the atomic num- 
ber of the first member of the sequence (2 in this 
case), when used as the frequency factor in (4) 
requires the value k;=1.576 to agree with (5). 
Data for excited 1S states in He and Li* calcu- 
lated from (4) using k,=1.58 are shown in 
Table II. The significance of the factor Z’Zo/Z 
which is unity for the first member of a sequence 
is shown below. 

It is of interest to test (4) in various other 
series of He-like atoms. Table III lists calculated 
and observed data for the next three series. The 
values of the constants k; and k which are given 
in the table were obtained from the term values 
shown in parentheses. The agreement appears to 
verify that the form of (4) is essentially correct. 

It is not readily apparent what values of FE, 
should be used in series of three-electron atoms, 
if indeed such methods as those used in He are 
applicable in three-electron atoms. However, we 
may define H, by taking Z’ equal to Z minus the 
number of electrons in the atom core. Also we 
find that the factor 1/n*-E,/H,-Z'/Z is only 
slightly different from the frequency factor used 
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in (4). In fact, it would be almost identical if a 
small constant were subtracted from Z. This 
factor amounts simply to expressing the ioniza- 
tion constants of He-like atoms by Z’Z Rhc where 
the product (Z’Z) is the. square of an effective 
atomic number. Although Z’(Z—.2) gives closer 
agreement, for use in the relative frequency 
factor we may neglect the small constant in Z 
and adjust the constants of (4) to correct for the 
error introduced, thus eliminating the necessity 
of introducing a third constant in (4). We thus 
have 


En,=H,tki(Z—k)Rhc/4-1/n®-En/Hn-Z'/Z, (6) 


which we shall use to express the term values of 
the various series in lithium and like atoms. 
Equation (6) holds as well in He and like series as 
(4) if the constants shown in Tables II and III 
are replaced by the followiny values: 1S: ki; = 1.8; 
k=.725; (P: ky=—.54, R=1.695; 3S: ky =4.14, 
k=.895;°P: ki =1.53, R=1.35. Table IV gives a 
comparison of data calculated from (6) and 
observed data for the 2S and ?P series of Li and 
Li-like atoms. As before, k; and k were evaluated 
from two experimental points in each series. The 
general agreement indicates that the empirical 
frequency factor is essentially correct, and shows 
that the methods used in two-electron atoms are 
valid in three-electron atoms as well. 

To illustrate by a second example the state of 
affairs with coupled electrons, we note that the 
ground states of Be and Be-like atoms are given 
within +.1 ev by the equation 


E=((4(2’) —1)?+6.94(Z —1.724) E/H,,-2'/Z 
— .84(Z —4.0) |Rhc/32. (7) 


The first term in (7) will be recognized as the 


ideal He term for n=2 and Z’ =Z—2. The second 
term is the resonance term for interaction of the 
two 2s (coupled) electrons with the two 1s 
(coupled) electrons. This term was simply taken 
over from the three-electron ground states, the 
effect of the interaction of two instead of one 2s 
electron being accounted for in E/H,. The last 
term evaluated empirically accounts for the 
interaction of the two 2s electrons with each 
other. Note that the ionization potentials of the 
Be sequence are obtained from (7) by subtracting 
the corresponding ionization potentials of the 
lithium sequence. 


DISCUSSION 


It is of significance irrespective of the nature 
of the wave characteristics of matter that one 
can set up simple relations based on classical 
models for expressing the electronic energies, thus 
lending additional support to the belief that the 
Rutherford-Bohr models when properly under- 
stood are valid in atomic and molecular structure. 
That the deviations in many electron atoms from 
the ideal Bohr structure may be expressed by a 
term which is proportional to the orbital fre- 
quency as determined by classical models seems 
quite evident from the results of this study. That 
is, the agreement obtained for higher quantum 
states than that from which the constants were 
evaluated is independent of any advantage which 
may have been gained by the use of empirical 
constants. Finally, one will discover that the 
methods used in two, three, and the ground 
states of four-electron atoms are not limited to 
these special cases although, naturally, compli- 
cations must arise as the number of orbital 
interactions increase. 
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Tables of the thermodynamic functions C,°, H®, S°, and F° are given for 1,3-butadiene and 
the normal butenes for temperatures from 298.16°K to 1500°K for the substances in the ideal 
gas states at one atmosphere pressure. These have been prepared using reliable spectroscopic 
and molecular data, together with calorimetric entropies, gaseous specific heats, and heats of 
formation. Equilibrium constants are given for reactions in thesystem 1,3-butadiene, n-butenes, 
and n-butane. Available experimental data are compared with the calculations. The calori- 
metric data for 1,3-butadiene furnish strong evidence for the existence of two geometric 
(cis-trans) isomeric forms of 1,3-butadiene in appreciable concentrations at room temperatures. 
The two forms differ in energy by 2.3 kcal. mole, and are separated by a C—C rotational 
barrier 2.6 kcal. mole above the lowest energy level of the cis, or higher, energy form. 





I. INTRODUCTION 


a various earlier papers'~‘ results of measure- 
ments of thermodynamic properties of n-bu- 
tane, cis-2-butene, 1-butene, and 1,3-butadiene as 
determined at The Pennsylvania State College 
and the National Bureau of Standards have been 
published. For trans-2-butene the entropy has 
been determined calorimetrically by Guttman 
and Pitzer® at the University of California. In 
the present paper are presented tables of thermo- 


* This paper is a revision of a report “The heat content 
and free energy of butadiene 1,3 and isoprene above 25°C. 
Equilibria in the dehydrogenation of n-butane and n-bu- 
tenes to butadiene-1,3’’ by John G. Aston and George 
Szasz, Pennsylvania State College and F. G. Brickwedde, 
National Bureau of Standards. Submitted to the Office of 
the Rubber Director, July 23, 1943. 

1J. G. Aston and G. H. Messerly, ‘‘The heat capacity 
and entropy, heats of fusion and vaporization, and the 
tideo pressure of n-butane,” J. Am. Chem. Soc. 62, 1917 

1 . 

2R. B. Scott, W. Julian Ferguson, and F. G. Brick- 
wedde, ‘‘Thermodynamic properties of cis-2-butene from 
15° to 1500°K,”’ J. Research Nat. Bur. Stand. 33, 1 (1944). 

3 J. G. Aston, H. L. Fink, A. B. Bestull, E. L. Pace, and 
G. J. Szasz, “The heat capacity and entropy, heats of 
fusion and vaporization and the vapor pressure of butene-1. 
The zero point entropy of the glass. The entropy of the 
gas from molecular data,”’ J. Am. a ce Soc. (in press). 

*KR. B. Scott, R. D. Rands, 7 < Meyers, F. G 
Brickwedde, and N. Bekkedahl, ‘‘Thermodynamic proper- 
ties of 1,3-butadiene in the solid, liquid, and vapor states,”’ 
J. Research Nat. Bur. Stand. 35, 39 (1945). 

5L. Guttman and K. S. Pitzer, ‘‘Trans-2-butene. The 
heat capacity, heats of fusion and vaporization, and vapor 
pressure. The entropy and barrier to internal rotation,” 
J. Am. Chem. Soc. 67, 324 (1945). 


dynamic properties of 1,3-butadiene and the 
three normal butenes in the ideal gas state based 
on spectroscopic data and accepted structural 
constants of the molecules, and such experi- 
mental values of entropy and heat capacity as 
are available. Tables of equilibrium constants 
for the system 1,3-butadiene, m-butenes, and 
n-butane are included together with experi- 
mental equilibrium data where available. 

Values of the physical constants used in the 
calculations are: 


R=1.98714 cal. mole deg.—, 
C2=hc/k =1.4384 cm deg., 

No =6.0228 X 107% mole“, 

Atomic weight of hydrogen = 1.0080, 
Atomic weight of carbon = 12.01. 


The calorie used is the artificial calorie defined 
by 1 calorie =4.1833 international joules. 


II. HEATS OF FORMATION AT 0°K 


The heats of formation of 1,3-butadiene and 
the normal butenes have been obtained from the 
the heat of formation of n-butane and the heats 


‘of hydrogenation of 1,3-butadiene and the 
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butenes to n-butane. 
The heat of formation of n-butane in the 
ideal gas state at 0°K, that is Eo, was obtained 
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by combining enthalpy differences for n-butane,® 
C (graphite),’ and He,’ with the value — 29812 
cal. mole for the heat of formation of n-butane 
from its elements at 298.16°K as given by Prosen 
and Rossini.® 

Values of Eo for 1,3-butadiene and the butenes 
were calculated using: (1) the enthalpy of 
n-butane at 355°K given by the use of Pitzer’s 
tables for (7°—E,°) for n-butane® with the Eo° 
evaluated as previously described, (2) values for 
heats of hydrogenation at 355°K determined by 
Kistiakowsky, Ruhoff, Smith, and Vaughan®?° 
corrected to the new atomic weight of carbon, 
(3) enthalpy differences of 1,3-butadiene and the 
butenes from the tables of the present paper, and 
(4) the enthalpy of He from reference 7. 


III. n-BUTANE 


In view of the greater complexity of the 
n-butane spectrum compared with the spectra 
of the other C, hydrocarbons and the resulting 
greater uncertainty of its interpretation, Pitzer’s 
1944 tables’ of thermodynamic functions of 
n-butane were adopted since they are in agree- 
ment with the calorimetric data on S®° and C,°. 
We have used these tables unchanged except for 
values for 600°K which seemed in doubt on the 
basis of other table values, and have inter- 
polated where necessary. A portion of his table 
values together with interpolated values is given 


in Table I. 


6K. S. Pitzer, ‘Thermodynamics of gaseous paraffins. 
Specific heat and related properties,” Ind. Eng. Chem. 36, 
829 (1944). 

7H. W. Woolley, M. Moskow, F. G. Brickwedde, and 

G. Aston, ‘ ‘Thermodynamics ‘of synthetic rubber and 
related substances. VI. Equilibrium constants of some 
reactions involved in the production of 1,3-butadiene,”’ 
Report to the Office of the Rubber Director, March 24, 
1944. To be published as a paper. 

8 E. J. Prosen and F. D. Rossini, ‘‘Heats of combustion 
and formation of the paraffin hydrocarbons at 25°C,” J. 
Research Nat. Bur. Stand. 34, 263 (1945). 

°G. B. Kistiakowsky, J. R. Ruhoff, H. A. Smith, and 
W. E. Vaughan, “Heats of organic reactions. IV. Hydro- 
genation of some dienes and of benzene,”’ J. Am. Chem. 
Soc. 58, 146 (1936). 

1G. B. Kistiakowsky, J. R. Ruhoff, H. A. Smith, and 
W. E. Vaughan, “Heats of organic reactions. II. Hydro- 
genation of some simpler olefinic hydrocarbons,”’ J. Am. 
Chem. Soc. 57, 876 (1935). 

Tt In the tables of this paper the following conventions 


have been adopted: For enthalpy, H®, and for free energy, . 


F®, the quantities are given as values ‘above assumed zeros 
for the elements in their standard states at 0°K. AH®, AP*, 
and K are, respectively, the enthalpy of formation, the 
free energy of formation, and the equilibrium constant of 
formation of the compound from the elements in their 
standard states at the temperature JT. * 


WOOLLEY, AND BRICKWEDDE 


TABLE I. Free energy and related quantities of n-butane 
in the ideal gas state at 1 atmosphere pressure. AF° and K 
are the free-energy change and equilibrium constant, 


























respectively, for the reaction forming m-butane from -its 
elements in their standard states at T°K (Eo = —23340 cal. 
mole“). 
F° —E,° 

ys -( - ) —F° AF? logio K 
—AF°/RT 
cal. deg. —- 

“— mole~! cal. mole™! cal. mole 2.3026 
298.16 58.51 40780 — 3750 2.747 
300.00 58.61 40920 — 3590 2.613 
400 63.49 48730 +5440 —2.971 
500 67.93 57300 14880 — 6.503 
600 72.06 66570 24590 — 8.957 
700 75.97 76510 34490 — 10.768 
800 79.70 87100 44510 — 12.160 
900 83.28 98290 54590 — 13.257 
1000 86.72 110060 64720 — 14.144 
1100 90.02 122360 74870 — 14.876 
1200 93.19 135160 85060 — 15.491 
1300 96.24 148450 95240 — 16.012 
1400 99.18 162190 105430 — 16.459 
1500 102.03 176380 115620 — 16.846 

IV. 1,3-BUTADIEN E 

Bradacs and Kahovec" investigated the 
Raman spectrum of 1,3-butadiene and con- 


cluded that the structure of this molecule at 
room temperature is ¢rans with no appreciable 
concentration of a cis modification. This con- 
clusion was based on a consideration of the cor- 
respondence between the butadiene spectra and 
the Raman and infra-red activities of the 
normal modes of vibration characteristic of a 
trans molecule, including polarizations of Raman 
lines. Earlier, Schomaker and Pauling” stated 
on the basis of the electron diffraction pattern 
of butadiene that the configuration is essentially 
coplanar and trans although they also stated 
that ‘‘the electron diffraction data are not in- 
compatible with a mixture of cis and (pre- 
dominantly) trans molecules. .. .” Price and 
Walsh!® did not find evidence of a second form 
in ultraviolet absorption and concluded that 
butadiene seemed to be mainly trans. Mulliken" 

11K, Bradacs and L. Kahovec, ‘‘Studies in the Raman- 
effect. Communication 121: Butadiene, homogeneous and 
in solution,” Zeits. f. physik. Chemie B48, 63 (1940). 

2. Schomaker and L. Pauling, “The electron dif- 
fraction investigation of the structure of benzene, pyridine, 
pyrazine, butadiene- 1,3, cyclopentadiene, furan, iasined, 
and thiophene,”’ J. Am. Chem. Soc. 61, 1769 (1939). 

13. W. C. Price and A. D. Walsh, “The absorption spectra 
of conjugated dienes in the vacuum ultraviolet,” Proc. Roy. 
Soc. A174, 220 (1940). 

14 R. S. Mulliken, “Structure and ultraviolet spectra of 


ethylene, butadiene, and their alkyl derivatives,” Rev. 
Mod. Phys. 14, 265 (1942). ‘ 
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indicated on the basis of quantum-mechanical 
calculations that the cis and trans varieties 
should differ by about 2500 calories per mole and 
that the trans variety should be the more stable, 
but he also expressed a doubt as to the relia- 
bility of the estimate of the energy difference. 
Rasmussen, Tunnicliff, and Brattain'® have 
claimed that the cis variety has the lower energy 
and predominates at dry ice temperatures but 
that at room temperature the trans variety pre- 
dominates, basing their view on the existence of 
some coincidences of frequencies in the Raman 
and infra-red spectra and on a rapid change of 
the ultraviolet absorption coefficient with tem- 
perature. 

It seems unlikely, however, that the frequency 
distribution would be so different for cis and 
trans forms of butadiene as to cause the relative 
abundance of the two forms to be overwhelmingly 
reversed in so limited a range of temperature. 
On the contrary, it would seem probable that 
many frequencies of the cis variety would be of 
about the same value as the corresponding fre- 
quencies of the trans variety. 

It may also be pointed out that changes of 
absorption with temperature near the limit of a 
band are not unusual in substances which have 
only a single form and hence do not necessarily 
indicate a rapid change of equilibrium propor- 
tions of cis and trans butadiene. From these 
considerations it is concluded that the observa- 
tions of Rasmussen, Tunnicliff, and Brattain do 
not show that the cis variety is the lower energy 
form. The thermodynamic measurements, in- 
cluding the entropy determined by Scott, 
Rands, Meyers, Brickwedde, and Bekkedahl* 
and the gaseous specific heats determined by 
Aston, Moessen, Hardy, and Szasz,'* Templeton 
and Davies with Felsing,” and Scott and 
Mellors,'!® furnish conclusive evidence for the 
existence of a second variety of 1,3-butadiene 


15R. S. Rasmussen, D. D. Tunnicliff, and R. R. Brat- 
tain, ‘“The configuration of 1,3 butadiene,”’ J. Chem. Phys. 
11, 432 (1943). 

16 J. G. Aston, G. W. Moessen, H. C. Hardy, and G. J. 
Szasz, ‘“‘The heat capacity of gaseous butadiene-1,3,”" J. 
Chem. Phys. 12, 458 (1944). 

17D. H. Templeton and D. D. Davies with W. A. Fel- 
sing, ‘The heat capacity of gaseous 1,3-butadiene from 0 
to 100°,” J. Am. Chem. Soc. 66, 2033 (1944). 

16R. B. Scott and J. W. Mellors, “Specific heats of 
oe 1,3-butadiene, isobutene, styrene, and ethyl- 

enzene,”’ J. Research Nat. Bur. Stand. 34, 243 (1945). 


having an energy of the order of 2000 cal. mole 
higher than the more stable form, as it is impos- 
sible to fit the measured thermodynamic proper- 
ties on the basis of any reasonable frequency 
assignment for a single variety. A paper failing 
to consider the cis-trans equilibrium has recently 
been published by Zeise,!® but the values for the 
thermodynamic properties listed therein are 
unsatisfactory. In the present paper the trans 
has been adopted as the more stable form since 
there is more evidence for this choice than 
against it. 

In calculating the thermodynamic functions 
for butadiene, the cis and trans varieties are 
treated as two chemical compounds in equi- 
librium with each other. The product of the 
principal moments of inertia for over-all rotation, 
the reduced moment of inertia for internal rota- 
tion, and the effective distribution of frequencies 
other than the single bond torsion were taken 
the same for the cis as for the trans form. When 
this is done the equilibrium between cis and 
trans varieties is governed solely by the potential 
energy function assumed to form the barrier for 
the two varieties. 

Contributions to the thermodynamic functions 
coming from internal rotation due to the single 
bond torsion were calculated for each of the two 
varieties using the methods and tables of Pitzer 
and Gwinn.”® For each variety the quantity 
QO; =2.7934 (1088 J,.47)?/n was calculated with 
1/n equal to that fraction of a complete revolu- 
tion of internal rotation to be ascribed to the 
variety in question. For the lower energy form 
the torsional free energy function is obtained 
from Pitzer’s table of (F°—F,°)/T using the 
formula 


(Fane a Eo trans) tors 
, 





Pe (FY F/! 
= ——=R In Q,-———, 
T T 


where the subscript f refers to hypothetical free 


19H. Zeise, “The thermodynamic functions of gaseous 
butadiene and the equilibrium constant for the decom- 
position of butylene into butadiene and hydrogen,” Zeits. 
f. Elektrochemie 50, 113 (1944). 

2K. S. Pitzer and W. D. Gwinn, “Energy levels and 
thermodynamic functions for molecules with internal 
rotation. I. Rigid frame with attached tops,” J. Chem. 
Phys. 10, 428 (1942). 
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rotation. For the higher energy form a similar 
quantity can be calculated 


(F° cis “a, Eo°cis) tors 
T 





F°® (F°—F 0) 
= ——=R In Q,-——_—, 
T 


, 


using constants appropriate to the cis variety. 
To treat the equilibrium between the two forms, 
the quantity 


(Fvis = Eo trans) tors 


T 





(F°—F,) AE, 
T T 





=RIln O;- 


is to be evaluated, where AEo® = Eo cis — Eo trans, 





H®—E,° TH cis— Eo cis 





r-: nt fT 


the estimated difference in energies of the two 
forms in their lowest quantum levels. The equi- 
librium constant for the formation of the cis 
form from the trans form is given by 


K =X cis/X trans 








with 
AF° (F°.is— Eo trans) tors 
Rin K=-—=- 
a T 
(F° trans — Eo" trans) =| 
T ’ 


wherc the X’s refer to mole fractions. The ther- 
modynamic functions for the mixture are ob- 
tained from the complete functions for the two 
varieties using the formulas 





Eo°cie— Eo trans HP trans— Eo" trans 
% trans ’ 


T 











F°— F)° Tr F°.4— Eo cis Eo cis — Eo trans PP? nena Log wane 
-( a )=Xe — oe ]+-Xemd -( )| 
> 4 T T i 


—RXecis In Xcist+X trans In p gn | 





S° _ X cioD cis +X waneS” trans a R[Xeis In X cis +X trans In X wml 
XcisX tranef H cis — Eo cis Eo'cis— Eo" trans BP rene Eee teens r 
+ . |. 
a T T 





Cc. = K ctboe cts +X wensCg* wens + 





The spectroscopic data for 1,3-butadiene normal modes of vibration of the trans variety 





include the Raman spectra of Bradacs and 
Kahovec,! and of Dadieu and Kohlrausch,” and 
the infra-red spectra of Bartholomé and Karweil,” 
and of the American Petroleum Institute.?* Data 
in a more recent paper by Yao*‘ do not lead to 
any change in the assignment. The funda- 
mental frequencies, in cm@', assigned to the 


21 A, Dadieu and K. W. F. Kohlrausch, ‘‘Raman effect 
and constitutional problems,” Ber. d. d. chem. Ges. 63, 
1657 (1930). 

2 E. Bartholomé and J. Karweil, ‘The infra-red spec- 
trum and the fundamental vibration bands of the hydro- 
carbons biacetylene, butadiene and vinylacetylene,”’ Zeits. 
f. physik. Chemie B35, 442 (1937). 

3 American Petroleum Institute, “Infra-red absorption 
spectrograms,” issued by the American Petroleum In- 
stitute Research Project 44 at the National Bureau of 
Standards. 

*Y. T. Yao, “Absorption bands in the spectrum of 
butadiene,” paper presented at the Columbus meeting of 
the American Physical Society, June 15, 1945. 


of the butadiene molecule (point group C2) are 
as follows, where the symbols have the same sig- 
nificance as in Kohlrausch’s book.?5 The symbols 
(p) and (dp) denote polarized and unpolarized 
Raman lines while (i.r.) denotes an infra-red 
band. The values given in parentheses were not 
observed spectroscopically but were selected on 
other considerations. 


Carbon skeleton frequencies 
Ag: w1, 1635 (p); we, 1200 (p) ; ws, 513 (p) 
By: ws, 1602 (t.r.) ; ws, (326). 


y(C—H) rocking frequencies 
A,: 520 (4.r.), 910 (z.r.), 1014 (4.r.) 
B,: (667), 899 (dp), 914 (dp). 
25K. W. F. Kohlrausch, Der Smekal-Raman Effekt 


(Verlagsbuchhandlung, Julius Springer, Berlin, 1931-1937), 
Supplement, pp. 44 and 69. 
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5(C—H) bending frequencies 
A,: 1276 (p), 1303 (p), 1437 (p) 
B,,: 1282 (i.r.), 1385 (i.r.), 1490 (.r.). 


v(C—H) stretching frequencies 
A,: 3000, 3088, 3000 
B,,: three 3000. 


Except for several frequencies which are here 
mentioned explicitly, the assignment given is 
that of Bradacs and Kahovec." Values for the 
torsional frequencies associated with the double 
bonds have been selected on the basis of the 
corresponding value for ethylene. The estimates 
were 587 cm! for the A, mode and 667 cm for 
the B, mode. The 587 cm=! was shifted to the 
observed 520 cm~ and the calculated value 667 
cm~ has been retained unchanged. Force con- 
stant calculations were also made for the 
skeletal deformation vibrations. They indicated 
that for w;=513 cm™, ws would be expected to 
be about 345 cm. This value was lowered to 
326 cm in fitting the experimental thermo- 
dynamic data. Some of the frequencies used for 
the trans molecule may belong to the cis variety. 
A few observed frequencies were not used but 
it is not difficult to account for them in a way 
that seems plausible. The weak infra-red line at 
665 cm~ is presumed to belong to the cis form. 
The very weak Raman line at 1782 cm-! may 
be the overtone of the 899 cm fundamental and 
the weak Raman line at 1025 cm is regarded 
as the overtone of w3, 513 cm~!. A weak Raman 
line occurring at 340 cm~ is taken to be caused 
by an overtone of a frequency at 170 cm~ for 
the single bond torsion for the trans form. Calcu- 
lations of the type recommended by Lord” using 
low temperature specific heat data for crystalline 
butadiene favored the use of a frequency of this 
general magnitude for the lowest molecular 
frequency. Preliminary adjustments to fit S° 
and C,° at 300°K also indicated essentially this 
value. 

The values used for bond lengths and angles 
were selected on the basis of the work and 
recommendations of Schomaker and Pauling’ 
and of Galloway and Barker ;2”7 C—H, 1.09A; 

26R. C. Lord, Jr., “The heat -capacities of molecular 
lattices. I. Introduction,” J. Chem. Phys. 9, 693 (1941). 

27W. S. Galloway and E. F. Barker, “The infra-red 


absorption spectra of ethylene and tetra-deutero-ethylene 
under high resolution;’’ J. Chem. Phys. 10, 88 (1942). 


C-—C, 146A; <C—C=C, 124°; <H—C=C, 
115°; <H—C—H, 120°. These lead to the value 
ABC=0.8186 X10-""* g* cm® for the product of 
the principal moments of inertia and to 10.282 
X10-*° g cm? for the reduced moment of inertia 
for internal rotation for the trans variety of 
butadiene. Thus the expressions for the entropy, 
free energy function, and heat capacity of trans- 
lational and molecular rotational motions in cal. 
deg.—! mole are 


S,,,= 18.3022 logio 7+14.835, (1) 


F°—EF,° H—E,)°® 
(EF) nO) 
T t+r T t+r 


= Re ead 4R 
= 18.3022 logis T+6.886, (2) 
C°» t4,= 7.949. (3) 


The addition to the above Eqs. (1) and (3) of 
the intra-molecular vibrational . contributions 
apart from the undetermined w; and the single 
C—C bond torsion yields values which, when 
subtracted from the calorimetric values of S° and 
C,°, are the entropy and specific heat con- 
tributions of ws; and the torsional oscillation of 
the two vinyl groups about the middle C—C 
bond. A simultaneous fit of entropy, specific 
heat, and change of specific heat with tem- 
perature was obtained, working with ws and the 
parameters Vo, trans, Vo, cis, and 


1/Neis =1-— (1/n trene) 


of barriers for internal rotation of the trans and 
cis forms about the C—C bond where the bar- 
riers are given by the equation 


V=0.5 Vol1—cos n@]. 


The resulting shape of the barrier is shown as 
the continuous curve in Fig. 1. The constants 
found were 1/n=0.406 and V»p=5000 cal. mole 
for the trans barrier and 1/n=0.594and Vo=2575 
cal. mole for the cis barrier while ws; was 
adjusted to 326 cm~. Thus, 40.6 percent of the 
angle of rotation about the central C—C bond 
was assigned to the trans form and 59.4 percent 
of the angle to the cis form. When allowance is 
made for zero point energy, the lowest quantum 
level of the cis form is approximately 2300 cal. 
mole higher than the lowest quantum level 
of the trans form. While significance is attached 
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ANGLE OF INTERNAL ROTATION 


Fic. 1. Potential energy as a function of the angle of 
torsion about the central C—C bond in 1,3-butadiene. The 
continuous curve represents the barriers used in the cal- 
culations of Tables II and III: 


Vtrans = 2.500[.1—cos (0/0.406)] in kcal. mole! 
and 


Veis = 2.425 + 1.288[1 —cos ((@— a ]/0.594) ] in kcal. mole. 


The dashed curve represents a modified potential energy 
function giving the same thermodynamic properties. 


to the general features of the cis-trans barrier as 
given here, it must be pointed out that the exact 
details are quite uncertain. Thus, in Fig. 1, the 
dashed curve has been so constructed as to 
provide essentially the same over-all distribution 
of levels with regard to energy (and therefore the 
same thermodynamic functions) as occurs for 
the continuous curve. For the dashed curve, 50 
percent of the angle of rotation belongs to the 
trans form and 50 percent to the cis form, as 
might be expected from the theory of bond con- 
jugation. With the dashed curve, the concen- 
tration of the cis form would be slightly less 
than with the continuous curve. 

Calculated and experimental values of the 
specific heat of 1,3-butadiene in the ideal gas 
state are plotted in Fig. 2. In Tables II and III 
are given the calculated thermodynamic func- 
tions of 1,3-butadiene and the mole fraction of 
its cts form. 


V. TRANS-2-BUTENE 


The Raman and infra-red spectra of trans-2- 
butene have been photographed by Gershinowitz 
and Wilson.*8 

28 H. Gershinowitz and E. B. Wilson, Jr., ‘Infra-red and 


Raman spectra of polyatomic molecules. II. Cis- and 
trans-butene-2,”’ J. Chem. Phys. 6, 247 (1938). 
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Fic. 2. Specific heat of gaseous 1,3-butadiene in the ideal 
gas state. 


@ Scott and Mellors 

@ Templeton and Davies with Felsing 

Oo, ©, o Aston, Moessen, Hardy and Szasz, data 
variously corrected to zero pressure. 

The curve represents the calculated values of this paper. 


This molecule, like butadiene, belongs to the 
symmetry class C2, which allows no coincidences 
between the Raman and infra-red spectra. There 
are (excluding internal rotation) twenty-eight 
modes of vibration to which frequencies must be 
assigned. With the same notation used for 
butadiene, the assignment adopted was as 
follows: 


Carbon skeleton frequencies: 


Ag: w1, 1681 (p); we, 870 (p); ws, 507 (p) 
B,: ws, 978 (i.7.); ws, (313) (4.7.) ; we, (270). 


Frequencies due to CH groups: 


Ag: v1, 2739 (p); 61, 1309 (p) 
Au: ¥1, 1080 (7.r.) 

B,: Y¥2, 1043 (dp) 

By: ve, 2427 (4.r.); 62, 1326 (4.7.). 


Frequencies due to CH; groups: 


A,: v3(m), 2867 (p) ; 63(), 1309 (p) ; va(a), 2965 
(p); 64(o), 1309 (p); vs, 870 (p) 

Ax: vs(c), 2976 (i.r.) ; 55(¢), 1466 (2.7.) ; ya, 1072 
(z.r.) ; internal rotation 

B,: ve(a), 2926 (dp); d6(c), 1455 (dp); v5, 746 
(dp) ; internal rotation 

B,,: v(m), 2976 (i.r.); 57(ar), 1560 (2.7.); va(o), 
2976 (i.r.); 53(0), 1488 (t.r.); ye, 1214 (7.r.). 
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The frequency ws was not observed but filled in 
by analogy with other molecules. The frequency: 
ws was too low to be observed in the infra-red 
and the value 313 cm obtained by the normal 
coordinate treatment of O. Burkard*’ was chosen. 
The symbols z and oe are used to denote vibra- 
tions in which both methyl groups vibrate 
essentially as they do in methyl chloride vibra- 
tions of the w and o type, respectively. These 
symbols in addition to the symmetry class 
uniquely specify the types of vibration. The 
symbols (p) and (dp) denote polarized and 
unpolarized Raman lines while (i.r.) denotes an 
infra-red line (or band). 

The distances and angles used were: C—C, 
1.54A; C—H, 1.09A; C=C, 1.36A; and <C— 
C—H, 109°28’. Thus the product ABC of the 
principal moments of inertia is 1.449X10-"4 g3 
cm®, and the reduced moment of inertia Iyeq of a 
methyl group is 5.0X10- ¢ cm’. 

The entropy and free energy of the translation 
and external rotation of the molecule together with 
free internal rotations of the methyl groups are 


Se4r4fe= 22.875 log T+11.239, (4) 


F°— E,° 
™ (—) = 22.875 log T+1.305. (5) 
i t+r+fr 


TABLE II. Heat capacity at constant pressure, heat 
content and related quantities of 1,3-butadiene in the ideal 
gas state, and the mole fraction X.:; of the cis form of 1,3- 
butadiene in the mixture of cis and trans forms. AH® is 
the enthalpy of formation of 1,3-butadiene from its ele- 
ments in their standard states at the temperature T 
(Eo? = 30200 cal.mole“). 











H® —E»° 
i Xeis H° AH? Cp® 
T 
cal. deg.~! cal. deg.~! 
°K mole cal. mole~! cal. mole“! mole 





298.16 0.04 12.16 33830 26750 19.01 
300.00 04 12.20 33860 26730 19.11 


400 .09 14.59 36040 25830 24.29 
500 16 16.97 38690 25120 28.52 
600 abe 19.18 41710 24540 31.84 
700 27 ~~ ~21.19 45040 24080 34.55 
800 32 23.01 48610 23710 36.84 
900 35 24.66 52390 23420 38.81 
1000 38 26.16 56360 23210 40.52 
1100 40 27.54 60490 23060 42.02 
1200 42 28.80 64760 22980 43.32 
1300 44 29.96 69150 22940 44.47 
1400 A5 31.03 73650 22950 45.47 
1500 46 32.02 78240 22980 46.34 





Adding to S°.,,4;, the contributions made by 
the internal vibrations, an entropy value was 
obtained which was used with Pitzer’s experi- 
mental value,® 69.12 cal. deg.-'! mole, of the 
entropy at 274.1°K to calculate the potential 
restricting rotation of the methyl groups. The 
value found is 1900 cal. mole. 

The calculated heat capacities are compared 
in Table IV with the calorimetric data of Kis- 
tiakowsky and Rice.* 

Values of the thermodynamic functions for 
trans-2-butene are summarized in Tables V and 
VI. 

VI. CIS-2-BUTENE 


The thermodynamic functions of cis-2-butene 
are from reference 2 using a new value for E,° 
and are based on spectroscopic, molecular, and 
calorimetric data discussed in that publication. 
The change of value for Eo° is caused by changes 
in the values used for the heat of formation and 
thermodynamic functions of n-butane. 

Values of thermodynamic functions of cis-2- 
butene are given in Tables VII and VIII. 


VII. 1-BUTENE 


The entropy of 1-butene in the standard ideal 
gas state at the normal boiling point, 266.91°K, 


TaBLE III. Entropy, free energy, and related quantities 
of 1,3-butadiene in the ideal gas state at 1 atmosphere 
pressure. AF’ and K are the free energy change and equi- 
librium constant, respectively, for the reaction forming 
1,3-butadiene from its elements in their standard states at 
T°K (Eo°=30200 cal. mole). 


F0 —Eo° 
T -( - ) S° —Fo AF logio K 











—~AF°/RT 





cal. deg.~! cal. deg. 


"i. mole~! mole~! cal. mole“! cal. mole 2.3026 








298.16 54.46 66.62 —13960 36430 —26.707 
300.00 54.44 66.74 —13840 36490 —26.587 


400 58.38 72.97 —6850 39890  —21.793 
500 61.89 78.86 +740 434909 —19.009 
600 65.18 84.36 8920 47210 —17.198 


700 68.29 89.48 17600 51040 —15.934 
800 71.24 94.25 26790 54920 —15.003 
900 74.05 98.71 36440 58830  —14.287 
1000 76.72 102.88 46520 62780 —13.721 
1100 79.28 106.82 57010 66740  —13.260 
1200 81.73 110.53 67880 70720 —12.880 
1300 84.09 114.05 79110 74690 —12.557 
1400 86.35 117.38 90680 78670 —12.281 
1500 88.52 120.54 102580 82660  —12.043 





22K. W. F. Kohlrausch, Der Smekal-Raman Effekt 
(Verlagsbuchhandlung, Julius Springer, Berlin, 1938), 
Supplement, p. 10. 








 %G. B. Kistiakowsky and W. W. Rice, “Gaseous heat 
capacities. II,” J. Chem. Phys. 8,610 (1940). 
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TaBLe IV. A comparison of the calculated and experi- 
mental values for the heat capacity of trans-2-butene in 
the ideal gas state. 











Yt Cp® (calc.) Cp® (expt.) 

= cal. deg.~! mole cal. deg.~! mole 
298.6 20.90 20.98 20.81 
332.9 22.61 22.69 22.61 
371.5 24.51 24.53 24.51 








8 For the ideal gas state as reported by Kistiakowsky and Rice using 
Roper’s equation of state (see reference 30). 

b For the ideal gas state using experimental data of Kistiakowsky 
and Rice and the equation of state for cis-2-butene (see reference 2). 


has been determined calorimetrically* as 70.87 
+0.21 cal. deg.-! mole. Preliminary values of 
the specific heat of gaseous 1-butene in the ideal 
gas state as determined by Scott and Wacker*! 
are 21.46 cal. deg.-! mole at 313.56°K and 
24.18 cal. deg.—! mole at 363.26°K. With these 
data it is possible to evaluate the thermodynamic 
properties from the incomplete spectroscopic 
data with more certainty than originally obtained 
by Pitzer.” 

The symmetry class to which this molecule 
belongs depends on the relative position of the 
ethyl and vinyl groups. If all the atoms are 
coplanar it is Cy, = C,, otherwise, it is C;. In either 
case, all lines should appear in the Raman 
spectrum and in the infra-red, but in the former, 
Raman lines caused by modes symmetrical with 
respect to the plane of symmetry would be 
polarized. Unfortunately, no polarization meas- 
urements have been reported for the Raman 
spectra so that this important clue concerning 
the form of the molecule is lacking. 

The Raman and infra-red spectra indicate 
essentially only one form. A calculation of the 
barrier for the rotation of vinyl and ethyl groups 
with respect to each other based on hydrogen 
repulsions alone* would indicate two non-planar 
stable configurations which are optical isomers 
lying on either side of the planar trans con- 
figuration. Each of these non-planar forms 
occupies about one third of a revolution. The 


aR. B. Scott and P. F. Wacker, Nat. Bur. Stand. 
unpublished data. 

®K. S. Pitzer, ‘‘Thermodynamics of gaseous hydro- 
carbons: ethane, ethylene, propane, propylene, m-butane, 
isobutane, 1-butene, cis- and trans-2-butenes, isobutene, and 
co jo (tetramethylmethane),” J. Chem. Phys. 5, 473 
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Raman spectrum does not indicate the presence 


-of an appreciable concentration of a second or 


planar form. 

The Raman spectrum of this compound has 
been studied by Kohlrausch and Stockmair,** 
while the infra-red data have been obtained from 
American Petroleum Institute Spectrograms.”? 

From these results the following frequency 
assignment has been made: 

cm 
C=C stretching 1636 
C—C stretching : 850 and 1068 
CH; internal angle change 1450(2), 1376(1) 


skeleton angle change 320, 437 
= CH: internal angle change 1416 

= CH; torsion 623 

CH wagging in plane 1293 

CH wagging out of plane 908 

CH; wagging 972, 1019 
— CH; twisting 1257 

CH2 wagging (940), 1315 
CH: internal angle change 1450 

= CH: wagging 995, 1166 
CH stretching 2950(8) 


In making this assignment, the assignments 
made for propane by Pitzer*‘ and for propylene 


TABLE V. Heat capacity at constant pressure, heat 
content and related quantities of trans-2-butene in the 
ideal gas state. AH® is the enthalpy of formation of trans-2- 
butene from its elements in their standard states at the 
temperature T (Eo®=2520 cal. mole). 














H°—E,® 
r Ho AH® Cp® 
= 
cal. deg.-! cal. deg.~1 
. 3 mole! cal. mole=! cal. mole! mole! 
298.16 14.03 6700 — 2400 20.88 
300.00 14.07 6740 — 2420 20.97 
400 16.40 9080 — 3860 25.89 
500 18.78 11910 — 5090 30.57 
600 21.11 15180 —6120 34.78 
700 23.32 18850 — 6950 38.48 
800 25.42 22860 — 7580 41.68 
900 27.39 27170 — 8050 44.47 
1000 29.22 31740 — 8380 46.89 
1100 30.92 36540 — 8580 48.99 
1200 32.51 41530 — 8680 50.80 
1300 33.97 46690 — 8690 52.36 
1400 35.32 51970 — 8660 Sa.74 
1500 36.61 57430 — 8530 54.89 








33K. W. F. Kohlrausch and W. Stockmair, “Studies in 
the Raman effect. Communication 44. Raman effect and 
free rotation. III,” Zeits. f. physik. Chemie B29, 292 (1935). 

34K. S. Pitzer, ‘‘The molecular structure and thermo- 
dynamics of propane. The vibration frequencies, barrier to 
internal rotation, entropy, and heat capacity,” J. Chem. 
Phys. 12, 310 (1944). 
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TasBLe VI. Entropy, free energy, and related quantities 
of trans-2-butene in the ideal gas state at 1 atmosphere 
pressure. AF° and K are the free energy change and equi- 
librium constant, respectively, for the reaction forming 
trans-2-butene from its elements in their standard states 
at T°K (E.°=2520 cal. mole). 


F° —Ep° 
T -( ; ) S90 —Fo AFo logio K 











—AF°/RT 
cal. deg.~! cal. deg.~! ———- 
mole~ mole cal. mole cal. mole 2.3026 


298.16 56.91 70.95 14450 §=15300 —11.217 
300.00 57.00 71.07 14580 15410  —11.228 





400 61.38 77.78 22030. 21570 —11.784 
500 65.29 84.07 30120 28080  —12.274 
600 68.91 90.02 38830 34810 —12.681 
700 72.33 95.65 48110 41710  —13.022 
800 75.60 101.02 57960 48700 —13.304 
900 78.69 106.08 68310 55770 —13.543 
1000 81.68 110.90 79160 62880 —13.742 
1100 84.54 115.47 90480 70010  —13.910 
1200 87.30 119.81 102250 77160 —14.053 
1300 89.96 123.94 114430 84310  —14.174 
1400 92.54 127.87 127040 91440  —14.275 
1500 95.01 131.62 140000 98620 —14.369 








by Wilson and Wells** were considered. The 940 
cm~! frequency which has been used for a CHe 
wagging motion was not found in the 1-butene 
spectra but has been used to be consistent with 
the schematic frequency system for a CH: group 
as indicated by Pitzer. A 1-butene infra-red line 
at 1130 cm was under consideration for this 
motion but it was found in only one of several 
infra-red investigations and, if not due to an 
impurity, may well have been a difference 
frequency, since 1450 —320=1130. 

The bond lengths and angles used were: C—C, 
1.544; C=C, 1.36A; C—H, 1.09A; <C—C—C, 
109°28’; <C—C—H, 109°28’. The average of 
the products of moments of inertia about prin- 
cipal axes of the planar and non-planar forms is 
then 

ABC=1.53X10—" g? cm®. 


The reduced moment (using axes through the 
centers of gravity of the groups parallel to the 
bond) for the mutual rotation of the ethyl and 
vinyl groups is 


Trea = 13.37 XK 10-” g cm?, 
while that for the internal rotation of the methyl 


% E. B. Wilson, Jr., and A. J. Wells, “Infra-red and 
Raman spectra of polyatomic molecules. XIV. Propylene,”’ 
J. Chem. Phys. 9, 319 (1941). 


TaBLe VII. Heat capacity at constant pressure, heat 
content and related quantities of cis-2-butene in the ideal 
gas state. AH” is the enthalpy of formation of cis-2-butene 
from its elements in their standard states at the tempera- 
ture T (Eo? =3710 cal. mole). 














H°—E»° 
i H° AH? Cp® 
r 
cal. deg.1 cal. deg.~1 

°K mole~! cal. mole~! cal. mole~! mole~! 
298.16 13.50 7740 — 1360 19.17 
300.00 13.53 7770 — 1400 19.27 
400 15.63 9960 — 2970 24.56 
500 17.92 12670 — 4330 29.57 
600 20.24 15860 — 5450 33.96 
700 22.48 19450 — 6340 37.74 
800 24.60 23390 — 7050 41.00 
900 26.59 27640 — 7580 43.83 
1000 28.44 32150 — 7970 46.28 
1100 30.16 36890 — 8230 48.41 
1200 31.76 41820 — 8390 50.25 
1300 33.24 46930 — 8450 51.85 
1400 34.62 52180 — 8450 53.23 
1500 35.90 57560 — 8400 54.44 








TABLE VIII. Entropy, free energy, and related quantities 
of cis-2-butene in the ideal gas state at 1 atmosphere 
pressure. AF’ and K are the free-energy change and equi- 
librium constant, respectively, for the reaction forming 
cis-2-butene from its elements in their standard states at 
T°K (E.°= 3710 cal. mole). 











F°—Ep° 
 ¥ -( > ) S° —Fo AFo logio K 


—AF°/RT 
cal. deg.~! cal. deg. sntmiamennatin 
a mole! mole“! cal. mole! cal. mole=! 2.3026 











298.16 58.32 71.82 13680 16070  —11.782 
300.00 58.41 71.94 13810 16180  —11.787 
400 62.59 78.22 21330 22270 = —12.170 
500 66.32 84.24 29450 28750  —12.567 
600 69.79 90.03 38160 35480 —12.923 
700 73.07 95.55 47440 42380 —13.231 
800 76.21 100.81 57260 49400 —13.494 
900 79.22 105.81 67590 56480 —13.716 
1000 82.12 110.56 78410 63630 —13.906 
1100 84.91 115.07 89690 70800 —14.066 
1200 87.61 119.37 101420 77990 —14.204 
1300 90.21 123.45 113560 85180 —14.320 
1400 92.73 127.35 126110 92380 —14.421 
1500 95.16 131.06 139030 99580 —14.509 
group Is 


Trea = 5.0K 10-” g cm. 


A potential of 2700 cal. mole“! was assumed 
to hinder the internal rotation of the methyl 
group. The expression for the entropy of trans- 
lation and total rotation (including free internal 
rotation) in cal. deg.—! mole, is 


Stress int rot = 22.875 log T+15.048- 
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TABLE IX. A comparison of calculated and experimental 
specific heats and entropies for 1-butene in the standard 
ideal gas state. 




















x Cp® (calc.) Cp® (expt.) 
“——. cal. deg.~! mole! cal. deg.~! mole™! 
313.56 21.52 21.46 
363.26 24.17 24.18 
_ = * tT) 7 S°® (calc.) S° (expt.) 
°K cal. deg.~1 mole7! cal. deg. inole™ 
200.00 65.88 66.24 
230.00 68.18 68.49 
266.91 70.87 


70.87 











At the normal boiling point (266.91°K) the 
experimental entropy of the ideal gas is 70.87 
cal. deg.! mole. The vibrational entropy, cal- 
culated from the above spectrum, is 2.75 cal. 
deg.—! mole while the entropy due to translation 
and total rotation without correction for re- 
striction is 70.55 cal. deg.! mole. The re- 
stricting correction for the methyl group rotation 
obtained from the tables of Pitzer and Gwinn” 
is 1.33 cal. deg.—! mole. The correction for the 
restriction of the ethyl group rotation must thus 
be 1.10 cal. deg. mole. A barrier of about 
2000 cal. mole gives essentially this correction. 
In order also to fit the experimental gaseous 
specific heats, further adjustment of details of 
the barrier was required. A plot of the potential 
energy calculated on the basis of hydrogen re- 
pulsions’ suggests that in the region of the planar 
trans form there is so little change of energy with 
angle that this region may be considered as an 
energy plateau. This plateau-like region was 
taken as occupying 12 percent of a revolution 
with each of the two non-planar forms occupying 
30 percent of a revolution and having potential 
energy valleys of 1970 cal. mole. The calcula- 
tion procedure for a mixture of different forms 
of molecule has been indicated in the discussion 
of the calculations for butadiene. 

In the paper by Aston et al.* the entropy is 
given not only for 266.91°K but also for 200°K 
and 230°K. As will be seen from Table IX, it is 
found that when the experimental specific heats 
at 313.56°K and 363.26°K and the entropy at 
266.91°K are fitted, the agreement of the calori- 
metric and calculated entropies for 200°K and 
230°K_ is not exact. From one standpoint it 
seemed preferable to fit the value reported for 
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266.91°K since in this temperature region the 
reported thermal data satisfy the Clapeyron 
equation better than at the two lower tem- 
peratures. 

Tables X and X]I list the thermodynamic prop- 
erties of 1-butene calculated as outlined in this 
section. 


VIII. EQUILIBRIA 


1. Equilibrium Constants for the Isomerization 
of 1-Butene(g) to 2-Butene(g) 


From the data discussed in the preceding 
sections the equilibrium constants have been 
calculated for the reactions 


1-butene(g)—#trans-2-butene(g), 
1-butene(g)—cis-2-butene(g). 


The sum of the equilibrium constants for these 
two reactions is the equilibrium constant for the 
reaction 


1-butene(g)—equilibrium 2-butene(g). 


The results for the three reactions are sum- 
marized in Table XII. 

In Fig. 3 the results for the reaction 1-butene— 
equilibrium 2-butene are shown in comparison 
with experimental data determined by McCarthy 
and Turkevich,** Frey and Huppke,*” Runge and 

TABLE X. Specific heat at constant pressure, heat con- 
tent and related quantities of 1-butene in the ideal gas 
state. AH® is the enthalpy of formation of 1-butene from 


its elements in their standard states at the temperature T 
(E,° = 5370 cal. mole). 











(H° —Eo") 
ri = — - H? AH? Cp? 
cal. deg.7! cal. deg.-! 
“=. mole~! cal. mole! cal. mole! mole! 
298.16 13.60 9425 325 20.69 
300.00 13.64 9460 300 20.79 
400 16.09 11810 — 1130 26.07 
500 18.59 14660 — 2340 30.90 
600 20.99 17960 — 3340 35.08 
700 23.27 21660 — 4130 38.68 
800 25.39 25680 — 4760 41.79 
900 27.37 30000 — 5220 44.50 
1000 29.20 34570 — 5550 46.85 
1100 30.90 39360 — 5760 48.90 
1200 32.48 44350 — 5860 50.68 
1300 33.94 49490 — 5890 52.22 
1400 35.29 54780 — 5850 53.57 
1500 36.56 60210 — 5750 54.74 








36 W. W. McCarthy and J. Turkevich, “Butene-1 and 
butene-2 equilibrium,” J. Chem. Phys. 12, 405 (1944). 

37F, E. Frey and W. F. Huppke, “‘Equilibrium dehy- 
drogenation of ethane, propane, and the butenes,” Ind. 
Eng. Chem. 25, 54 (1933). 
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Mueller-Conradi,* and Hurd and Goldsby*® 
plotted as logio K against a scale linear in 1/T. 
The agreement of the calculated curve with the 
data of McCarthy and Turkevich, which appear 
more reliable than the other equilibrium data, 
would be improved appreciably if for 1-butene 
the calorimetric value of the entropy had been 
fitted at 200°K rather than at 266.91°K but it is 
thought that other uncertainties in both the 
experimental and calculated data are sufficient 
to account for the differences between observed 
and calculated equilibrium constants. Thus the 
fit of the experimental and calculated equilibrium 
constants could be improved by small adjust- 
ments in the cis- and trans-2-butene data. Com- 
parisons of calculated equilibrium constants 
with reliable experimental data for reactions 
having small heats of reaction, as the isomeriza- 
tion reaction in question, provide stringent tests 
for the accuracy of thermodynamic data. 


2. Equilibrium Constants for the Dehydro- 
genation of n-Butane(g) to n-Butenes(g) 
and 1,3-Butadiene (g) 


From the tables of the preceding sections the 
equilibrium constants (K, for pressures in at- 


TABLE XI. Entropy, free energy, and related quantities 
of 1-butene in the ideal gas state at 1 atmosphere pressure. 
AF® and K are the free-energy change and equilibrium 
constant, respectively, for the reaction forming 1-butene 
from its elements in their standard states at T°K 
(Eo = 5370 cal mole“). 




















Fo —F,° 
i -( ) So —F° AF® logio K 
r 
—AF°/RT 
cal. deg.~! cal. deg. - 

a 4 mole! mole! cal. mole! cal. mole! 2.3026 
298.16 59.47 73.07 12365 17390 — 12.745 
300.00 59.56 73.20 12500 17490 — 12.745 
400 63.82 79.91 20160 23440 — 12.808 
500 67.67 86.26 28460 29740 — 12.999 
600 71.28 92.27 37400 §=©©36240 — 13.201 
700 74.69 97.96 46910 42910 — 13.396 
800 77.94 103.33 56980 49670 — 13.570 
900 81.05 108.42 67580 56500 — 13.720 

1000 84.03 113.23 78660 63380 — 13.851 
1100 86.89 117.79 90210 70280 — 13.963 
1200 89.64 122.42 102200 77210 — 14.062 
1300 92.30 126.24 114620 84120 — 14.143 
1400 94.87 130.16 127450 91040 — 14.212 
1500 97.34 133.90 140640 97980 — 14.275 








38F, Runge and M. Mueller-Conradi, U. S. Patent 
1,914,674, June 20, 1933. 

3C. D. Hurd and A. R. Goldsby, ‘Rearrangement 
during pyrolysis of the butenes,”’ J. Am. Chem. Soc. 56, 
1812 (1934). 
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_ Fic. 3. Logarithm of the equilibrium constant for the 
isomerization of 1-butene to 2-butene. Scale of abscissas 
is linear in 1/T. 


T - . , . , 
° McCarthy and Turkevich, with their estimated range 


of error 
© Frey and Huppke 
@® Runge and Mueller-Conradi (I. G. Farben) 
e@ Hurd and Goldsby. 
The curve represents the calculated values of this paper. 


mospheres) have been calculated for the four 
reactions: 


n-butane(g)—trans-2-butene(g) + H2(g), 
n-butane(g)—cis-2-butene(g) + H2(g), 


_ TABLE XII. Equilibrium constants K for the isomeriza- 
tion of 1-butene to 2-butenes: 1-butene—2-butene (trans, 








i tvans-2-butene cts-2-butene mixed-2-butenes 
"x. 
298.16 33.7 9.18 42.9 
300.00 32.9 9.08 42.0 
400 10.6 4.34 14.9 
500 5.31 2.70 8.01 
600 3.32 1.90 5.22 
700 2.37 1.46 3.83 
800 1.85 1.19 3.04 
900 1.50 1.01 2.51 
1000 1.28 0.881 2.17 
1100 1.13 -789 1.92 
1200 1.02 Bey 4 1.74 
1300 0.931 -665 1.60 
1400 .865 .618 1.48 
1500 .805 .583 1.39 
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n-butane(g)—1-butene(g) + H2(g), 
n-butane(g)—1,3-butadiene(g) +2H2(g). 


The results are summarized in Table XIII. 

Experimental values for the n-butane-n-butene 
equilibrium constants were determined by Frey 
and Huppke.*” A comparison with their values is 
shown in Table XIV. At the highest temperature, 
723°K, the agreement is a little poorer than 
would be expected on the basis of the experi- 
mental accuracy originally claimed by Frey and 
Huppke. At the lowest temperature, 623°K, the 
experimental values of K, are higher than the 
calculated by approximately the factor two. Con- 
sidering the small concentrations of components 
that had to be measured in determining the 
experimental values it is not surprising that 
there is a discrepancy between measured and 
calculated values. 


3. Equilibrium Constants for the Dehydrogena- 
tion of the n-Butenes(g) to 1,3-Buta- 
diene (g) : 


Equilibrium constants calculated from our 
data for the three reactions 


trans-2-butene(g)—1,3-butadiene(g) + H2(g) 
cis-2-butene(g)—1,3-butadiene(g) + H2(g) 
1-butene(g)—>1,3-butadiene(g) + H2(g) 


are presented in Table XV. 

The equilibrium constants for the reaction 
mixed butenes(g) (in equilibrium)—1,3-buta- 
diene(g) +H2(g) involved when a mixture of cis- 
and trans-2-butene and 1-butene is in equilibrium 
with 1,3-butadiene are given in the last column 
of Table XV. This constant is calculated from 


TABLE XIII. Equilibrium constants for the gas phase 
dehydrogenation of n-butane to butenes or to 1,3-buta- 
diene: n-butane—butene (or 1,3-butadiene) and hydrogen. 











x Kp in atmos. Kp in atmos.? 
4 trans-2-butene cis-2-butene 1-butene 1,3-butadiene 
298.16 1.09 X10-14 2.96 K10715 3.22 K10716 3.52 K10-30 
300.00 1.44 X1074 3.98 K107 4.39 K10716 6.31 K10~30 
400 1.54 X1079 6.32 K10710 1.46 X10710 1.51 X10719 
500 1.69 X1076 8.63 X1077 3.19 K1077 3.12 K10713 
600 1.89 X1074 1.08 X10~4 5.70 K1075 5.74 K10~9 
700 5.57 X10-3 3.44 X1073 2.40 X1073 6.82 X10-6 
800 7.18 X1072 4.63 X1072 3.89 X1072 1.44 K10-3 
900 0.518 0.347 0.344 9.33 K1072 

1000 2.52 1.73 1.96 2.65 
1100 9.25 6.46 8.18 41.3 

1200 27.4 19.4 26.9 408 
1300 68.9 49.2 74.0 2850 
1400 153 109 177 15100 
1500 300 217 372 63500 
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Fic. 4. Equilibrium constants for the dehydrogenation 
of butenes to butadiene. Scale of abscissas is linear in 1/T. 

Curve A—1-butene(g)—1,3-butadiene(g)+ Ho(g) 

Curve B—cis-2-butene(g)—>1,3-butadiene(g)+ H2(g) 

Curve C—trans-2-butene(g)—1,3-butadiene(g)+H2(g) 

Curve D—mixed normal butenes(g)—1,3-butadiene(g) 

+H2(g) 

o—Dementjeva, Frost, and Serebriakova (nature of 
butenes unspecified). Circles with arrows are reported as 
upper limits 

e—Balandin, Zelinsky, Bogdanova, 
(1-butene, compare with curve A) 

@—Universal Oil Products (mixed butenes compare 
with curve D). 


and Sceglova 


the formula 
1 1 1 1 


Kix K wene-2 





Keis—2 Ky, 
where the K’s are the equilibrium constants for 


TABLE XIV. A comparison of calculated and experi- 
mental*® equilibrium constants for the gas phase dehydro- 
genation of n-butane to butene: n-butane—butene+hy- 
drogen. 











T Kp in atmos. 
°K trans-2-butene cis-2-butene 1-butene 
623 obs. 0.00083 0.00052 0.00045 
calc. 0.00045 0.00026 0.00015 
673 obs. 0.0039 0.0025 0.0022 
calc. 0.0025 0.0015 0.0010 
723 obs. 0.014 . 0.0087 0.0075 
calc. 0.011 0.0067 0.9048 








a F, E, Frey and W. F. Huppke, Ind. Eng. Chem, 25, 54 (1933) 
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the separate reactions above. This relation is 
readily obtained from the fact that the partial 
pressure of the mixed m-butenes is the sum of the 
partial pressures of the trans-2-butene, the cis- 
2-butene and the 1-butene. 

The published experimental data for reactions 
forming 1,3-butadiene from butenes include those 
of Dementjeva, Frost, and Serebriakova,” 


‘ Balandin, Zelinsky, Bogdanova, and Sceglova," 


TaBLE XV. Equilibrium constants for the gas phase 
dehydrogenation of the individual and mixed normal 
butenes to 1,3-butadiene: butene—1,3-butadiene+hy- 
drogen. 











ri Kp in atmos. 

a cis-2-butene érans-2-butene 1-butene mixed butenes 
298.16 1.19 K10745 3.24 K10716 1.09 X10-™ 2.48 X10716 
300.00 1.59 X10715 4.37 K10716 1.44 X10-4 3.35 K10716 
400 2.38 K10719 9.79 X10 1.04 X10-9 6.50 X10711 
500 3.61 K1077 1.84 X1077 9.77 X1077 1.08 X10- 
600 5.31 K1075 3.04 K10-5 1.01 X1074 1.62 X1075 
700 1.98 K1073 1.22 X1073 2.90 X10-3 6.00 X1074 
800 3.10 X10 2.00 X1072 3.69 X1072 9.14 K10-3 
900 0.269 0.180 0.271 7.72 X1072 

1000 1.53 1.05 1.35 0.426 
1100 6.40 4.47 5.05 1.73 
1200 21.1 14.9 15.2 5.54 
1300 57.9 41.4 38.6 14.8 
1400 138 98.6 85.3 34.4 
1500 292 212 171 71.4 








 M. I. Dementjeva, A. V. Frost, and E. K. Serebria- 
kova, ‘‘Equilibrium dehydrogenation of n-butylenes to 
butadiene -C,Hs<>C,He+Ho,” Comptes rendus (Dok- 
lady) Acad. Sci. U.R.S.S. 15, 141 (1937). 

“t A. Balandin, N. Zelinsky, O. Bogdanova, A. Sceglova, 
“Production of divinyl by catalytic dehydrogenation of 
butylene,”’ J. App. Chem. U.S.S.R. 14, 435 (1941). 


and the Universal Oil Products Company.” 
These data and the calculated equilibrium con- 
stants are shown for comparison in Fig. 4 in 
which values of logio K, are plotted against a 
scale linear in 1/T. 

It is interesting to note that the experimental 
data are in best agreement with the calculated 
equilibrium constants for the dehydrogenation 
of 1-butene to 1,3-butadiene. In the work of 
Dementjeva et al. the composition of the butene 
was not given; in the work of the Universal Oil 
Products Company it was described as a mixture 
of a- and 8-butylene (1- and 2-butenes); in the 
work of Balandin et al. it was stated to be a- 
butylene (1-butene). 

It is believed that the calculated values of the 
1,3-butadiene-butene equilibria are more reliable 
than the experimental values. Experimental de- 
terminations of equilibrium constants should be 
made by approaching an equilibrium condition 
from opposite sides as was done by McCarthy 
and Turkevich* for the 1-butene-2-butene equi- 
librium (see Section VIII-1). Apparently this 
was not done by the above-mentioned inves- 
tigators of the 1,3-butadiene-butene equilibria. 
Chemical analyses of butadiene-butene mixtures 
are difficult, and there is a question whether the 
investigators made satisfactory analyses. 


42 Universal Oil Products Company, ‘Dehydrogenation 
of aliphatic hydrocarbons,” Chem. Zentralbl. I, 4262 
(1938). 
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A molecular theory of relaxing media is presented which gives an expression for the stress in 
terms of the strain history. At any given time the strain history produces a distribution in 
internal strains which for mechanical properties can be characterized by a limited number of 
internal strain parameters. The second law of thermodynamics is used to define dissipation of 
energy at constant temperature and explicit expressions for dissipation of energy for any 
strain history are obtained. Inasmuch as relaxation during straining causes an essential reor- 
ganization of structure which is in fact the cause of dissipation, the kinetic theory of elasticity 
is extended to non-isotropic polymeric networks. A tensor expression for the stress-strain-time 


relations is thereby developed. 





ANY attempts have been made to develop 

a general theory which would describe the 
mechanical behavior of bodies which exhibit 
simultaneous elastic and viscous behavior. The 
classical relaxation theories of Maxwell, Voigt, 
and Wiechert on the one hand, and the super- 
position principle of Boltzmann on the other, 
provide a substantial mathematical framework 
for such theories, but admittedly the original 
classical ideas must be regarded as lacking a 
modern molecular theory foundation. Recent 
researches,* particularly with regard to high 
polymeric materials, have been concerned mainly 
with attempts to interpret the classical mathe- 
matical framework in terms of the new develop- 
ments of molecular structure theory and chemical 
rate theory. Although these researches have 
thrown much useful light on many aspects of 
the general problem, particularly with regard to 
the nature of the local relaxation process,’ it is 
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the purpose of this paper to give a complete 
detailed geometrical picture of the molecular 
changes that occur in the body as a whole 
during the processes of deformation which in- 
volve relaxation. 

For polymeric networks such as rubber the 
kinetic theory of elasticity® provides a useful and 
self-consistent molecular theory of the relation 
between stress and strain under various types of 
deformation in the absence of relaxation. It 
provides a point of departure for the molecular 
model of hysteretic materials that we shall here 
present, although the important ideas can no 
doubt be generalized to substances other than 
polymeric networks. 

We shall deal here with materials for which no 
permanent changes are occurring with time; in 
other words the internal processes of relaxation 
that are occurring in the material are character- 
ized by equal rates of breaking and reforming of 
the bonds which constitute the network chains 
(either along the chain or at the network junc- 
tures). For simplicity we shall treat materials 
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exhibiting only one internal relaxation time, 
although the extension to systems characterized 
by a distribution of internal relaxation times is 
obvious. 

A class of substances which approximately 
obeys the conditions enumerated above is the 
polysulfide rubbers. When held at constant 
elongation the stress decay can approximately 
be represented as: 


f=fo exp (—’'t), 


where f is the stress at time #, fo the initial stress, 
t the time, and k’ a relaxation rate. For sufficiently 
small elongations this decay law is independent 
of the elongation. 

On the other hand, if these rubbers are main- 
tained in a relaxed position and only occasionally 
and momentarily extended to a fixed elongation, 
it is found that there is no change with time of 
the stress required to attain that fixed elongation. 
In other words, there is no modulus change with 
time (except at very high temperatures). This 
does not mean that no internal processes of 
breaking of bonds occurs in the relaxed position, 
but merely means that the bonds that break are 
reforming, and that the rates of breaking and 
reforming are exactly equal. This must be so 
according to the kinetic theory of rubber elas- 
ticity which states that the modulus is propor- 
tional to the concentration of network chains. 

Since the decay of stress at constant extension 
is approximately independent of the extension 
(at least for small values of the extension), it is 
not possible to explain this phenomenon in terms 
of an increase in the rate of bond breaking with 
elongation of the network. Nor is it reasonable 
to assume that the rate of bond reformation is 
affected appreciably by the elongation of the 
sample. The only reasonable explanation for the 
decay of stress at constant extension which 
occurs in spite of the fact that bond breaking 
and bond reformation are proceeding at equal 
rates is that the reforming of bonds occurs in 
such a way they do not contribute to the stress 
at constant extension,’ whereas the breaking of 


®M. D. Stern and A. V. Tobolsky, J. Chem. Phys. 
14, 93 (1946). 

7A. V. Tobolsky I. B. Prettyman, and J. H. Dillon, J. 
App. Phys. 15, 380 (1944). 

R. D. Andrews and A. V. Tobolsky, paper on permanent 
set soon to be published. 


bonds relaxes the stress in chains that are in a 
stretched condition. It is indeed very reasonable 
to suppose that formation of bonds produces 
relaxed chains no matter what the external 
dimensions of the sample or what stresses are 
acting on the sample. 

Mathematically stated we can classify the 
chains of the polymeric network into chains of 
type 7 which were formed when the rubber 
sample was held at length /; or at shear strain oj. 
According to the kinetic theory of rubber elas- 
ticity and the assumption that chains formed at 
length /; were formed in a relaxed state, the 
relation between stress on attained section and 
deformation of the sample is given by 


f/ckhT => s,(l?/l?2—1,/l) tension, (1) 
/ckT =>, s;(o—o;) shear, (2) 


where f is the stress, k is Boltzmann’s constant, 
T the absolute temperature, / the actual macro- 
scopic length of the sample, and o the actual 
macroscopic shear of the sample. s; is the number 
of chains that were formed when the sample had 
a length /; or shear o; and c is a constant defined 
in Appendix A. For purposes of illustration we 
shall also consider the following law for tension 
(which is only an approximate law, but which 
is much simpler than Eq. (1)): 


f/ckT => s\(l/l;—1). (3) 


Equations (1)—(3) represent the instantaneous 
relations between stress and strain arising from 
the instantaneous chain length distribution. It 
is useful to introduce the following definitions 
corresponding to the quantities that appear in 
Eqs. (1)-(3). 


> 5;/l2=s9/m?; > sdi=son, (1a) 
Le Sidi=Sooe, (2a) 
LD $i/li=So/le, (3a) 


where in each case so is the number of network 
chains of all kinds per cc. 

In order to derive the general relation between 
stress-strain and time it is necessary to know the 
rate of “breaking” and the rate of forming of 
network chains. It will be assumed that the rate 
of breaking of any given set of chains s (due to 
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breaking of bonds) obeys the law The stress at any subsequent time ¢ owing to 
—(1/s)(ds/dt) =R’. the chains of the network which have remained 

; unbroken since t=0 is obtained by multiplying 

The number of chains per cc. of all types re- the right-hand side of Eqs. (4)—-(6) by exp (—’t). 
forming in the time interval d?’ is therefore sok’dt’. }{owever, a contribution to the stress at time ¢ 








At any reference time #=0 we may write also arises from chains that have been broken 
f/ckT =soL2?/mo?—no/1], (4) and reformed in the time interval dt’ (between 0 
a 3 and #) and have not subsequently relaxed. The 
f/ckT = soo —o00], (5) integral equations between stress, strain, and 
f/ckT =so{l/leo—1]. (6) time are therefore of the following form: 
f W(t) \? 1’) 
——=5,) exp (— vo[—-"|+ sok’ exp [—k’(t-t 7 ( ) -— lw. (7) 
ckT me ue) U2) 
f t 
——— = So CXD (—k’t)[o—o0 |+ Sok’ exp [ —k’(t—?’) JL o(t) —a(t’) jdt’, (8) 
ckT 0 
: (—wof- 1 ‘a L—h'( ye (9) 
—=5,exp (—h’t -- }+fs “exp [—h'(t-?’ | r. 
ckT leo _— I(t’) 


The integrand in Eqs. (7)—(9) represents the contribution to the stress of those chains which have 
formed in the interval between ?¢’ and ¢’+d?’ (namely, sok’dt’) and have still remained unbroken at 
time ¢ (namely, sok’dt’ exp [—kR’(t—?’) }). 

The integral Eqs. (7)—(9) can be reduced to the following forms: 




















f A. 
——_=—--, (10) 
csokT =m? 1 
where 
1 1 texp (k’t’) ’ 
—=exp (— H'| +f ——we . 
m? Mo” 0 /(t’) : 
and 
t 
, 
n=exp (—¥'| not f exp (k’t’)1(t’)k'dt’ 
0 = 
t 
=o—exp (—H0| oat f exp Wrho(eyear' | (11a) 
csokT™ 0 
j ; 
=o—Ge, (11b) 
CSokT 
where 
t 
oe= exp (-¥) (cot f exp WV)o(vWdt ), 
0 
f l 
=——1, (12) 
csokT 1, 
where 


1 meer), 
"eel vo(— af i 1) Wal’). 


It is to be particularly noted that Eq. (11a) which is an alternative form for expressing the integral 
equation for shear is equivalent to the equation proposed by Boltzmann in his paper on the super- 
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position principle in elastic-viscous behavior, if 
the Boltzmann recollection function be taken as 
exp [—k’(t—?’) ] and oo is zero. 

The parameters o, and /, of Eqs. (11b) and 
(12) have a particularly simple physical signifi- 
cance. As can be readily seen, they are the values 
of the shear and of the length to which the rubber 
network would return if the external stresses 
were removed. They can be interpreted as being 
the values of ‘‘internal shear’ and ‘internal 
length,” respectively. It is a consequence of the 
fact that we are dealing with a system having 
but one relaxation time that upon removing the 
external stresses there is an instantaneous change 
in deformation, but no changes with time in the 
external dimensions thereafter. The physical 
meanings of m and n defined in Eq. (10) are not 
quite so suggestive as in the case of ¢, and /,. 

Equations (10)-(12) can be converted to 
differential equations by differentiations with 
respect to time of 1/m?, n, o,, and /, followed by 
further algebraic manipulation. 





d(1/m?) 1 1 dn . 
=" (—-—); =H d-m, (13) 
dt 1? m? dt 
do,./dt=k'(s—«,), (14) 
W(i/l,)’ ft 1 
= ~=¥(---), (15) 
dt lL i, 


Differential equations corresponding to (10)- 
(12) can all be obtained by differentiation of 
these equations with respect to time followed by 
substitution of (13)—(15). Inasmuch as the differ- 
ential equation corresponding to (10) is a lengthy 
second-order equation we shall not write it here. 
The equations corresponding to (11) and (12) 
are (if we let G=csokT) 


do idf k’ 
tie el, (16) 
dt Gdt G 


L 1 . h’ 
(Gate)/G4) 
L- G a 
It is to be particularly noted that Eq. (16) is 
identical to the equation proposed by Maxwell 
to explain elastic viscous behavior. However, 
Eq. (17), based on a slightly different (non-linear) 
stress-strain law, gives an entirely different 





equation. The creep behavior under constant 
load is strikingly different, for example, for these 
two cases, Eq. (16) giving a linear creep with 
time whereas Eq. (17) yields an exponential 
creep. 


DISSIPATION OF ENERGY 


In thermodynamics the properties of a system 
in equilibrium are specified by certain param- 
eters, temperature, and density, for instance, 
which are said to determine the state of the 
system. When the state parameters are given, 
all macroscopic physical quantities are thereby 
fixed. The systems which we are considering also 
have parameters which may be called state 
parameters because all macroscopic physical 
quantities are determined when they are given. 
For the case of the hysteretic material under- 
going shearing strain they are the actual shear 
strain o, and o, the shear strain at which the 
shearing force would be zero. In the case of the 
material undergoing a tensile strain the state 
parameters are /, the actual length of the sample, 
and the quantities m and n. 

The states determined by these parameters are 
not thermodynamic states because they are not 
in general equilibrium states. However, it will 
be profitable to apply some of the ideas of 
thermodynamics to them. In particular it will 
be profitable to see what the consequences of the 
second law of thermodynamics are for transfor- 
mations between these states. The transforma- 
tions which are being considered all occur at 
constant temperature. The second law of thermo- 
dynamics states that if a body undergoes a 
cyclic process during which it exchanges heat 
with a single heat reservoir at a fixed temperature 
the net work done on it must be positive. The 
transformations usually considered in the appli- 
cations of this principle are processes in which 
the initial and final state are the same thermo- 
dynamic equilibrium state. In this discussion 
this principle will be assumed to apply to cyclic 
processes in which the initial and final state are 
the same state, where states are defined by their 
values of o and o,, whether these are thermo- 
dynamically equilibrium or non-equilibrium 
states. 

The internal energy of the material being 
considered is a function of the temperature alone. 
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according to the assumptions of the kinetic 
theory of elasticity, and does not change during 
any process which occurs at constant tempera- 
ture. If an amount of work dW is done on the 
system, an amount of heat dQ equal to dW must 
be given up by the material to the surroundings. 
This represents an increase of entropy of the 
surroundings by an amount dW/T. The net 
effect of the cyclic process is that work has been 
done, the surroundings have been heated, and 
the entropy of the surroundings increased, while 
no change has occurred in the material itself. 
Moreover it is not possible by any other cyclic 
process to get back any of the energy which has 
been delivered to the surroundings. It is natural, 
therefore, to define the heat energy delivered to 
the surroundings in a cyclic process as the energy 
dissipated in the process. 

It is somewhat more difficult to define the 
dissipation of energy in non-cyclic processes for 
some of the energy delivered to the surroundings 
during the process can be recovered in another 
non-cyclic process. However, it is true that if a 
non-cyclic process is completed so as to make a 
cyclic process the total work done on the body 
and the total heat delivered to the surroundings 
are positive. In symbols, if W(AB) is the work 
done on the material in a process whose initial 
state is A and whose final state is B and W(BA) 
is the work done in some cycle-completing proc- 
ess W(AB)+W(BA) 20, —W(AB)< W(BA). If 
we consider a definite process AB and all possible 
cycles which will complete it, BA, it is seen that 
there is one completing cycle (BA), for which 
W(BA), is the smallest, for the work done in 
any one of these processes is always more than 
or equal —W(AB). The dissipation of energy in 
the non-cyclic process AB can then be defined 
as the work done in AB plus the smallest amount 
of work which can take the body from B back 
to A again. In symbols 


Dissipation = W(AB)+ W(BA),. 


If the process (AB) has once been allowed to 
occur, at least this much heating of the surround- 
ings must take place before the state A can be 
re-established. If it should happen that there is 
a process which goes from A to B (not neces- 
sarily the one which we have been considering) 
and another which goes from B to A, say (BA), 
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such that W(AB)+W(BA)=0, then W(BA) 
= W(BA),, the smallest amount of work which 
can be done on the system to re-establish A. In 
particular if there is a reversible process between 
A and B, —W(AB),= W(BA),=W(BA)-. 

We shall use the conception of the second law 
of thermodynamics which has just been discussed 
to develop simple expressions for the dissipation 
in any process occurring in media which obey 
the equations developed here. We consider first 
the work done in shearing a one centimeter cube 
of hysteretic material. In the shearing process 
two parallel faces of the cube are moved parallel 
to each other. The area of these faces does not 
change nor does their distance apart. The tan- 
gential stress on one of the faces is given by 


f 
CSokT 





=T—d¢. (11b) 


The tangential distance the face moves is o so 
that the work done between the times 0 and ¢ is 


é do 
Wook [ (o—o,)—dt. 
0 dt 


By integration by parts this expression can be 
transformed to 


Wo? ‘ ‘ de, 
=“ —oe.| +f o—dt ; 
CSokT 2 0 0 dt 





by adding the identity 
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we have that the work done in any process is 
given by the expression 


W (e-¢,)? 
csokT 





|+ f k'(o—«,)%dt. (18) 


This expression shows that the ideas developed 
here are consistent with the second law of 
thermodynamics for if (0) = a(t) and ¢.(0) =c.(t), 
the work done is positive as required by the 
second law. The dissipation is given by 


t 
csukT J k’(o—o,) dt. 
0 


This expression for the dissipation holds good 
even if the process is not cyclic. If (01, 0-1) are 
the state parameters for the initial state of such 
a process and (¢2, ¢.2) those for the final state, 
there is a reversible process which goes from 
(o2, Ge2) to (1, 1). This process can be repre- 
sented on a two-dimensional diagram in which 
@ is plotted against o, as shown in Fig. 1. Physi- 
cally the process is this. The cube of material is 
allowed to snap back from its actual position o2 
to its equilibrium position at o.2. If the mass of 
the material is negligible, this can take place in 
such a short time that o, does not change during 
the process. The work done in this process is 


CSokT 





(2 — Oe2) 2. 


By exerting a very slight shearing force o and o,, 
which are now both equal, can be made to change 
until they are both equal to o.1. This process 
takes a very long time. No work is done in it. 
Finally o can be changed from ¢,1 to o1 by doing 
an amount of work csokT(o1—o-¢:)?/2. The 
smallest amount of work which can be done to 
re-establish (01, o-1) from (o2, oe2) is then 


which when added to the actual work done in 
going from 1 to 2 gives for the dissipation 


t 
csukT J k'(o—o,)*dt. 
0 


The dissipation per unit time per unit volume is 


csokTk'(o—o,)? and is a function of the state 
parameters o and o,. It is zero if k’ is zero, that 
is to say if there is no breaking of bonds. The 
dissipation is also zero if the process is carried 
out so slowly that o=oe, throughout. 

This method of deriving the rate of dissipation 
can be applied to the other cases of dissipative 
strain. A more general point of view can be 
obtained by considering the problem as one in 
Pfaffian differential expressions. We will apply 
this method to the simplified equation for tensile 
stress. The work done in a small change in 
length of the sample is 


dw l 
-(--1)a, (19) 
csokT A i. 





where A is the area of the cross section of the 
sample. This is a Pfaffian differential expression 
in which the coefficient of di, is zero. In a time 
dt according to Eq. (15) 7, changes in such a way 


that 

dl, 1 1 

—+#(-——)at=0. (20) 
L,? ; & 


This Pfaffian expression is the equation of motion 
of 1,. The problem of demonstrating consistency 
with the second law and determining the rate of 
dissipation can be considered to be equivalent to 
the following : Find a function of / and /,, X (J, /.), 
such that the Pfaffian expression 


l dl, i 1 
(-- 1 )at-x—4-x¥' (-—— )at (21) 
le l,? L |, 


can be broken up into two parts one of which is 
the differential of a function of /.and /, and the 
other a positive function D of / and /, times the 
differential of the time. Then 


dw 
CSokTA 


W t t 
-F| +f Dae 
CSokT A 0 0 


and the confirmation of the second law can be 
carried through as before. 

In order that ((//1,) —1)dl+-Xdl,./l?2 should be 
a perfect differential, X must satisfy the equation 


=dF+Ddt, 





and 
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X = —(P/2)+g(l-) where g(l.) is some function 
of J.. g must be chosen so that D is a positive 
function of / and 7,. A possible choice is g(/.) 
=)2/2. 

Then 





l 12-2? 
aF=(- ~ 1 )al+ dey 
L, 2 


dF is seen to be the differential of the function 
F=(1,/2)(l/l.—1)? plus a constant, and D is 
positive (or zero) for all / and /,. 

Any choice for g(/.) will make the differential 
((1/l.) —1)dl+-Xdl,/l2 a perfect differential, but 
it is easy to see that only the choice g(l.) =12/2 
will make the coefficient of dt a function of / and 
1, which is always positive. If g(/.) is written as 
(12/2)+¢(l.) where g(/.) is another arbitrary 
function of /,, D must be replaced by 


[(1/2)(P-1?— (Le) ((1/Le) — (1/1) JR’, 


and if g(l,) is different from zero for any value 
of /, the second parenthesis can always be made 
to have a different sign from the third by a 
suitable choice of 1. The only choice therefore 
which will make D always positive is ¢(/.) =0. 
It is seen that this way of expressing the 
differential of work is unique and therefore a 
definite work function and rate of dissipation 
function are determined by the integral equation 
for stress. This case is analogous to the preceding 
one in that a very simple reversible process can 
be found between any two states (J1, 11) (ls, ¢2). 
Moreover the work done in this reversible 
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process is simply 
l final state 
csokT Al,{ ——1 }. 
Ll, initial state 


P—f2y 71 1 
esskTAn'(——)(——-) 
2 I 1 


actually is the rate of dissipation of energy. 
The case of the unsimplified equation for 
tension is somewhat more complicated because 
there are three rather than two parameters as 
before. The stress is given by Eq. (10). 
f ye 8 


csekT m2 | 


so that 





The area of the material changes in such a way 
that the volume remains constant so that if the 
volume is V the area is V// and the force is 
given by the equation 


force ( [? ~) 
csokTV I m? ld 
The differential of work is 
dw isi n 
- --(—- ja 
csokRTV 1l\m? 1 


The differential expressions for the equations of 
motion are 








(22) 


dm 1 1 
o=2——+8'(-—— at, (23) 
ms l? m? 
and 
0= —dn+k'(l—n)dt. (24) 


By using the multipliers X = m?/2—F/2 for (23) and Y=1/n—1/I for (24) the differential for the 


work can be put into the form 


dw 1s? n 1s? 1/n 
--(—-")ai-—( -1)dm+—(“~1 
csokRTV 1l\m? 1 m\m? n\ 1 








)an 


1 1 . 4 
+{a(m—1)(—— )+(--=) nn |an (25) 
1? m? n 1 





In very much the same way as before it can be proved that this is the only expression for the work 
for which D is always positive. The work done in any process is then 
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This is enough to show that the model is con- 
sistent with the second law of thermodynamics 
for tensile strain for in a cyclic process the first 
term certainly vanishes and the work done is 
positive. However, there is no simple reversible 
process between any state (J1, m1, 1) and any 
other (lz, m2, 22). For such a path the integral 
term would be zero. This could only happen if 
the process takes a very short time or if l=m=n. 
If the material is allowed to snap back so that 
(?/m?) —(n/l) =0 there is ne contribution to the 
integral because the process occurs so rapidly 
that no relaxation occurs. m and n, however, are 
still different and cannot be made to be the 
same except in a slow process in which case there 
is a contribution to the integral. 

There is no way of accomplishing a reversible 
process between any two states by simply ad- 
justing the length of the sample. Nevertheless, 
it is conceivable that there is some system which 
would interact differently with the rubbery 
material when / and m and m are the same and 
when they are different. In this case the two 
systems together might be in thermodynamic 
equilibrium, even if 7 and m and nm are not 
equal. By adjusting the external parameters of 
the auxiliary system a reversible process could 
then be found between any two states. In the 
differential of Fthe quantity — (1/m)[(2?/m?) —1] 
would be interpreted as the generalized force cor- 
responding to a change in m, and (1/n)[(n/l) —1] 
the generalized force corresponding to a change in 
n, and inasmuch as the work done in a reversible 
path would then be given by the first term of 
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Eq. (26), the dissipation would actually be given 
by csokT f Ddt as in the previous cases. 

It is useful to consider a few examples of 
processes in which dissipation occurs for the case 
of shearing strain. We consider first creep at 
constant load in which a constant tangential 
stress is applied to the surface of a parallelopiped 
of hysteretic substance. The material is assumed 
to be originally cubical and to have been left 
alone for a long time so that the o and go, are 
the same in the unstretched position. The 
cubical shape will be taken as zero strain. When 
the load is first applied, the cube immediately 
assumes a sheared position which slowly changes 
as time goes on. If f is the constant shearing 
stress the first position assumed after the stress 
is applied is given by the equation 


f/(csokT) =a, 


where gp is the first sheared position. Thereafter 
the shear strain is determined by the equation 


f/csokT =o —0,=0. 


A differential equation for o can be obtained 
by substitution f=constant in Eq. (16). This 
becomes 

da/dt=k'so, 


o=(1+k’'t)oo, 
o.-=o—a9=k' too. 


On a (oc, ¢,.) diagram the process can be plotted 
as in Fig. 2. At time ¢ the work done on the 
material is simply the work done in the first 
quick extension csokTo.?/2 plus the work done in 
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actual creeping which is f(¢—«oo) =csokTk’to,?. 
The total work is csokT((1/2) +k’t)oo?. The dissi- 
pation function is csokTk’o,?. (Note that for this 
case the rate of dissipation is constant.) The 
dissipation is csokTk’too. The work function 
change is csokToo?/2. It is seen that all the work 
done in creeping is dissipated while all the work 
done in the first quick extension goes into 
changing the work function. 

The second dissipative process which will be 
considered is stress relaxation at constant shear 
strain. In this process the cube of material is 
sheared to a shear strain o) and kept there. The 
tangential force is originally csokToo and is 
csokT(o—«o-.) thereafter. By its definition 


t 
omen (—¥D f exp (k't’)o(t’)dt’, 
0 


t 
=o exp (ep f exp (k’t’)dt’, 
0 


=¢9(1—exp [—&’t]}). 


a. is originally zero but increases. The force 
decreases until after a very long time o becomes 
equal to o, and the force becomes zero. The 
process can be represented on a (¢, o-) diagram 
as in Fig. 3. The ultimate outcome of the process 
is that the cube takes up the stressed position oo 
permanently. The work function change is zero 
because o equals ¢, in the initial and final states. 
Work is done on the rubber in the first quick 
extension of amount cspkToo?/2. No work is 
done in the actual stress relaxation. The dissi- 
pation function is csokTk’o? exp (—2k’t) and the 
total dissipation csokTo,?/2 just equal to the 
work done in the first quick extension. This is 
an example of a process in which dissipation 
goes on although no work is being done, and 
consequently no heat is being delivered to the 
surroundings. 

It is of interest to note that for the case of 
actual tension, it is possible to have a condition 
for which the sample is under no external stress 
(and hence maintains a constant value of length 
since no recovery is exhibited by a medium 
with a single relaxation time) which nevertheless 
is dissipating energy according to Eq. (26) 
because the values of m and 7 are not equal. 
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The Viscosity Law 


The limiting case of elastic-viscous behavior, 
when the relaxation rate of the medium is very 
rapid compared to the change of rate of strain 
but the times of observation are long compared 
to the relaxation time, can be studied (in the 
case of shear) by substituting the Taylor expan- 
sion of 


a(t) =0(0)+0'(0)t+30"(0)H2+--- 


in the integral Eq. (11). If terms containing 
exp (—k’t) and o’’/(k’)’, etc., are regarded as 
approaching zero, the equation reduces to the 
Newtonian expression 


f=(G/k')o’. (27) 


This relation would, of course, also be deduced 
from the Maxwell equation. The constant G/k’ 
must be interpreted as the coefficient of viscosity 
of the medium. 

If the same procedure (expansion into a 
Taylor series and substitution into the integral 
Eq. (10)) be applied to actual tension, the 
following is obtained: 


f= or ri (28) 


which can be interpreted as the viscosity law for 
tension. 

It is interesting to note that we have derived 
these expressions for viscosity on the basis that 
the force has no effect in altering the rate of 
relaxation. If couched in terms of Eyring’s formu- 
lation of relaxation processes‘ we are dealing here 
with a unit process for which the free energy 
barrier is very: asymmetric, such that the local 
force does not increase the rate of traversing the 
barrier in the forward direction (to the relaxed 
state) and the rate of jumping the barrier in the 
backward direction toward the strained state is 
negligibly slow. In the theory developed here 
the dissipation of the medium (which is mani- 
fested as viscosity) cannot be explained in terms 
of a single representative unit process of relaxa- 
tion, but arises from the fact that the internal 
distribution of strain lags behind the externally 
imposed strain. 

At this point it might be asked to what extent 
these results are applicable to all types of ma- 
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terials, and to what extent they are applicable 
only to polymeric elastic substances for which 
the kinetic theory of elasticity applies. The 
answer is that the theory developed here applies 
to media for which the following hold true: 
(1) The- external stress can be computed as a 
sum of stresses due to the straining of individual 
elements; (2) the relaxation process is such that 
the number of structural elements does not 
change with time; (3) each effective relaxation 
process consists in the breaking of the con- 
straints on a strained element and its subsequent 
reformation in an unstrained condition. If these 
conditions are met, and the stress-strain law for 
the purely elastic behavior of the medium is 
known, the integral equation relating stress- 
strain and time can be written down immediately 
in the manner shown previously. 

The further particular assumptions of this 
paper have been: (1) There is only one kind of 
relaxation process, and one single relaxation rate ; 
(2) the relaxation rate is independent of the 
stresses on the elements. In the case of media 
exhibiting several relaxation processes character- 
ized by a distribution in relaxation rates, the 
generalization of our equations is very simple. 
The result can be stated as follows: The ex- 
ternal stress is a sum of stresses which can each 
be attributed to a definite relaxation process and 
for each of which the simple integral Eqs. (8) 
can be written. This leads to a formulation 
identical to the Maxwell-Wiechert theory in the 
case of a linear stress-strain relation. 

The case that the relaxation rate depends on 
the stress on the element according to the law 


hk! = ky'eFf 
can be accounted for formally in the integral 


equations. In the case of relaxation of stress at 
constant length this leads to the formula 


f f 
log ko’t= —Bfo—+1n (in =), 
fo f 
where fo is the initial stress. 
APPENDIX A 


The kinetic theory of rubber elasticity as 
developed by a number of authors® is based on 
the following assumptions: (a) The structure of 


the materials for which this theory applies are 
assumed to be polymeric networks with rela- 
tively fixed juncture points and with compara- 
tively free motion of the portions of the mo- 
lecular chains connecting network junctures, (b) 
the internal energy and the volume of the ma- 
terial do not change with strain (no compressi- 
bility), (c) the stresses arising from a given strain 
are due to a change of entropy alone. The 
entropy is calculated by means of Boltzmann’s 
relation on the assumption that the number of 
configurations available to a network chain with 
a given vector distance r, between its ends is 
exp (@n—B,7fn”), where a, and 6, are constants 
which depend on the structure of the chain. 

It is the purpose of this section to develop the 
kinetic theory of elasticity for general three- 
dimensional strains (Appendix A) and to extend 
this theory to take into consideration changes 
with time that occur due to the breaking and 
reforming of network chains (Appendix B). 

The general expression for the free energy F 
per unit volume of a polymer network obeying 
conditions (a) and (c) is 


F=Fo—kT Son QntkT Lin Br2tn?, (At) 


where F» is a function of the volume and temper- 
ature, and the summation is over all chains in 
the unit volume. 

The equilibrium distances between network 
junctures in the unstretched state can be ob- 
tained by minimizing F with respect to all 
variations of the r, which do not break any of 
the junctions of the network. It is clear that 6 
alone is insufficient to specify these distances and 
that a detailed knowledge of Fo is required. It 
will become apparent in the development below, 
however, that it is merely necessary to know the 
mean value of 8,’r,2 in the unstretched state (if 
the unstretched state is isotropic) to obtain the 
equation of state. 

According to assumption (b) Fo and the vol- 
ume do not change under any allowable strain. 
The part of the free energy which changes in a 
straining process is therefore 


F=kT ¥.n Ba?t pn”, (A2) 


where r,’ is the new chain vector in the strained 
state. 
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Assuming that all the vectors r are transformed 
by the same tensor T which is therefore the same 
as the macroscopic strain tensor, we can write 
r'=Tr and r?=trace (T(rr)7T) where (rr) is the 
tensor which transforms a vector x into r(r-x), 
and T* is the transposed tensor of T. According 
to formula (A2), the free energy per unit vol- 
ume in the strained state F’=kT)., 8,” traces 
(T(r,r,) 7) where the summation is over all chains 
in the unit volume of the strained state. If we 
classify all chains in the unstrained states into 
groups which have the same 8 and r and if s; is 
number of chains in the /th group this sum can 
be written 


Dez $817(Titi) \ 
SsokT trace (7(———) r), (A3a) 


So 


sokT trace (TET), (A3b) 


because the number of chains per unit volume of 
any type is the same in the strained and un- 
strained state. The tensor £ is a local property of 
the unstrained state and depends on a detailed 
knowledge of the equilibrium chain lengths and 
6’s in the unstrained state. It is always sym- 
metric and positive definite. If the unstrained 
state is isotropic as would presumably be the 
case in a sample of material which had not been 
stretched for a long time, all off-diagonal terms 
of would be zero and all diagonal terms would 
be equal. would then be 


1d s6Pr? 
3 So 2 


where I is the unit tensor. The free energy per 
unit volume of the strained state would then be 


P 
=sokTtr(TT). (A4) 
This formula is precisely the same as Treloar’s® 
except that there appears an additional param- 
eter c which is not necessarily one as it is in 
Treloar’s formulation. 

In order to calculate the stresses in a substance 
which has a strain energy function given by 
(A3b) we consider the change in the strain energy 


* Throughout this paper the T’s are tensors and would 
normally be boldface but for difficulties of type setting. 
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due to a small additional strain and equate this 
to the work done, taking into account the fact 
that the small additional strain must not change 
the volume. If the positions of material points 
of the substance in the strained state are denoted 
by R’ and in the unstrained state by R the 
strain tensor -T=R’Vr where R’Vr is the 
tensor whose zjth component is dx;//dx; and 
where {x,/} are the components of R’ and {x;} 
those of R. If each material point in the strained 
state is made to undergo a small displacement 
6R’ the additional strain 6T=6R’Vr. This ex- 
pression for the additional displacement is some- 
what inconvenient because the gradient is taken 
with respect to the unstrained position of the 
displaced point. Expressed in terms of gradients 
with respect to the strained positions 


6T = (5R’VR’)T. 


If F is the free energy of a volume V of the 
material, 








= [ wctetav, 
SokT V 
dF . 
=2f w(eteT)av (A5) 
SokT V 
=? f ir((6R’VR’)TET dV. 
V 
By means of the identity 
div [(TéT)6R’]=¢tr[(6R’Vr’)TET | 
+[div (TET) ]-6R’, (A6) 


the expression for dF can be transformed to 


dF 





= f aS-[(Te7)6R’] 
2sokT Ss 
— [ (oR’-Caiv (TéT) ]}dV. (A7) 
V 


The first term is a surface integral over the 
surface S of the volume V. The change in work 
function must be equal to the work done by the 
external forces, which is 


i) [6R’-(PdS) ], 


where P is the stress tensor, since (PdS) is the 
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force and 6R’ is the displacement. This is equal to 
S sdS-[PéR’] because P is symmetric. 





fas-p-on’= [as-(reT)-aR’ 
2sokT Ss Ss 


- i) SR’ -div (TET)dV. (A8) 


J 


This equation is to be satisfied for all 6R’ which 
are consistent with the constancy of volume. 
The condition on 6R’ is that div 6R’ is always 
zero. If this condition is multiplied by a Lagrange 
multiplier g and integrated over the volume of 
the material and added to the above equation, 
the resultant equation must be satisfied by all 
variations 6R’. 
The equation then is 


1 
2sokT 





f4s-(eer']= f a8-[(teT)R) 
1 f dS- esR’— f 6R’-div (TET) dV 
Ss V 


~ f 6R’-grad gdV, (A9) 
: 


where {¢ div 6R’dV has been transformed by 
an obvious vector identity. This equation implies 
that the stress tensor is 


P=2s.kT(TET) + ¢'l, (A10) 


where I is the unit tensor and ¢’ is an arbitrary 
scalar function, which is to be determined from 
the boundary conditions of the problem. For the 
case that the unstrained state is isotropic the é 
tensor is (c/2)I and the stress tensor is 


P=cskT(TT)+¢'l. 


For tensile strain the tensor T has the matrix 
form 


L/l, 0 0 
T=| 0 (.,/D) 0 : 
0 0 (l/lu)} 


where / is the length of the sample and /, the 
length of the sample in the unstrained state. 


The stress tensor has the matrix form 
P/l,? 0 0 ¢g 0 0 

P=csokT | 0 L/l O|;+/0 ¢ O 
0 0 L/l 0 


Under actual experimental conditions the ma- 
terial would be in the form of a rectangular 
parallelopiped with the direction of stretching 
perpendicular to one of the faces. g’ is determined 
by the condition that the stress on all faces 
parallel to the stretching direction is zero. ¢’ 
must then be —csokT1,,/1 and the stress tensor is 


?/l,2—1l,/l 0 0 
P=csokT 0 0 OO}. (A111) 
\ 0 0 0 


0 ¢g) 


It is seen that two samples of rubbery material 
with the same concentration of chains can differ 
with respect to their elastic properties only 
insofar as the tensors é differ. If both samples 
are isotropic in the unstrained state, they can 
differ only insofar as the constants c differ. 


APPENDIX B 


With the aid of the above discussion of the 
kinetic theory of rubber elasticity for arbitrary 
strains it will be very easy to extend the ideas 
developed in the main part of this paper to the 
most general three-dimensional strain history. 
It will also be possible to give a more detailed 
account of the molecular state of a hysteretic 
material which has undergone a strain history. 

The chain vectors in a sample of rubbery 
material which has had no force act on it for a 
long time will be distributed isotropically and 
the tensor é will be (c/2)I (c is a constant which 
depends on the chemical nature of the material). 
When the material is strained the distribution of 
chain vectors becomes anisotropic. However, the 
processes which go on in the material tend to 
make the distribution of chain vectors isotropic. 
Any chains which have broken down in a small 
time interval dt’, then, will reform again iso- 
tropically, and with the tensor £ equal to (c/2)I. 
The reformed chains are attached to the re- 
maining unbroken structure and when this struc- 
ture is strained further the newly precipitated 
chains will again become anisotropic. 
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The strain tensor which determines the contri- 
bution to the stress tensor of a given group of 
chains is the relative strain from the time the 
group was formed to the actual time. If T(¢) is 
the strain tensor calculated from some definitely 
chosen state as the unstrained state, the relative 
strain tensor T,(t, t’) is T(#)T-(t’), where ?’ is 
the time when the chains were formed and ¢ the 
actual time. This relative strain is simply the 
further strain which must be applied to the 
situation at time ?¢’ to get the situation at time ¢. 

The arbitrary hydrostatic pressure in the 
stress tensor does no work since the volume does 
not change. It will therefore be neglected in the 
discussion of dissipation. The contribution to 
the stress of a group of chains s; per unit volume 
formed at time #;’ is simply 

cs:kT(T,(t, t’)T+(t, ti’)), 
eskT(T()T Ut) Tt )T()). 


or 


The contribution of all chains formed before 
t=0 to the stress at t=0 is simply 


| ( (Lis. THW)TW)) 
ckT sof T —1). 


So 





The summation is the analog of the parameters 
Mo, No in the general theory. It will be called to~ 
because it is symmetric positive definite tensor 
and has a reciprocal. At time ¢=¢ the stress will be 


esokT| exp (—k’t)(Tro?T) + J exp [—R’(t—t’)] 


XBT) T-W)TO)Ha | (B1) 


where as before the first term represents the 
contribution of chains which existed before t=0 
and the second term the contribution of chains 


which were created thereafter. 
With the definition 


t 
7? =exp (—¥) (124 f exp (k’t’) 
0 


xT) T)Hal), (B2) 
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the stress becomes 


csokT(T rT). (B3) 


This is precisely the force which would be 
derived from the strain free energy function 


c " 
~sokTtr(Tr27). (B4) 
2 


The work done in a small change of strain dT is 


9 


simply the change in this function keeping 7~* 
constant 


C . _ 
dW =—sokTtr((dTr?T)+(Tr-dT)). (BS) 
2 


The equation of motion of 7~? in differential 
form is 


dr? —k'(T-T!"' — rd. 
If this equation is multiplied by 


csokT 


(7?— TT), 





and the trace of the product added to the above 
equation, the resulting expression has the form 


2dW 7 
= d(tr(Tr-*T) —log determinant (7~?)) 





cSokT 


+k trl (TT —7)(r2-TT-) jdt, (B6) 


where as before the first term is a perfect differ- 
ential and the second a positive function times 
the differential of the time. Equation (B1) can 
be reduced to the differential equation 


1 dP P dT 
aes | a 
G dt G dt G G dt 


3k’ 
“F \t (B7) 
trP-! 








1 ee 


G trP-! 
(The rate of change of the arbitrary hydrostatic 
pressure is determined from the condition that 
the determinant of T is a constant.) 





33) 


be 


34) 


T is 


9 


B5) 


tial 


Ove 


B6) 
fer- 


nes 
can 


B7) 


atic 
hat 








THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 14, 


NUMBER 2 FEBRUARY, 1946 


Stress-Time-Temperature Relations in Polysulfide Rubbers 


M. D. STERN AND A. V. TOBOLSKY 
Frick Chemical Laboratory, Princeton, New Jersey 


(Received November 14, 1945) e 


Polysulfide rubbers of various internal structures have 
been investigated by measurements of continuous and 
intermittent relaxation of stress and by creep under 
constant load at temperatures between 35°C and 120°C. 
Continuous stress relaxation measurements indicate that 
these rubbers approximately obey the simple Maxwellian 
law of relaxation of stress, which indicates that one definite 
type of bond in the network structure is responsible for 
stress decay. The activation energy for the relaxation 
process in each of the polysulfide rubbers is nearly the 
same, indicating that the same type of bond is responsible 


for the relaxation behavior of all the polysulfides investi- 
gated. In contrast to the hydrocarbon rubbers, oxygen is 
not the cause of high temperature relaxation in polysulfide 
rubbers, nor does heating in air at moderate temperatures 
for times comparable to the relaxation time produce 
changes in physical properties as determined by modulus 
or by appearance of the samples. Several possibilities 
regarding the mechanism of the relaxation process and 
the type of bond involved are considered in the light of 
the experimental results. 





EASUREMENT of relaxation of stress for 

various hydrocarbon rubbers has been 

used to investigate the relationships between the 

physical and structural-chemical characteristics 

of these materials.1 The present investigation is 

an extension of these studies to the polysulfide 
rubbers. 

The structures of rubber-like materials have 
been considered to be three-dimensional networks 
constructed from long chain molecules linked at 
infrequent points. The point where two chains 
are linked is referred to as a network juncture, 
and the portion of the long polymer molecule 
between two junctures is referred to as a network 
chain. In the unstrained state the network chains 
tend to assume randomly coiled configurations, 
the degree of randomness being limited by steric 
considerations and by interactions with neigh- 
boring chains. In the stretched state the network 
chains of the rubber are forcibly uncoiled, and 
the stress is calculated from probability con- 
siderations in the kinetic theory of elasticity? by 
the equation 


peuf({)} 


where f=stress on attained section, s is the 
number of strained elements (network chains) of 


1A. V. Tobolsky, I. B. Prettyman, and J. H. Dillon, J. 
App. Phys. 15, 380 (1944); A. V. Tobolsky and R. D. 
Andrews, J. Chem. Phys. 13, 3 (1945). 

2E. Guth and H. M. James, J. Chem. Phys. 11, 455 
(1943); F. T. Wall, J. Chem. Phys. 10, 132 (1942); 10, 485 
(1942); L. R. G. Treloar, Trans. Faraday Soc. 40, 59 
(1944); P. J. Flory, Chem. Rev. 35, 51 (1944). 
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the structure per unit volume, kT is the thermal 
energy, and //l, is the ratio of stretched to 
unstretched length. 

The network structure must be regarded as a 
dynamic structure. Strained elements tend to 
relax their strain because of slippage or rupture 
at the network junctures or by scission or rupture 
along the chains. On the other hand, new struc- 
tural elements (network chains) are always being 
formed by reformation of bonds or by the 
formation of new bonds between the chains. 

In previous papers! it was argued that bond 
reformation or cross-linking occurs in general in 
such a way that the network elements thereby 
formed are relaxed at the moment of formation 
and do not contribute to the stress in any con- 
figuration of the rubber at which they are 
formed unless the sample is further strained. 
For this reason studies of the decay of stress at 
constant extension isolate the various processes 
of interchain slippage or rupture at juncture 
points and/or bond scission along the network 
chains. On the other hand, occasional and mo- 
mentary measurements of the stress required to 
attain a definite extension, the sample being 
maintained in an unstrained state between meas- 
urements, provide a means of measuring the 
net rate of bond breakage plus bond reformation 
through Eq. (1). These experiments are referred 
to as continuous relaxation of stress measure- 
ments and intermittent relaxation of stress meas- 
urements, respectively. The techniques of these 
measurements have been previously described. 
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Fic. 1. Continuous and intermittent stress relaxation of 
Butyl rubber: 130°C, 50 percent elongation. 0 continuous, 
@ intermittent. 
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Fic. 2. Continuous and intermittent stress relaxation of 
Butaprene NM: 130°C, 50 percent elongation. © continu- 
ous, @ intermittent. 


Figures 1 and 2 show intermittent and con-— 


tinuous relaxation of stress data at 130°C for 
Butyl (polyisobutylene-isoprene) rubber and for 
Butaprene NM (butadiene-acrylonitrile) rubber, 
respectively, the data being. plotted as stress 
divided by initial stress against logarithmic time. 
These gum vulcanizates have been described in 
a previous paper.! It is apparent from these data 
that processes of bond scission and bond forma- 
tion (or cross linking) are occurring simultane- 
ously in both of these vulcanizates. In the case 
of Butyl rubber, scission is more rapid since the 
intermittent curve shows stress decay, whereas 
in Butaprene NM cross-linking predominates be- 
cause the intermittent curve shows a stress rise. 
Both scission and cross-linking are caused by the 
presence of molecular oxygen. 

Figure 3 shows similar data for H-11 at 60°C. 
It is apparent that, inasmuch as the intermittent 
curve is perfectly flat, the rates of bond breaking 
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Fic. 3. Continuous and intermittent stress relaxation of 
H-11: 60°C, 50 percent elongation. © continuous, @ 
intermittent. 
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Fic. 4. Continuous stress relaxation of H-11 at various 
elongations, and plot of exp (—k’t) at 60°C. ------ exp (—’t), 
e 10 percent elongation, © 30 percent, @ 50 percent. 


and bond formation are exactly equal. This is 
generally true for all the polysulfide rubbers 
studied over a wide range of temperature. 
Furthermore, the bond breaking reaction appears 
independent of the presence of molecular oxygen 
from experiments run in high vacuum. From 
Fig. 3 it is apparent that heating in air for times 
comparable to the relaxation time produces no 
change in modulus, nor are any other physical 
or chemical changes observable on visual and 
manual inspection. Finally, the initial stress 
required to attain 50 percent elongation was 
found not to vary appreciably over a wide 
temperature range. For example, for H-11, the 
initial stress required to attain 50 percent 
elongation over a range of temperatures from 30 
to 100°C did not vary by more than 1 percent, 
which is within the limits of experimental error. 

Figure 4 shows the continuous stress relaxation 
curves of H-11 at 10 percent, 30 percent, and 
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50 percent extension and at 60°C with stress 
divided by initial stress being plotted against 
logarithmic time. A dashed curve representing 
a plot of exp (—k’t) versus log time is also shown 
where k’ has been chosen to give the best fit of 
the experimental relaxation data. It is apparent 
that the effect of elongation on these relaxation 
curves is very small in this range of elongation, 
and that the experimental curves closely obey 
the law 


f/fo=e*™*. (2) 


It follows from Eq. (1) that the rate of stress 
relaxation caused by slippage at the network 
junctures or by rupture of bonds along the chain 
is given by 


—(1/s)(ds/dt) =k’, (3) 


where k’ may contain, as a multiplicative factor, 
the number of bonds along the chain susceptible 
to scission. This simple stress decay law indicates 
that a single type of bond (with one single rate 
constant) is responsible for relaxation. 


CREEP UNDER CONSTANT LOAD 


The behavior of H-11 (which is representative 
of all the polysulfide rubbers studied) as exempli- 
fied in Figs. 3 and 4 would appear to indicate 
that we were dealing with an almost perfect 
Maxwellian body which obeyed Maxwell’s law 


ds/dt=(1/G)(df/dt) + (k’f/G), (4) 


where ds/dt is rate of strain, f is stress, G is an 
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Fic. 5. Creep curves of H-7 at 60°C. ]=extended length 
at time ¢, /,=unextended length, Jo=initial extended 
length: k’ is obtained from stress relaxation data at 60°C. 
Logio vs. k’t, O lo=25.1 percent, @ 1o=16.9 percent, 
@ 1o6=6.5 percent. 
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Fic. 6. Creep curves of H-11 at 60°C. 1=extended length 
at time ¢, /,=unextended length, J)=initial extended 
length: k’=rate constant obtained from stress relaxation 
data at 60°C. logio vs. k’t, O lb =32.8 percent, @ /o=20.8 
percent, @ /)>=9.4 percent, @ /»=3.2 percent. 


elastic modulus, and k’ a relaxation rate. Equa- 
tion (4) would quantitatively explain the results 
shown in Figs. 3 and 4 and Eq. (2). However, 
since Eq. (1) indicates that the elastic stress- 
strain law is not linear (non-Hookean spring) it 
is clear that Eq. (4) could not be an exact 
representation of all the elastic viscous behavior 
of polysulfide rubbers, even were we to overlook 
the slight disparity in the fit of the relaxation 
curve to e~*’'. Results on creep, obtained by a 
previously described method,* bear this state- 
ment out very strikingly. Whereas the Maxwell 
Eq. (4) would predict a linear creep with time, 
the experimental data (Figs. 5 and 6) indicate a 
creep whose time dependence is even stronger 
than exponential. This cannot be accounted for 
in terms of changing cross section alone, but 
rather must be interpreted in terms of Eq. (1). 
A complete discussion of the theory of creep will 
appear in an accompanying article. 

Observations were also made on the mechani- 
cal behavior after removal of load during creep 
experiments. An initial rapid contraction fol- 
lowed by a small, slow creep recovery was 
noticed. If we had a system that could be strictly 
characterized by one reaction rate, as in Eq. (2), 
we could expect no creep recovery. Actually 
about 5 percent recovery is observed which isn’t 
surprising in view of the slight deviation from 
e—*’t evinced in Fig. 4. 

The theory previously derived which appar- 


3R. D. Andrews, R. B. Mesrobian, and A. V. Tobolsky, 
Trans. A.S.M.E. 67, 569 (1945). 
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64H 





Organic tetrasulfide 

Organic tetrasulfide 

Organic disulfide cross-linked 

Organic disulfide ? 

Organic disulfide tetrafunctional 
Organic disulfide 

Organic disulfide with 40% monosulfide 
Zinc oxide 

Pelletex . 
Micronex 

Stearic acid 

Altax 

DPG 

Tuads 

Lead peroxide 


Cure 





H3 


100 


10 
25 
50 


.10 
20 


50 X 287 


H5 H8 





30 X 287 


100 


10 — 
60 60 


30 X 287 


H9 


30 X 287 


H10 


100 


1.50 


30 X 287 
450 


Hil 


100 


10 
60 


50 
.30 
10 


50 X 298 
515 


H1i2 


100 


60 


S 


0 


ie 
1.50 
30 X 287 





Modulus 100 


300 — _ 
500 —_ _ 


450 
400 


Tensile 730 
Elongation 370 
Set a 


Hardness 74 70 


540 310 
— 1200 


340 
1475 1185 — 


1185 985 
300 280 
- ; 3 10 5 

73 68 74 76 64 


1475 
300 


1160 
190 420 
5 








ently satisfactorily explains the creep of Butyl 
and Hevea rubbers at elevated temperatures is 
strictly true only in the case that there is a 
scission reaction only. The satisfactory agree- 
ment for Butyl and Hevea occurs because in 
these rubbers the scission reaction predominates. 
This theory does not apply, nor can it be 
expected to apply, to the polysulfide rubbers 
where the rate of bond formation in the network 
structure is equal to the rate of bond breaking. 


RELAXATION STUDIES OF VARIOUS POLYSULFIDE 
RUBBERS OVER A RANGE OF TEMPERATURE 


In order to obtain some further insight into 
the nature of bonds whose rupture is responsible 
for the relaxation of stress in the polysulfide 
rubbers, a series of these rubbers prepared for 
us by the Thiokol Company were studied over a 
range of temperatures. The rubbers studied 
included two tetrasulfide polymers H-3 and H-S, 
two linear disulfides, a disulfide containing a 
fairly large percentage of monosulfide linkages, 
one disulfide containing two percent of trifunc- 
tional agent, and one disulfide containing one 
percent of tetrafunctional agent. Table I details 
these polymers and the compounding and curing 
recipe for the rubbers obtained therefrom. The 
relaxation curves for all these rubbers at 60°C 


is shown in Fig. 7. It is apparent that the 
relaxation rate varies by a factor of about one 
thousand from the fastest to the slowest of these 
compounds. The tetrasulfides relax the most 
rapidly; the inclusion of polyfunctional agents 
apparently slightly retards the rate of relaxation. 

The temperature dependence of the relaxation 
rates for these rubbers is shown in Figs. 8-14. 
In these graphs, stress divided by initial stress 
is plotted against logarithmic time for a series 
of different temperatures. Inasmuch as all the 
decay curves are approximated by the decay 
law e*’', a very simple method for obtaining 
the activation energies for stress decay was 
employed. As shown in Fig. 8 a horizontal line 
was drawn across the graphs at a value of 
f/fo=.368=1/e. The intersection of this hori- 
zontal line with the decay curves at the various 
temperatures gives the time for which k’ =1/t. 
Since in these curves time is already plotted on 
a logarithmic scale, a scale of reciprocal absolute 
temperature is plotted on the right-hand side of 
the figures along the ordinate scale. By locating 
the time for which ¢=1/k’ against the corre- 
sponding temperature on the reciprocal temper- 
ature scale, a plot of log 1/k’ vs. 1/T is obtained. 
The figures indicate that straight lines are ob- 
tained for all the rubbers studied and the slopes 
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Fic. 7. Continuous stress relaxation of various polysulfide 
rubbers at 60°C, 50 percent elongation. @ H-10, o H-2, 
e H-7, a H-11, © H-3, ® H-5. (See Table I.) 
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Fic. 8. Continuous stress relaxation of H-3 at various tem- 
peratures and 50 percent elongation; plot of log k’ vs. 1/T. 
O 35°C, @ 40°C, @ 50°C, © 60°C. 
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Fic. 9. Continuous stress relaxation of H-5 at various tem- 
peratures and 50 percent elongation; plot of log k’ vs. 1/T. 
oO 35°C, @ 40°C, @ 50°C, 60°OC. 


of the straight lines give a measure of the heat 
of activation. The heats of activation for the 
various rubbers are tabulated in Table II as are 
the entropies of activation. It is clear that 
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Fic. 10. Continuous stress relaxation of H-8 at vari- 
ous temperatures and 50 percent elongation; plot of 
log k’ vs. 1/T. © 60°C, @ 70°C, @ 80°C, @ 90°C, © 110°C. 
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Fic. 11. Continuous stress relaxation of H-9 at vari- 
ous temperatures and 50 percent elongation; plot of 
log k’ vs. 1/T. © 60°C, @ 70°C, & 80°C, S 90°C, © 100°C, 
© 110°C, @ 120°C. 
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Fic. 12. Continuous stress relaxation of H-10 at vari- 
ous temperatures and 50 percent elongation; plot of 
log k’ vs. 1/T. 0 60°C, @ 70°C, @ 80°C, © 90°C, © 100°C, 
© 110°C, @ 120°C. 


although the absolute values of the relaxation 
rates are very different (see Fig. 7), the heats of 
activation for all the polysulfide rubbers do not 
differ markedly. 








98 M. D. STERN AND A. V. 


357 
4032 


soy! 


6030 


70 55 














Fic. 13. Continuous stress relaxation of H-11 at vari- 
ous temperatures and 50 percent elongation; plot of 
log k’ vs. 1/T. 0 40°C, @ 50°C, © 60°C, @ 70°C, @ 80°C, 
© 90°C, @ 100°C. 


A pair of curves suffices to present the com- 
plete stress-time-temperature relations for a 
given rubber. For example, Fig. 15 shows how 
this can be done for H-11. A “universal stress- 
decay curve’’ is plotted against log k’t inasmuch 
as the decay curves are approximately super- 
posable by translation along the logarithmic time 
axis. A plot of log k’ against 1/T is shown on the 
same figure where the same logarithmic scale is 
used as abscissa, and the 1/T scale is plotted as 
the ordinate scale on the right-hand side of the 
figure. The values of k’ at any temperature 
for this graph can be readily obtained from 
Fig. 13 by taking the reciprocal of the time at 
the given temperature at which the value of 
f/fo=e". Figure 15 may now be used to give 
the complete stress-time behavior for relaxation 
at any temperature. The log-time scale at any 
temperature always has the spacing shown along 
the abscissa and may be located along the 
abscissa by using the straight-line locus of 
log k’ vs. 1/T. 


DISCUSSION 


The simple relaxation behavior (given by Eq. 
(2)) of any single polysulfide rubber studied here 
is strong evidence that a single type of bond 
in the structure is involved. Furthermore, the 
approximate equality of the energy of activation 
for stress decay in all the polysulfide rubbers is 
an indication that the relaxing bond is in each 
case closely similar, perhaps identical, even in 
the case of the disulfides and_tetrasulfides. 
Moreover, the facts indicate that the total num- 
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Fic. 14. Continuous stress relaxation of H-12 at vari- 
ous temperatures and 50 percent elongation; plot of 
log k’ vs. 1/T. 0 90°C, @ 100°C, @ 110°C, © 120°C. 


ber of structural elements does not change with 
time or temperature, so that the bonds involved 
reform at exactly the same rate that they break 
and the equilibrium between breaking and re- 
forming of bonds does not change with temper- 
ature. 

It is not possible to make a final decision at 
this point as to precisely which bond of the 
structure is responsible for relaxation. It is 
conceivable that a primary chemical bond along 
the chain is breaking and reforming, possibly 
by the simple expedient of bond interchange 
between two chains. On the other hand, it is 
also possible that the cross-linking bond between 
chains is in fact responsible for the thermoplastic 
behavior. It was hoped that the use of tetra- 
functional and trifunctional cross-linking agents 
would provide a means of tying the molecules 
into a three-dimensional structure. If these 
cross-linked polymers showed a limited relaxa- 
tion instead of a complete relaxation to zero 
stress, it would be possible to state unequivocally 
that the relaxation in the linear disulfides and 
tetrasulfides was due to a slippage of interchain 
bonds. On the other hand, there is no conclusive 
evidence that the tetra and trifunctional agents 











TABLE II. 
Polymer ss AHI (kcal) ASF (cal) 
H-3 25.41 +3.3 
H-5 20.24 — 10.0 
H-8 24.28 — 13.1 
H-9 25.70 —8.2 
H-10 25.08 —11.5 
H-11 22.92 —8.9 
H-12 24.47 —11.5 
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Fic. 15. ‘‘Universal curve” of H-11 and plot of log k’ vs. 1/T. 


used actually produce three-dimensional struc- 
tures, 

Considerations of the various possibilities put 
reasonable limits to the types of bonds responsible 
for polysulfide behavior. There are at least three 
bonds along the network chains, any one of 
which may be responsible for the observed phe- 
nomena. Whatever the bond, it must (1) have 
a low activation energy, (2) be formed as rapidly 
as it is broken throughout the temperature 
range, 35°-120°C, and (3) be insensitive to 
whether it is heated in vacuum or air. The three 
bonds along the polysulfide network chains are 
the C—C, C—S, and S—S bonds. Pauling‘ gives 
the C—C bond energy (per mole) as 58.6 kcal., 
the C—S as 54.5 kcal., and the S—S bond 
energy as 63.8 kcal. The C—S bond thus appears 
to be the weakest and the S—S bond the 
strongest of the three. 

The experiments on creep in vacuum and 
stress relaxation in oxygen-free nitrogen indicate 
that the relaxation behaviors of the hydrocarbon 
rubbers are due to oxidative attack, since the 
rate of scission is decreased a thousandfold in 
vacuum at 10-5-mm pressure.* Experiments on 
relaxation of polysulfide rubbers in vacuum at 
10-5-mm pressure indicate that oxygen is not 
responsible for their relaxation behavior (Fig. 
16). On the other hand, temperature does appear 
to have a primary effect on the relaxation rates 
of polysulfides both in vacuum and air as con- 
trasted to its secondary effect on hydrocarbon 
rubbers in vacuum.* Furthermore, the AH* of 

*L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1939.) 


* Unpublished experiments of W. C. Scheider and M. 
Magat and also unpublished work of present authors. 
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Fic. 16. Continuous stress relaxation of H-11 in vacuum 
and air at 70°C, 50 percent elongation. oO Air, @ 
Vacuum. 











Fic. 17. Continuous stress relaxation of H-9 and H-12 
at 120°C, 50 percent elongation. @ H-9 (disulfide), Oo H-12 
(monosulfide and disulfide). 


hydrocarbon rubbers of 30.4+2.0 kcal. as con- 
trasted to the AH* of 24+1.5 kcal. in poly- 
sulfides indicate that essentially one type of 
chemical reaction involving very similar chemical 
bonds is responsible for the relaxation behavior 
of all the hydrocarbon rubbers and that this 
reaction is not the same as the one responsible 
for polysulfide behavior. Thus it appears reason- 
able to assume that the hydrocarbon portion of 
the structure is not responsible for polysulfide 
behavior, although the possibility does exist that 
a C—C bond adjacent to a sulfur atom in a 
network chain may be sufficiently weakened by 
the presence of the sulfur atom as to satisfy the 
conditions stated above. The accessibility of 
d-orbitals in sulfur to partial double-bond forma- 
tion makes such a weakening of adjacent C—C 
bonds possible, and poses a problem for further 
research. 
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The relaxation experiments on H-12, a rubber 
consisting of both monosulfide and disulfide links 
along the network chains, indicate that there is 
little difference in its relaxation rate from that of 
a disulfide polymer (Fig. 17). The fact that there 
is little difference in the relaxation rates of the 
two rubbers indicates that the C—S bond is at 
least as strong, if not stronger, than the S—S 
bond. Furthermore, as pointed out in Gilman,‘ 
the formation of monosulfides (the C—S bond) 
from mercaptans occurs at high temperatures in 
the presence of catalysts such as metallic sulfides; 
the hydrogen attached to the a-carbon in mono- 
sulfides is more acidic than in corresponding 
ethers (pointing up the tendency of sulfur to 
expand its valence shell due to the accessibility 
of its d-orbitals, producing partial double bonds,® 
thus strengthening the C—S bond and weakening 
the C—H bond on the a-carbon); organic mono- 
sulfides do not react readily with elementary 
sulfur to form higher sulfides, whereas disulfides, 
trisulfides, and tetrasulfides readily add sulfur 
to form higher sulfides; conversely di and tri- 
sulfides may also be readily formed from higher 
sulfides by distillation. All these facts indicate 
that the S—S bond is more labile and less stable 
than the C—S bond, in contrast to their re- 
spective bond energies as given by Pauling. A 
further instance of this is provided by Powell 
and Eyring’ who have indicated an S—S bond 
energy of 27.5+5 kcal. as compared with the 
Pauling value of 63.8 kcal. 

In addition to bond scission along the primary 
valence chains of the polysulfide rubbers, it is 
also possible that relaxation is due to a slippage 
or rupture of bond between the polysulfide 
chains. It is not known at present whether the 
vulcanization process introduces primary cross 
links into the structure, and if so what the 
nature or energy of these cross bonds might be. 
It is conceivable that the cross linking in these 
rubbers is not due to the vulcanization process 


5H. Gilman, Organic Chemistry (John Wiley and Sons, 
Inc., New York), Vol. I, pp. 853-866. 

§ J. Rothstein, J. Chem. Soc., 1550, 1553, 1558, 1560 
(1940). 

7R. E. Powell and H. Eyring, J. Am. Chem. Soc. 65, 648 
(1943). 
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but instead occurs through formation of rela- 
tively strong bonds between sulfur atoms in 
adjacent chains. These bonds, which may pos- 
sibly involve the unshared electron in the sulfur 
valence shell, would be stronger than ordinary 
secondary bonds between chains, but weaker 
than most primary chemical bonds. The relaxa- 
tion process may consist in the breaking and 
reformation of these interchain bonds. 

The relaxation may also occur by a process of 
chain interchange to give two new chains which 
are chemically indistinguishable from the old 
ones and from each other. Such a process of 
chain interchange is illustrated below for the 
case of the disulfide linkage 


R—S—S—R R—S—S—R 

R—S—S—R R—S—S—R 

“ae he R—S . 
| KS 


R—S:S—R R—S S—R 


A mechanism of this or similar type would 
account for all of our observed phenomena. 

A final possibility exists, namely that some 
chemical agent introduced into the rubber in 
the emulsion or during compounding may be 
involved in breaking the polysulfide chains. 

It is apparent that a complete identification 
of the bond responsible for relaxation in poly- 
sulfide rubbers would be of great interest, and 
that this identification will accompany an en- 
hanced knowledge of the chemistry of mono and 
polysulfides. 
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<4 Near Ultraviolet Absorption of Pyridine Vapor 
in 
0S- H. SPONER 
lfur Duke University, Durham, North Carolina 
ary AND 
ker HILDEGARD STUCKLEN 
xa- Sweet Briar College, Sweet Briar, Virginia 
und (Received October 29, 1945) 
; of The absorption spectrum of pyridine at 3000-2500A has been studied in the first order of a 
sof 3-m grating spectrograph. The band system represents an electronic transition A:—B, (using 
symmetry C2, for pyridine). Several progressions of totally symmetric vibrations are observed. 
old The band 34769 cm is taken as 0,0 band. Besides carbon ring vibrations the occurrence of 
of hydrogen vibrations is established and discussed. 
the 
INTRODUCTION arc furnished the comparison spectrum. Eastman 
INCE the investigation of the near ultraviolet spectroscopic plates I 0 and II 0 were used. 
absorption spectrum of pyridine vapor by Exposure times varied from 20 minutes to 3 
a Henri and Angenot,! much theoretical progress hours. The measured wave-lengths of the bands 
has been made in the interpretation of polyatomic agree well with those of Henri and Angenot. 
-R spectra. As the analysis of the pyridine spectrum Wherever possible, the edges on the violet side 
is important as basis for an interpretation of the Were measured. The appearance of the bands will 
ild spectra of pyridine derivatives and as the latter be discussed later. 
are of physical, chemical, and biological interest, Table I contains our measurements, visual 
me it was felt that a reinvestigation of the pyridine intensity estimates and interpretations. We have 
in spectrum was well justified. A preliminary report observed much fewer bands than Henri and 
be of this analysis has been published recently.? Angenot. This might be owing to the circum- 
stance that their measurements cover a larger 
on EXPERIMENTAL variety of pressure and temperature changes than 
y- As in previous work the spectrum was photo- ours; for example we have not repeated a study 
nd graphed in the first order (dispersion 5.54A per of the temperature dependence which brings out 
n- mm) of a 3-m grating spectrograph (Eagle bands from higher vibrating levels in the ground 
id mounting) with a hydrogen discharge lamp of the state. A number of these bands, not measured by 
ordinary type as light source. The pyridine wasa_ ourselves, are included, as they fit into our 
commercial Merck product. It was purified* by progressions. 
distilling it twice under ordinary pressure and The intensity estimates refer to particular ex- 
re then a third time after letting it stand over perimental conditions which are stated in the 
ty potassium hydroxide for some time. The sub- column under “Intensity.”” Henri and Angenot 
a- stance was finally carefully introduced into used a more detailed classification of intensities 
s, quartz absorption tubes of 15- and 75-cm lengths based on a more precise study of the conditions 
ig by vacuum distillation. During exposures the under which the bands make their first distinct 
ar pyridine was kept in a side arm at temperatures appearance. Although we have made use, for our 
to varying from —60°C to +32°C, covering pres- interpretation, of their intensity data besides our 
4 sures from a few hundredths to 30 mm. Aniron own, we give in Table I only our own estimates. 
I. ” $¢, elton © Angenot, J. de Chim. Phys. 33, 641 Figure 1 represents the absorption spectrum of 
id (1936). See their previous references. pyridine vapor. 
‘We noe oto pain pote Fe hh ng og G. P LYSIS 
a of the Biochemistry Department of Duke Medical ey Cre an 
school and to Dr. M. B. Hal of she Physics Department The pyridine molecule belongs to symmetry 
Wilmington, Delaware) for the purification of the sample. group C2, like many substituted benzenes. But 
101 




















102 SPONER AND H. STUCKLEN 
TABLE I. Absorption bands of pyridine. 
Separation , Separation 
Wave- from 0—0 Intensity Wave- from 0—0 Intensity 
number band 4H-A This research Assignment number band 4H-A This research Assignment 
*32488 —2281 ew 0—1218—1063; 777. —- —992 vw 0—99 
0—1218— 2601+(139); 0- 398) — —601. 
0—2xX1141. 
$24 —945 w 0—945. 
* 503 —2266 ew 0—2X<601—1063; 
0—4X601+-(139). * 843 —926 vw 0+-(139) — 1063; 
—2X164—601. 
* 520 —2249 ew 0—1031—1218. 
* 860 -—909 w 
* 557 —2212 ew 0—995—1218, 
t=31°C 
* 606 -—2163 ew 0—945—1218. p=28.5 mm 
tube 15 cm 
* 677 —2092 ew 0—891—2601; 
0—1031—2601+(139). 878 —891 w 0-891: 
0+(139) —1031. 
* 704 —2065 ew 0—2X1031. 
901 —868 w 0—(395) —601+-(139)? 
* 799 -—1970 ew 0—1372—601; 
0—1031—945? * 957 —812 ew 0—2405? 
* 935 —1834 ew 0—891—945; 34 003 —766 w 0—601—164. 
0—1031+-(139) —945. = 
* 010 —759 ww 0—2X378? 
* 951 —1818 ew 0—1218—601. n 
057 —712 vw 0-712. 
* 970 —1799 ew 0—3X601. 
073 —696 vw 
* 986 —1783 ew 0—2X891. £5 
091 —678 vw 0—1218+-542. 
*33027 —1742 ew 0—601—1141. a 
111 —658 mw 0—2X6014+542. 
* 101 —1668 ew 0—1063—601; 
0+(139)—3X601. 120 —649 m vw 0—649. 
* 137 —1632 ew 0—1031—601. 168 —601 st m 0—601; 
0—1141+-542. 
* 163 —1606 ew 0—945—601—59. e 
191 —578 vw 0—1063—601+2 542; 
* 225 —1544 ew 0—945—601. 0+(139) —3 X601+-2 X542. 
* 278 —1491 ew 0—891—601; * 201 —568 ww 0+(139) —712. 
0—1491; 0—2X747? ss 
° 3 —558 ww 0—1141+581. 
* 310 —1459 ew 0—(395) —1063; " 
0—712—747? * 226 —543 ww 0—378—164?; 
0+ (672) —1218? 
* 397 —1372 ew 0—1372?; . 
0—995—374? 250 —519 w vw 0—1063+-542; 
0+(139) —2601+-542 
— 262 —507 
tube 75 em 282 —487 m w 0—1031-++542. 
* 410 —1359 ew 0—712—649?; 305 —464 mst mw 0+-(139) —601. 
0—2X945+542? 
322 —447 w 0—992+-542; 
* 496 —1273 w 0—1218—59. 0—1031+-581. 
* §12 —1257 ow 0—2X601—59. 364 —405 ew 0—405; 
' 0—945+-542. 
* §25 —1244 w —2X891+-542. 
374 —395 mw 0—(395); 
551 —1218 w 0—1218. 0-+(672) —1063? 
567 —1202 w 0—2X601; * 391 —378 ww 0—378. 
0—1141—59. 
408 —361 wd vu 0+-(672)—1031. 
* 587 —1182 ow 0—1031—151? 
427.  —345 mst mw 0—891+-542. 
* 628 —1141 w 0—1141; 
0- (395) — —747? 446 —323 mw w 0—1218+-899; 
0—2X164; 
* 654 —1115 ww 0—1063—59; 0+-(672)—992. 
0+-(139) —3 X601+-542. 
462 —307 vw ew 0—2X151? 
* 665 —1104 ww 0—(395) —712?; tee 
0—945—164? 501 —268 wd vw 0—2405+-542; 
0+-(672) —945? 
* 681 —1088 w 0—1031—59. 
516 —253 wd vw 0—1218+-968. 
* 692 —1077 0+(139) —1218. 
“ + 549 —220 vw 0—1218+4+-995; 
706 —1063 w 0—1063; 0—164—59; 
0+(139) — —2xX601. 04.672) 801. 
738 —1031 mw 0—1031. * 566 —203 w 0—1218+1017. 
* 758 —1011 ww 0—945—659. 592 —177 w 
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Separation Separation 
Wave- from 0—0 Intensity Wave- from 0—0 Intensity 
number band 4H-A This research Assignment number band 4H-A This research Assignment 
t=—16°C 350 581 vw 0+581; 
p=1 mm 0—1218+2 899. 
tube 15 cm 
t=-19°C 
605 164 w mw 0—405+ 241; p=.9 mm 
0—1031+-864; tube 15 em 
0—1063+899? 
393 624 w rw 0+(139)+2 542 —601; 
* 618 -151 w 0—378+-226. 0—(395)+-1017? 
* 636 —133 ww 0—1031+-899; 407 638 w vw 0—9924+-3 542. 
0—1218+2X542. 
441 672 mst m 0+-(672); 
657 —112 vwd wd 0—649+-542?; 0—2X164+-995. 
0—2X659? 493 724 vwd 
* 693 —76 ww 0—945+-864?; 521 752 wd 
0+-(672) —747? 
544 775 vwd 0—1218+2 995. 
711 —58 u m 0—601+-542. 
567 798 vwd 0+864—59? 
721 —48 w m 0—945+899; 
0—1063+-1017. * 607 838 wd 0+-(672)+542 —378?: 
0+995—151? 
735 —34 vw 0—1031+-995. 
633 864 mst m 0+864. 
769 0 mst st 0-0. 
668 899 ow ew 0+899. 
779 +10 eu vu 0—891+899? 
677 908 vw 0+-968—59; 
* 793 +24 w 0—1063+-2 542; 0+-(672)+241 
0+(672) —649. 
695 926 w 0+2X542—164. 
812 43 vw w 0—1141+1184? 
737 968 mw w 0+-968. 
822 53 rw w 0—1031+2<542; 
0—945+-995. 764 995 mst mst 0+-995. 
831 62 vu u 0—1031+-968. * Ti 1002 ow 0+(139) +864. 
841 72 vw u 0—945+1017; 786 1017 mud wd 0+1017; 
0+-(672) —601. 0—601+4-3 X542? 
849 80 u mw 0+(139)—59; 826 1057 vwd vod 
0—891+-968. 
855 1084 st st 0+2X542, 
874 105 ew vw 0—891+-995; 
0+(139)+995—1031. * 872 1103 ww 0+-(139)+968. 
908 139 m mst 0+(139). * 889 1120 ww 0—601+2 864. 
922 153 ew 0—(395)+-542. 908 1139 wd vod 0+-(139)+-995. 
950 181 vw w 0—(395)+581?; * 929 1160 w 0+2X581; 
0+-(672)+-542—1031. 0+-(139)+3 542 —601 
975 206 w mw 944 1175 w w 0+1175? 
989 220 ew vw 953 1184 mw w 0+1184? 
* 993 224 ww 0+-(672)+542—992; 984 1215 mw w 0+-(672) +542. 
0+226. 
989 1220 mw w 0+-(139) +2542. 
35 032 263 ew w 0—601+-864; 
0+(672)—405; *36 030 1261 ww 
0+-(672)+542—945? 
052 1283 wd vwd 
047 278 ew wd 0—2405+2 542. 
073 1304 wd vwd 
075 306 vw wd 0—712+-1017? 
091 1322 wd vwd 0+(139)+1184. 
100 331 mw mst 
137 1368 wd wd 0+-995+-542—164. 
117 348 w m 0—649+-995. 
160 1391 m m 0+2X995—601. 
147 378 w mw 0--542—164. 
172 1403 m mud 0+-542+4+-864. 
165 396 mw mst 0—601+-995; 
0+542—161? 220 1451 wd 
190 421 mw mst 0—2X59-+-542. * 230 1461 ww 0+3X542—164. 
* 202 433 vwd 0—649+-2542. 257 1488 wd 
224 455 ew rw 0+2X226. 271 1502 wd 
254 485 mw mst 0—601+2542; 281 1512 vw 0+542+968. 
0+2X241. 
307 1538 st 0+542+-995; 
311 542 st vst 0+-(672)+-864. 
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TABLE I.—Continued. 








Separation 


Separation 





Wave- from 0—0 Intensity Wave- from 0—0 Intensity 
number band 4H-A This research Assignment number band 4H-A This research Assignment 
341 1572 w 0+-(672)+899. 933 2164 wd 0+4X542? 
353 1584 vw 942 2173 vud 0+4X542? 
397 1628 mstd 0+3X542. 971 2202 vud 
438 1669 mud 0+-(672)+-995. 999 2230 vwd 
491 1722 wd 0+2X864. 37 026 2257 wd 
* 5614 1745 w 0+3581. ~ 069 2300 md 0+(139)+4 542. 
539 1770 mud 0+-(139)+3 542? 151 2382 md 0+3X995—601? 
* 569 1800 w 0+2899. 209 2440 md 
* 579 1808 ww 226 2457 wd 
622 1853 md 0+995+864. 268 2499 wd 
669 1900 vwd 309 2540 md 0+542+2X995. 
682 1913 mud 0+2X542+-995—164. 394 2625 md 0+3X542+-995. 
697 1928 mwd * 494 2725 w 0+2X864+995. 
704 1935 wd 0+2X968; * 583 2814 m 
0+2X995—59. 
599 2830 md 
739 1970 vw 
672 2903 md 0+3X968. 
760 1991 m 0+2X995. 
758 2989 wd 0+3X995. 
847 2078 std 0+2X542+995; 
0+ (672)+542+-864. 836 3067 mud 0+2X542+2995 
* 885 2116 mud wd 0+-(672)+542+899. *38 374 3605 mw 0+3542+2 995. 
904 2135 vud 0+-(139)+2995. 











Norzs: Intensities are roughly estimated from plates taken under the stated conditions. Those from the paper of Henri and Angenot (denoted by H-A) have been 
adjusted to our scale. The letters have the following meaning: ew—extremely weak, rw—very weak, w—weak, mu—medium to weak, m—medium, mst—medium to strong, 
vst—very strong, d—diffuse. Numbers in italics mean 1—1 transitions. Brackets indicate that the separation from the 0—0 band cannot satisfactorily be given in 
terms of vibrational frequencies. Asterisks indicate bands measured by Henri and Angenot only. Doubtful assignments are denoted by a question mark. 


unlike these, the nitrogen atom in pyridine re- 
places a whole CH group whereas in the mono- 
derivatives only the outer H atom is substituted. 
This should result in an intensification of the 
pyridine spectrum over that of benzene itself. On 
the other hand, the electronic structure of 
pyridine is much like that of benzene because the 
mz electrons still have Ds, symmetry. Hence 
pyridine should ‘“‘remember’”’ the spectroscopic 
characteristics of benzene, namely, a weak ab- 
sorption system forbidden by symmetry but 
made allowed through interaction of a vibration 
of suitable symmetry. Indeed, the maximum 
molar extinction coefficient of pyridine* (in 
hexane solution) is only about 10 times as great 
as that of benzene. This is of the same order of 
magnitude as the intensification in phenol and 
aniline absorption. It was therefore interesting to 
analyze the pyridine spectrum with a comparison 
of the benzene spectrum in mind. 


4C. W. F. oer and J. P. Wibaut, Rec. Trav. Chim. 
Pays-Bas 56, 573 (1937). 


If we assume that we can treat the pyridine 
molecule as a substituted benzene, then the 
levels between which the near ultraviolet transi- 
tion takes place will have the symmetry which 
the benzene levels acquire upon transformation 
from Ds, symmetry to C2, symmetry. Thus the 
pyridine spectrum should represent an 14,—'B, 
transition. The essential symmetry elements, 
which refer to electronic and vibrational levels as 
well, are given in Table IJ. As usual, C2” denotes 
a twofold axis in the molecular plane passing 
through the nitrogen atom and ga, denotes the 
molecular plane. The + and — signs indicate 
symmetrical and antisymmetrical behavior when 
symmetry operations with respect to these ele- 
ments are performed. The ground state of the 
pyridine molecule is totally symmetrical, i.e., A1. 
It can be seen from Table II that the transition 
A,—B, is allowed with a transition moment lying 
in the molecular plane perpendicular to the CN 
axis in the x direction. Since the transition 
moment M has the transformation properties of 





ae 
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a translation, denoted as TJ, it can also be seen 
_ from Table II that the product A1B,M,=A1B,Ai1 
contains B, and that the product A,BiM, 
—_ =A,B,B, contains A». Hence excitation of vibra- 
tions of type’ 6; and a2 may bring out the y and z 
components of the electric moment, otherwise 
forbidden. Column 5 contains the selection rules 
for Raman and infra-red spectra. Column 6 indi- 
cates how many of the 27 possible pyridine 
vibrations belong to the different symmetry 
classes. 

We may summarize briefly: we expect the 
pyridine spectrum to have the characteristics of 
an allowed transition, but keeping in mind that 
because of the resemblance of the electronic 
structures in pyridine and benzene the absorption 
strength will not be much greater than that of 
benzene. The 0,0 band should be among the 
strong bands. Bands resulting from excitation of — 3ssss— 
totally symmetric vibrations will make up the — 3s764— 
main part of the spectrum. In addition, a number ae 
of weak bands may occur which correspond to — 35633— 
the y and z components of the electric moment 
_ made allowed through interaction of vibrations 
- i ail ae: 35441— 
“4 From these remarks it would seem that the 
pyridine spectrum should have much resemblance 
to the corresponding spectra of simple substi- 


—0 +995 +2 X542 
—0 +2 X995 


—0+3 X542 
—0 +995 +542 


—0 +2 X542 


—0 +995 
0 +968 


—0 +864 


—0 +672 


35311— —0 +542 


Absorption spectrum of pyridine vapor. 





me tuted benzenes of symmetry C2, and that the 

the analysis could be carried through in close analogy 

_ to the analyses of these spectra. However, the 

ich appearance of the pyridine spectrum reveals at 3490s— | —o+139 

_— once characteristic differences from the spectra of : 

“" the monohalogenated benzenes or of phenol, for — 347¢0— | —0-0 5 
i example. In one respect it is even different from - 

ts, any of the substituted benzene spectra: while the ¥ 

a oone— | 0-405 +241 

bands of these spectra show pronounced degra- 

tes dation to the red, this is not the case with the 

es pyridine bands. We will take up the matter of the 34497 —0 —891 +542 

the structure of the individual bands in a later 

we section and first discuss the other differences  34305— a eal 

ra which concern the vibrational structure of the 

spectrum. : 34168— a 

the 

A. OCCURRENCE OF RING VIBRATIONS 

“a We have taken the strong band at 34769 em— 

~N as the 0—O band of the pyridine spectrum, 

ion If the symbol refers to the symmetry character of a 

of vibrational instead of an electronic level, small Greek 


letters are used. 
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TABLE II. Selection rules for pyridine (symmetry C2»). 








Selection rules 


Essential symmetry 
Infra- Number of 


Symmetry elements 





class Cw 2 Transl. Raman red vibrations 
A 1 “ }- + i. p a 10 
A» }- _ dp 1a 3 
By - + 7; dp a uw) 
By -- _ i dp a 5 








whereas Henri and Angenot had considered the 
band 35762 cm~ as the origin of the system. 
Their formulas representing five progressions 
contain five different ground state frequencies 
and one upper state frequency. The progressions 
include bands coming from high vibrational 
levels in the ground state. Some of the given 
intensities seem too large for the corresponding 
Boltzmann factors. For example, the band 34336 
which appears at 63°C at a pressure of 59 mm ina 
40-cm tube with an intensity 45 (given intensity 
45 means 100/45=2.0 mm is the pressure at 
which the band would appear in a 100-cm tube) 
is interpreted as 35762+4X542—6X600. This 
means that 6 quanta of 600 cm™ are excited in 
the ground state. The Boltzmann factor for the 
excitation of one quantum of 600 cm~ is about 
0.05, that of 6 quanta about 410-7. Another 
band 37023, which appeared in 100-cm layer at a 
pressure of 0.86 mm at 18.5°C with an intensity 
1000, was interpreted as 35762—1144—5 600 
+10 542 giving a Boltzmann factor of 210°. 
It was mainly the difficulties in reconciling some 
Boltzmann factors with the given interpretations 
and the occurrence of irregularities in the in- 
tensity distribution within the progressions which 
called for a new analysis. We have shifted the 
0—0 band to longer waves, and we have used for 
the analysis many Raman and infra-red vibra- 
tions in the lower and upper states. Although we 
also have small Boltzmann factors for a number 
of bands, yet they do not seem incompatible with 
the observed intensities. For example, the weak 
band 33551 which involves the excitation of 1218 
cm! in the normal state gives, at a temperature 
of 32°C, a Boltzmann factor of 210-*. Another 
band at 32520 .cm~!, measured by the French 
authors and resulting from the simultaneous 
excitation of the vibrations 1218 and 1031 in the 
ground state, gives a Boltzmann factor of 
2.0 10-* for the used temperature of 308°C. 
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Progressions involving the excitation of totally 
symmetrical carbon vibrations in the 1000-cm— 
frequency region appear with only medium in- 
tensity while they represent the strongest bands 
in benzene and benzene derivative spectra. We 
interpret the bands 35737 and 35764 cm as 
transitions from the vibrationless® ground state 
to singly excited symmetrical carbon ring vibra- 
tions in the upper electronic state. The values in 
the ground state are known from Raman 
spectra’—” and infra-red spectra!*® as 1028 and 
990 cm~. The values in the excited state are 995 
and 968 cm-', respectively. Both vibrations form 
progressions which, although strong at higher 
pressures, do not stand out as clearly as do the 
corresponding progressions in the substituted 
benzenes or in benzene itself. The 995 progression 
is the stronger of the two. 

As in the monohalogenated benzenes there 
should be in pyridine another symmetric carbon 
ring vibration not far from the discussed fre- 
quency range. Indeed, a Raman line has been 
reported at 1485 which is considered as the 
symmetric component resulting from the splitting 
of the 1485 e,-— in benzene. The other component 
is probably a little lower and has been assigned to 
the strong infra-red band’*!® at 1440 cm—!. We 
have not been able to assign definite fre- 
quencies to these vibrations in the excited state. 
While we have observed the vibrations 1028 and 
990 (ground state) in absorption our temperature 
range did not allow observation of bands re- 
sulting from high frequency vibrations in the 
lower state but Henri and Angenot have meas- 
ured a very weak band displaced by 1491 cm! to 
the red from the 0—0 band. 


® We call the unexcited ground state with its zero-point 
vibrations the vibrationless state for simplicity’s sake. 

7R. Ananthakrishnan, Proc. Ind. Acad. 3, 52 (1936). 

8 P. Krishnamurti, Ind. J. Phys. 6, 543 (1931). 

9K. W. F. Kohlrausch and A. Pongratz, Ber. d. d. 
chem. Ges. 67, 1465 (1934). 

10 C, S. Venkateswaran, Phil. Mag. 15, 263 (1933). 

11K, W. F. Kohlrausch, Zeits. f. physik. Chemie B30, 
305 (1935). 

12E. Herz, L. Kahovec, and K. W. F. Kohlrausch, 
Zeits. f. physik. Chemie B53, 124 (1943). 

133 W. W. Coblentz, Investigations of Infrared Spectra 
(Washington, D. C., Carnegie Institute, 1905), Part I, pp. 
96, 145, 272, 274, and insets. 

14]. Lecomte, J. de phys. et rad. 9, 512 (1938). 

% J. Turkevich and P. C. Stevenson, J. Chem. Phys. 11, 
328 (1943). 

16 Ch. H. Kline and J. Turkevich, J. Chem. Phys. 12, 300 
(1944). 
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The most characteristic feature of the pyridine 
spectrum, noticed and stressed by all observers, 
is a progression with separations of 542 cm=!. A 
frequency of 500-520 cm has been found also 
in the substituted benzene spectra but there are 
no pronounced progressions with separations of 
this frequency. In fact, the occurrence of bands 
with the single excitation of the 500 frequency in 
these spectra was explained by a forbidden 
transition with the moment in the y direction 
made allowed through the interaction of this 
particular vibration.!” It is the analog to the 
vibration, the excitation of which is necessary for 
the production of the benzene spectrum. It is 
true that bands, which very probably correspond 
to two- and perhaps threefold excitations of this 
vibration, appear in brom-, chlor-, and fluorben- 
zene!’ but they are weaker or much weaker than 
the band with the single excitation, and the bands 
caused by twofold excitations seem to have a 
structure different from that of the band resulting 
from the excitation of one quantum. If a similar 
explanation were correct for pyridine, the in- 
tensity of the 542 progression should be less than 
that of the progressions involving the symmetrical 
carbon vibrations in the 1000 cm region, and 
the intensity should fall off very rapidly. In 
reality, the 542 progression is the strongest in the 
spectrum and at least 4 members of it have been 
observed with certainty. This indicates that the 
progression does not represent a forbidden transi- 
tion. We therefore interpret the 542-cm™ fre- 
quency in the excited state as the symmetric 
component a; which results for symmetry C2 
from the splitting of the 606 e,+ benzene vibra- 
tion. Its value in the ground state is known from 
the Raman line 604 cm while we find 601 from 
the absorption spectrum (band 34168 cm-). 
Excitations of two and three quanta of this 
vibration have been observed in the ground state. 
The 542 progression occurs in all combinations 
with the 995 progression. Besides these combina- 
tions involving vibrations in the upper state 
only, there are also present bands in which one 
or both vibrations are excited in the ground state. 
Combinations of the 542 (or — 601) with the 968 
vibration are less frequent. 


(1941. Sponer and S. H. Wollman, J. Chem. Phys. 9, 816 
18S. H, Wollman, J. Chem. Phys. (to be published soon). 


The other component #; produced in pyridine 
by the splitting of the 606-benzene vibration 
should have a value close to 604. A fairly strong 
Raman line was observed at 652 cm~. However, 
both the 604 and 652 lines are depolarized!® ” 
and hence cannot represent the a; and 6; com- 
ponent of the 606 vibration. It has been sug- 
gested,!'6 therefore, to consider the 604 as 
unresolved containing both components a; and 
B:. This is not surprising because of the similar 
masses of benzene and pyridine. The question 
arises whether this is also true in the excited 
electronic state. If not, one should expect the 
band at 35311 cm~ (0,0+542) to be a doublet 
with one component much weaker than the other. 
Our plates do not indicate the existence of such a 
doublet nor can it be derived from Henri and 
Angenot’s measurements. We do not think that 
the very weak band at 35350 can be regarded 
as a component of the 35311 band. Hence the 
conclusion appears justified to consider the 542 as 
unresolved within the dispersion of our instru- 
ment. The band is broader than some other bands 
but that does not necessarily suggest superposi- 
tion of two close bands. 

From this discussion the strong intensity of the 
542 progression appears as a result of the fact 
that it represents a totally symmetrical carbon 
vibration superposed on an allowed transition (as 
explained before) and also the unsymmetrical 
component which corresponds to the vibration 
which makes the electronic transition in benzene 
allowed. 

The band 34711 cm, which is 58 cm= from 
the 0,0 band, is probably the 1—1 transition of 
the discussed vibration (601—542=59). A few 
other weak bands shifted by the same amount 
towards the red from strong bands may indicate 
the superposition of this 1—1 transition on sym- 
metric ring vibrations. The explanation has to be 
given with some caution, however, because the 
shifted bands are not always present at the 
expected places. 

The Raman frequency 652 cm~ was assigned 
to an out-of-plane vibration, a CH vibration by 
the Kohlrausch group and a ring vibration by 
Kline and Turkevich. The latter authors con- 
sidered an interpretation as $2 vibration as 
probable, corresponding to Bs, in benzene. Ac- 
cording to selection rulesa 82 vibration should not 
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occur singly excited in the discussed absorption 
system. If it would appear in violation of the 
selection rules for C2, symmetry, it should repre- 
sent an extremely weak band. This is not true for 
the band 34121 cm™ displaced by 649 cm to the 
red from the 0,0 band and probably involving the 
same vibration as the Raman line 652. Super- 
posing 542 on this band gives 34663 whereas a 
diffuse band has been measured at 34657 and 
might fit the explanation. Adding one more 
quantum leads to a weak band observed by 
Henri and Angenot. Another band at 35117 may 
be explained as 0,0 —649+-995 =0,0+ 346. There- 
fore, it does not seem very probable that the 
—649 difference occurring in the ultraviolet 
absorption spectrum can be explained as a Be 
vibration. Hence, the band 34121 has either a 
different origin or the 649 vibration must be 
assigned differently. We shall take up these 
questions later again in connection with the 
assignment of hydrogen vibrations. 

It is plausible to consider the band 35350 as 
containing the corresponding vibrational transi- 
tion of the 649 in the upper state. This gives a 
value of 581 cm~ for this vibration. The band is 
weak, so are the bands 35889 and 36341 resulting 
from superpositions of 542 and 995 on the 35350 
band. 

The e,* benzene ring vibration of 406 cm™ 
splits in pyridine into the components 374 and 
404 cm™. In benzene strong bands were found on 
the long wave-length side of all prominent bands 
at a distance of 160 cm which was explained as 
resulting from a 1—1 transition of this vibration. 
Assuming a similar behavior in pyridine we 
tentatively interpret the band 34605 (—164 from 
the 0,0 band) as 1—1 transition of the 404 com- 
ponent. As may be seen from Table I bands can 
be found displaced by 164 to the red from 35311 
(1542), 35855 (2542), 36397 (3542) and 
from combination bands of the 542 and 995 
progressions. Strangely enough they do not occur 
as companions of the 995 bands alone. This is one 
of the numerous little puzzles still present in the 
spectrum. The band 34003 (— 766 from 0,0) fits 
the interpretation —601—164= —765. 

Henriand Angenot report two very weak bands 
at 34391 and 34364 displaced to the red from the 
0,0 band by 378 and 405, respectively. It is quite 
possible that these separations correspond to the 


Raman lines 374 and 404 cm. The symmetry 
type of these vibrations is B, and ae, respec- 
tively.’ 1° The single excitation of a82vibration is 
forbidden in our absorption system while the 
single excitation of an a2 vibration produces a 
transition moment perpendicular to the mo- 
lecular plane. We have an indication of the 404 
band on some plates but were unable to locate 
the 374 band. Using 405 and the difference 164, a 
value of 241 cm™ is obtained for the vibration in 
the upper state. Twofold excitations of 405 and 
241 lead to allowed transitions. Indeed, the band 
33957 observed by the French authors may be 
interpreted as —2X405. 

It should be mentioned that the medium strong 
band 34908 fits the explanation 0,0+542—403 
=0,0+139 but it is hard to believe that the 
combination of an a; and an a: vibration produces 
a band of the observed intensity, particularly 
when the transition 0,0—403 is _ practically 
absent. A 542 progression can be built up on the 
139 band and subtracting 601 from it leads to 
another prominent band. We shall discuss this 
band later again. 

In the search for an occurrence of the 1—1 
transition of the 374 vibration we only found a 
weak band at 34618 observed by Henri and 
Angenot and displaced by —151 from the 0,0 
band. Using 378 from the ultraviolet absorption 
spectrum 226 cm obtains for this vibration in 
the upper state. Bands can be found in the places 
0,0—2X151, 0,0—2 378, 0,0+2 226, but sub- 
traction of 151 from the strongest bands leads to 
new bands in a very few instances only. We 
emphasize that the suggested scheme of the 151 
difference must be considered as very tentative. 


OCCURRENCE OF HYDROGEN VIBRATIONS 


The identification of hydrogen vibrations in the 
spectrum of pyridine is of particular interest 
since, with a few exceptions, their appearance 
could not be proved unambiguously in the ab- 
sorption spectra of benzene and its derivatives. 
An examination of pyridine plates taken at high 
pressures reveals at once the presence of many 
bands on the long wave-length side whose 
separations from the 0,0 band coincide with 
Raman lines and infra-red bands, and which 
cannot be explained as ring vibrations or combi- 
nations of such. Omitting separations resulting 








~~, _ a Eh hU6ehle 


try 
eC- 
n is 
the 
sa 
no- 
404 
cate 
4ia 
1 in 
and 
and 

be 


ong 
403 
the 
ices 
rly 
illy 
the 
; $0 
his 


—1 


ind 
0,0 
ion 
in 
ces 


3 to 
We 
151 
ve. 


the 
est 
ice 


ab- 


igh 
ny 
ose 
ith 
ich 
bi- 
ng 








ABSORPTION OF PYRIDINE VAPOR 109 


from ring vibrations, the following frequency 
values in cm~! denote such separations on the 
long wave-length side from the 0,0 band: 345, 
395, 464, 649, 868, 891, 945, 1063, 1218. Weaker 
bands appear at distances of 678, 712, 1141, and 
some very weak bands at 1372 and at higher 
values. The differences 712, 891, 945, 1063, 1218, 
1141, 1372 agree with frequencies observed in the 
Raman and infra-red spectrum. Assignments of 
these frequencies to modes of vibration have been 
offered by Kohlrausch and co-workers, and by 
Kline and Turkevich. Both groups consider them 
as hydrogen vibrations, but disagree in the indi- 
vidual interpretations. Their proposals for the 
hydrogen bending vibrations are collected in 
Table III. 

The bands —891, —945, —1063 and —1i218 
are all of similar intensity. The respective separa- 
tions appear also in other bands in conjunction 
with the symmetric ring vibrations 542, 968, 995 
in the upper and 601, 992, 1031 in the lower state. 
Furthermore, they combine with each other and 
with other hydrogen vibrations. The close 
agreement of the four separations with Raman 
and infra-red frequencies suggests that they 
represent fundamental vibrations. However, a 
detailed discussion of the experimental data will 
require some modification of this conclusion. 

When studying the spectrum of pyridine two 
characteristic pairs of bands are easily spotted: 
the strong pair 0,0 and 0+139 and the weaker 
pair 0—601 and 0+139—601 = —462 (see Fig. 1). 
A third still weaker pair is found shifted by 
another 601 cm~ towards longer waves: 0,0—2 
X601 and 0+139—2 601. The separations of 
these bands from the 0,0 band are —1202 and 
—1202+139= —1063. Henri and Angenot have 
observed even a fourth pair of extremely weak 
bands when moving another 601 cm~ to the red. 
As discussed in the preceding section the 601 is 
the frequency value of a symmetric carbon vibra- 
tion. If now the meaning of the separation +139 
would be known, all bands of the different pairs 
would be explained at once. The 139 band will be 
discussed later in some other connection but it 
may be stated at this place already that it is one 
of the few bands for which a satisfactory explana- 
tion cannot be offered at present. Assuming for 
the moment that the +139 band be explained, 
then the —1063 band would represent a transi- 


tion in which two vibrational levels of the 601 
vibration are excited in the ground state on top 
of the vibrations already involved in the +139 
band. If the —1063 band contains also the 
excitation of the hydrogen vibration 1061 cm =, 
the intensity ratio between the bands — 1202 and 
— 1063 should differ from that of the bands —601 
and —464, for example. Indeed, the first two 
bands are of similar intensity while the intensities 
in the second pair are about 2:1. Hence it may be 
concluded that the —1063 band represents a 
superposition of the two transitions. We have 
noticed that the 1063 frequency, besides com- 
bining with 601, occurs also in superposition with 
542 and 2542 but not with 995 and 968 nor 
with the corresponding vibrations in the ground 
state. The assignments of 1063 to a vibrational 
type (from Raman and infra-red data) are shown 
in Table III. Although the ultraviolet data do 
not permit a definite decision between the two 
propositions, an interpretation as a; or 6, vibra- 
tion is more probable from our point of view 
since excitation of an a, vibration is very likely 
connected with a very small transition proba- 
bility. It is not impossible that the diffuse band 
35786 corresponds to the excitation of the dis- 
cussed vibration in the upper state (1017 cm~) 
although Table I gives also an alternative expla- 
nation, but different transitions involving the 
same energy are not infrequent. 

The —891 separation can likewise be explained 
in a twofold way. First, it may be noticed that it 
is displaced by +139 from the —1031 band and 
hence may be a combination of the +139 
separation and the carbon ring vibration 1031. 
Although the intensity ratio of —1031: —891 as 
taken from Henri and Angenot’s paper is about 
2:1, we think that excitation of the CH vibration 
884 cm! (average value from Raman and infra- 
red data) is partly responsible for the band. This 
may be concluded from the fact that more 
superpositions of ring and hydrogen vibrations 
can be found on the —891 band than on any 
other band on the long wave-length side from the 
0,0 band. Of the two suggestions (Table III), to 
assign the frequency 884 cm= to an a, or ae 
vibration we prefer the first, but regard an 
interpretation as $; vibration also as possible. 
The band at 35633 with a separation of 864 from 
the 0,0 band is tentatively considered as repre- 
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TABLE III. Hydrogen bending vibrations in pyridine. 








In the molecular plane Perpendicular to the molecular plane 














Type in Den (3) a2 (15) Bou (9) e,* (18) eu (11) au (5) Bx (10) eg (17) eu* 
Type in C2, Bi Bi a1 Bi a1 Bi Be 2 a2 Be a2 Be 
Kohlrausch et al. 1370 1290 1218 1139 1069 wm m 747 890 652 710 
Kline and Turkevich 1210 1139 1139 om 884 1037 | 669 941 747 710 1063 m™ 
| 
Numbers in brackets in the first row refer to Wilson's numerical notation. m =missing denotes hydrogen vibrations which disappear in pyridine 


as compared to benzene. 


senting the discussed vibration in the excited g vibrations. But it is peculiar that a 542 pro- 


electronic state. The 542 and 995 appear super- 
posed on it and 2864 has also been observed. 

No double interpretation seems possible for the 
—1218 and —945 bands. A medium strong 
Raman line with a depolarization factor between 
0.6 and 0.75 was found at 1218 cm~. The fre- 
quency occurs also in the infra-red. From this an 
assignment was suggested to a bending vibration 
in the molecular plane, namely a (1 vibration 
(Kline and Turkevich) or an a; vibration 
(Kohlrausch and others). Both explanations are 
compatible with the relatively strong intensity of 
the —1218 ultraviolet absorption band. No 
preference can be derived from ultraviolet ab- 
sorption data. Adding 542 to —1218 obtains 
—676 which explains another of the mentioned 
separations. It is not impossible that part of the 
intensity of the —676 band is caused by the 
superposition of a fundamental vibration sug- 
gested by Turkevich and Stevenson as 669 cm7 
from infra-red data. It is not certain whether the 
band 35944 or the band 35953 can be interpreted 
as transition to the single excitation of the 1218 
in the upper state. If so, the corresponding fre- 
quency would be 1175 or 1184 cm. Superposi- 
tions of 542 and 995 on 1175 or 1184 fall close to 
or on top of other more likely transitions. 

The fundamental frequency 941 cm™ is known 
from a very weak Raman line and from an infra- 
red band observed with medium intensity in the 
vapor spectrum. Kline and Turkevich assign it to 
a B2 vibration resulting from the benzene vibra- 
tion B2,. It is very difficult to reconcile the 
intensity of the ultraviolet band —945 with this 
interpretation. In fact, the intensity of this band 
would suggest rather an a or §; vibration. A 
number of bands may be explained by super- 
posing on —945 the ring vibrations in the 
1000-cm™ region and the —601 difference; there 
are also combinations of 945 with other hydrogen 


gression cannot be built up on it. —945+542 
gives a separation of —403 from the 0,0 band. 
There is an extremely weak band at this location 
which we have interpreted before as resulting 
from the excitation of an out-of-plane carbon 
ring vibration. The band 34374, which is of an 
intensity comparable to that of the —945 band, 
is displaced by only —395 from the 0,0 band. No 
satisfactory interpretation of this difference can 
be offered in terms of other frequencies. However, 
superposing 2542 on —945 leads to the strong 
band +139 with an intensity too strong for the 
Boltzmann factor of a 945-cm~ vibration. Actu- 
ally, this band appears almost simultaneously 
with the 0,0 band and hence must have a 
different origin. To explain it as resulting from a 
1—1 transition is not a likely possibility. The 
vibration would have, contrary to the rule, a 
larger frequency in the excited than in the 
normal state, and because of the intensity of the 
band, it would have to be a low vibration. Not 
only would 139 cm-! then be a relatively large 
change but it cannot obviously be connected 
with any of the few low frequency vibrations. 
Adding 542 on 139 gives again a strong displaced 
band (35441 instead of 35450, in Table I the 
separation of 35441 from the 0,0 band is used as 
(672)). The succeeding member is regular and 
weak. Attempts to explain the first displacement 
as caused by resonance with the 405 vibration 
lead to difficulties because of group theoretical 
selection rules. If 945 represents a 8; and 405 an 
a2 vibration resonance should not take place. If 
they were both of type az resonating interaction 
would be possible between 405 and —945+542. 
But in this case 945 should not occur in the infra- 
red whereas Turkevich and Stevenson report a 
medium strong band in the vapor. Also the 
intensity of the ultraviolet band —945 is high for 
an a» vibration. We feel that we should point out 
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TABLE IV. Pyridine vibrations occurring in ultraviolet absorption. 








Carbon ring vibrations 


Hydrogen vibrations 





Lower state 1485 1031 992 [649] 601 405 1372? 1218 1139 1063 945 891 712 676? [649]? 


Upper state 995 968 [581]? 542 241? 


Type ay ay ay ? ay ae 





1175? 
or 1017? 899? 864? [581]? 
1184? 
ay ay ay 
By or or Bi or ? ? ? 


Bi Bi Bi 








[ ] means uncertainty as to ring or hydrogen vibration. 


these difficulties but we must admit that we 
cannot offer a convincing solution of them at 
present.!® The frequency of the 945 vibration in 
the upper state is very tentatively suggested as 
899 cm calculated from the weak band 35668. 
Addition of 995 leads to a band of similar in- 
tensity, 2899 gives a band of double intensity. 
Superpositions of 542 could not be found. 

There remain the differences 712, 1141, 1372 
which agree also with Raman and infra-red fre- 
quencies. The band 34057 with a separation of 
712 to the red from the 0,0 band is very weak. 
This separation may correspond to the medium 
strong infra-red band 710 cm™ which has been 
explained as 82 from ¢€,~ or e,*. Kline and Tur- 
kevich assume the other component a, from e,~ 
at 747 cm whereas Kohlrausch and co-workers 
prefer an assignment to 82 coming from Bey. (See 
Table III.) A fundamental frequency of 747 cm= 
is not present in the ultraviolet absorption 
spectrum. We have not been able to settle with- 
out ambiguity the question of the nature of the 
712:and 747 frequencies. 

The 1141 difference represents a weak band 
and corresponds probably to the Raman and 
infra-red band 1139 which was interpreted as 
symmetric component originating mainly from 
e,*. The 1141 difference occurs combined with 
—601 and 542. 

The infra-red band 1374 was tentatively as- 
signed to a 8; vibration by the Kohlrausch group 

1? Tt would be interesting to investigate the absorption 
of a-picoline (2-methyl-pyridine CsNHsCH3) which is a 
likely impurity. But although our purification method was 
not extreme, such impurity bands should occur with very 
small intensity only. Henri and Angenot mention in their 
paper that they have photographed the a-picoline spectrum 
and they would have noticed coincidences in both spectra. 
Unfortunately, their results on a@-picoline have evidently 


not been published. Earlier investigation by E. Herrmann, 
Zeits. f. wiss. Phot. 18, 253 (1919). 


but considered as combination or impurity by 
Kline and Turkevich. In ultraviolet absorption it 
shows up in an extremely weak band (33397 cm~). 
Both possibilities are mentioned in Table I. . 

Before concluding the section on hydrogen 
vibrations we should take up once more the 
question of the origin of the 649 difference. The 
corresponding Raman line was assigned by the 
German authors to the optically inactive com- 
ponent, i.e., a2, of the e,+ CH vibration in 
benzene. They remark that the intensity of the 
line is unusually high for this interpretation. The 
fact that since then Turkevich and Stevenson! 
located a medium strong band at 650 cm™ in 
the infra-red vapor spectrum and the fair in- 
tensity of the corresponding ultraviolet band 
render the interpretation as a2 doubtful. Because 
of these and other difficulties mentioned before, 
no assignment is suggested for this frequency 
difference at present. 

Reviewing our results in comparison with those 
from Raman and infra-red spectra we found 
agreement on the interpretation of most ring 
vibrations but deviations in the interpretation of 
hydrogen vibrations. We made some effort to find 
an unambiguous set of hydrogen bending vibra- 
tions which ewill satisfy the results of all three 
spectra. This could not be achieved. It seemed 
preferable, therefore, to tabulate separately 
(Table IV) the vibrations which occur in ultra- 
violet absorption and to indicate their vibrational 
type as suggested from the ultraviolet spectrum 
without correlating the types with those of the 
benzene vibrations. 


APPEARANCE OF INDIVIDUAL BANDS 


While benzene and the known substituted 
benzene spectra show pronounced degradation of 
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the bands to the red, most pyridine bands have a 
line-like appearance. Although they too have 
edges on the violet side with a very narrow shading 
toward the red, a number of them show also 
fuzziness on the violet side of the edge. Under 
suitable pressure conditions fine structure can be 
distinguished in this fuzziness in some bands, as in 
the 0,0 band at 34769 cm—, the first two members 
of the 542 progression and a few others. The 
extension of the fuzzy region is larger on the 
violet side than on the red side, but both never 
become as pronounced as the shading (to the 
red) in the substituted benzene spectra. Super- 
ficially the bands look like more or less broad 
lines with a sharper edge on the violet side and at 
higher pressures they are accompanied by narrow 
continuous regions (which broaden with pressure) 
on either side. 

Such a structure might be understood if the 
moments of inertia are very nearly the same in 
both the lower and upper states. Considering in 
first approximation the pyridine molecule as a 
symmetrical top like the benzene, its rotational 
energy can be represented by 


i? 
E, leila gh 
where A is the principal moment of inertia about 
the z axis (perpendicular to the molecular plane). 
The main branches P, Q, and Rare determined by 
AJ=-—1, 0, +1. If A’ and A” (moments of in- 
ertia in upper and lower state) are approximately 
equal, the lines of the Q branch fall nearly in the 
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same place (vo). The Q branch is then intense and 
narrow while the heads of the P and R branches 
(in our case R branch) lie, for A’ =A’, at some 
distance from the zero line or the Q branch. This 
distance from the origin may be so great that for 
the corresponding J value the intensity of the 
lines may have fallen off so much to be too weak 
for observation. This would account satisfactorily 
for the fact that a number of bands look as if 
degraded to both sides with an intense head in 
the middle. That in some bands an indication of 
structure can be seen on the violet side’ may be 
explained by assuming non-overlapping of the R 
branch on this side by the Q branch. This may 
account partly for the different extension of the 
fuzzy regions. Superimposed on this coarse struc- 
ture is the K structure. Much narrower p and r 
branches according to AK=+1 or narrow gq 
branches according to AK = 0 occur superimposed 
on the main P, Q, and R branches. The bands 
AK = +1 correspond to changes of the transition 
moment in the x, y directions and the bands 
AK =0 to a change in the z direction. It is of 
course hopeless to observe the K structure, even 
with very large dispersion. In spite of this we 
know that most bands belong to the type 
AK=-+1 because the vibrational analysis has 
revealed that in all strong bands the transition 
moment lies in the molecular plane. 

One of us (H. Sp.) acknowledges with pleasure 
the opportunity of discussing some questions 
with Dr. A. L. Sklar. We also wish to acknowl- 
edge a grant-in-aid from the Rockefeller Founda- 
tion for carrying out this research. 
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HIS section will accept reports of new work provided 

these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 





Latin American Collaboration 


HARLOW SHAPLEY 


Chairman of the Committee on Inter-American Scientific Publication, 
Harvard College Observatory, Cambridge, Massachusetts 


November 8, 1945 


CIENTISTS in the United States are cooperating 

with their colleagues in Latin American countries by 
contributing files of scientific journals to a number of new 
and acitve institutions in these countries where lack of 
such facilities is greatly impeding research. This work is 
being carried on by the Committee on Inter-American 
Scientific Publication, headed by Harlow Shapley, and 
already the libraries of two such institutions have been 
assisted in this way. 

The committee is now collecting journals for a number 
of other institutions in Latin America and is particularly 
eager to receive contributions of unused files of chemical 
periodicals for this purpose. We hope that any persons 
who wish to assist in this very genuine and timely method 
of promoting hemispherical solidarity and the progress of 
science, and who have even short runs of unused chemical 
periodicals, will communicate with the Comité Inter- 
americano de Publicacién Cientffica, Harvard College 
Observatory, Cambridge 38, Massachusetts. 





Creep under Constant Load in Tension 


PAUL J. BLATZ AND ARTHUR V. TOBOLSKY 
Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 
January 11, 1946 

NE of the most commonly observed and measured of 

the phenomena pertaining to elastic-viscous behavior 

is creep under constant load. The simplest postulated 

elastic-viscous body, namely one which obeys Maxwell’s 

equation, should exhibit a creep that is linear with time. 

It has been shown, however, that elastic materials that 

obey the kinetic theory of elasticity and are characterized 

by a single relaxation time, although they obey Maxwell’s 

equation in shear, are governed by the following relation 

between stress, strain, and time for simple tensile 
deformations, 
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where f is the stress on attained section at time ¢, s is the 
number of network chains of the elastic polymer per cc, 
k is Boltzmann’s constant, T is the absoluté temperature, 
k’ is the relaxation rate, /(t) is the length at time #¢, l(t’) 
the length at time ?#’, J, the unstretched length, and ¢ a 
constant.! 

If, as is often true, the volume of the elastic material is 
unchanged during stretching, the following relation ob- 
tains for an experiment for which a constant load is sus- 
pended on the material (creep under dead load). 


VO} 


where Jy is the initial length obtained immediately after 
the load is placed on the material. In other words, Jo is 
the length attained by the sample under load due to the 
purely elastic response of the material before creep due to 
internal relaxation processes is observable. 

If Eq. (2) is substituted into Eq. (1), the resulting 
equation may be transformed into a differential equation. 
For simplicity the following changes in notation were 
made: 





rA=1(t)/l., Ao=lo/lu, r=k't. 


The resulting differential equation is: 
nr” Qn’ r\2 x’ 
in ae 3 
YU 1+ x 2(*) By (3) 


B =3/(Ao— 1/Xo?) 
and 0’ and 0” are first and second derivatives, respectively, 
with respect to ¢. 

The solution of this equation was obtained by Laplace 
transformations giving: 


na. we 21 — aftu_ bu 
exp ()=P(5--1)+0 5. +1) exp | a(;: *)I. 


where 
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Bro (4) 
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Graphs of this function in the form of In //l, plotte@ The dissociation constants for ¢-butylamine-trimethyl- 


against k’t are given (Fig. 1) for several values of Jo/l.. Itis 
particularly to be noted that the creep rate shows a very 
marked dependence on the load (as measured by the initial 
extension) even though the molecular relaxation is as- 
sumed independent of the stress. It should be remem- 
bered that these creep curves apply only to substances 
like the polysulfide rubbers? exhibiting a single relaxation 
time. Creep data obtained by Mr. Mochulsky on several 
polysulfide rubbers are in good agreement with the theory 
and will be published soon. 

We wish to acknowledge the important mathematical 
advice we received from Professor V. Bargmann of the 
Institute for Advanced Study. 


1M. S. Green and A. V. Tobolsky, J. Chem. Phys. 14, 80 (1946). 
2M. D. Stern and A. V. Tobolsky, J. Chem. Phys. 14, 93 (1946). 





Steric Strain in Compounds with Tertiary 
Butyl Groups; Steric Strain in 
Polyisobutylene 


HERBERT C. BROWN AND GERALDINE K. BARBARAS 


Department of Chemistry, Wayne University, Detroit, Michigan and — 
George Herbert Jones Laboratory, University of Chicago, Chicago, Illinois 


January 11, 1946 
Dye POUNDS containing two tertiary butyl groups at- 
tached to oxygen, nitrogen, or carbon, such as 
di-t-butylether, di-t-butylamine, and di-t-butylmethane, 
are characterized by difficulty of formation and relatively 
low stability. Examination of molecular models points to 
a considerable degree of steric interference (indicated by 
the arrow in I, II, and III) between the two tertiary butyl 

groups joined to the common atom. 


O NH 
fo, a 
Jas, LaF, 
c—C-—C CC-C c—C—C CCC 
| ] | | ] | 
. C C C 
I II 


This interference must result in a considerable increase 
in the bond angle @ (I-III) from its preferred value and 
thereby lead to marked strain. In order to obtain an esti- 
mate of the strain, we have measured, by the experimental 
method previously reported,! the heat of dissociation of 
t-butylamine-trimethylboron (IV), a molecule whose struc- 
tural skeleton and molecular dimensions are similar to 
those of the molecules under consideration. 


CH: NH: 
™, ~~ 
J a La, 
C—C—C C—C—C C—C—C C—B—C 
- | | | ] | 
’ Cc & Cc 
III IV 


It is not anticipated that the minor differences in the 


atomic radii of boron, carbon, nitrogen, and oxygen would 
bring about any major differences in the magnitude of 
the strain. 


boron between 50° and 90° are represented by the equation 
log K = —2850/T+8.618. 

The heat of dissociation of the compound is 13.0+0.2 

kcal. per mole. 

Previous results have shown that where steric hindrance 
is not a major factor, the heat of dissociation of the 
addition compounds of trimethylboron with primary 
amines varies between the limits 17.6 and 18.4 kcal. per 
mole, with an average value of 18.0 kcal. per mole.? The 
steric strain in ¢-butylamine-trimethylboron may therefore 
be estimated to be in the neighborhood of 5 kcal. per mole. 

This estimate receives interesting confirmation from the 
data reported by Evans and Polanyi? for the heat of 
polymerization of isobutylene to polyisobutylene. The 
experimental value, 12.8 kcal. per mole, is some 9 kcal. 
below the value of 21.4 kcal. per mole calculated for a 
strain-free molecule. Examination of a scalar molecule for 
polyisobutylene (represented schematically by V) reveals 
that the strain in this molecule is similar both in kind and 
in degree to that postulated for the simpler substances, 
I-IV, with one noteworthy difference. 


CH, CH: CH, CH:e CH: 
Pa. a al “a . 
c—C—C C—C-—C C—C—C C—C-—C 
t t t 
Vv 

In the simpler molecules containing two tertiary butyl 
groups united to a common atom, each tertiary buty!] 
group encounters steric interference on but one side (I-IV). 
In the polyisobutylene molecule, however, each unit cor- 
responding to these tertiary groups encounters steric in- 
terference on two sides (the end groups are neglected in 
this discussion). It follows that the over-all strain per 
tertiary butyl unit in polyisobutylene should be twice 
that per tertiary butyl group in the simpler molecules, 
I-IV. In view of this consideration, the value of 9 kcal. 
per mole arrived at by Evans and Polanyi for the steric 
strain in polyisobutylene is in excellent agreement with 
the value of 5 kcal. per mole that our studies indicate to 
be present in ¢-butylamine-trimethylboron and related 
substances. 

1 Brown, Taylor, and Gerstein, J. Am. Chem. Soc. 66, 431 (1944). 


2? Unpublished work of Brown and Taylor. 
3 Evans and Polanyi, Nature 152, 738 (1943). 





The Absorption Spectrum of Metal- 
Ammonia Solutions 


RICHARD A. OGG, JR. 
Department of Chemistry, Stanford University, California 
January 14, 1946 


EW experimental studies of the visible absorption 

spectrum of dilute solutions of alkali and alkaline 
earth metals in liquid ammonia reveal complications ap- 
parently not recognized by previous investigators. On 
cooling fairly dilute (some 10~* molar) solutions of sodium 
and magnesium from room temperature to about — 80°C, 
there is a very noticeable darkening, together with a 
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change in the quality of transmitted light.. These phe- 
nomena are the more pronounced, the greater the con- 
centration of the solution. The color change is rapid and 
completely reversible with respect to temperature change. 

Quantitative spectrophotometric measurements confirm 
the qualitative visual observations as to change in color 
distribution of transmitted light. The absorption co- 
efficient of a given solution for blue light is but slightly 
increased by the above temperature change, while that for 
red light is increased considerably—by a factor of two or 
more for sufficiently concentrated solutions. The change 
for light of intermediate wave-lengths falls smoothly be- 
tween these extremes. 

Previous studies’ ? have shown the absorption spectrum 
of very dilute metal-ammonia solutions to be continuous, 
with a pronounced maximum in the infra-red. Were this 
absorption band to undergo the normal broadening with 
temperature increase, the visible absorption spectrum 
would be affected in a fashion the reverse of that described 
above. The most probable interpretation of the present 
studies postulates the presence in the solutions of two 
different colored components, the equilibrium between 
these being dependent upon concentration and tempera- 
ture. The component favored by high concentrations and 
low temperatures (which will be termed B) would appear 
to display also a continuous absorption spectrum, but 
with a band maximum lying nearer the visible region than 
that of the other component (which will be termed A). 
This would account for the observed differential increase 
of total absorption toward longer wave-lengths on cooling 
a relatively concentrated solution. The experimental re- 
sults would indicate that the conversion of B into A is an 
endothermic process and must involve an increase in the 
number of particles. 

Comparison with solid solutions of alkali metals in 
alkali halide crystals* affords a clue as to the difference 
between the A and B components. In these systems the 
F centers and F’ centers (consisting, respectively, of single 
electrons and electron pairs trapped at vacant negative 
ion lattice sites) show different continuous absorption 
bands, the maximum of the F band lying at shorter wave- 
length. Absorption in both bands results in ejection of a 
conduction electron, the corresponding work function 
being greater for the F centers. Correspondingly, the F’ 
centers are thermally unstable with respect to dissociation 
into F centers. In the present studies, it is suggested that 
the A component consists of single electrons trapped in 
cavities in the liquid dielectric medium‘ while the B com- 
ponent consists of electron pairs similarly trapped. How- 
ever, the order of thermal stability of A and B components 
would appear to be reversed as compared with that of F 
and F’ centers. Correspondingly, the positions on the wave- 
length scale of the corresponding band maxima are in all 
probability also reversed. 

Theoretical treatment! of the electron pair trapped in a 
liquid medium of high dielectric constant predicts it to 
be stable with respect to dissociation into single electrons 
trapped in separate cavities, and that its ground state 
should be diamagnetic. The trapped single electron should 
have a paramagnetic susceptibility corresponding to elec- 
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tron spin. Experimental studies’ * of the concentration 
and temperature dependence of the magnetic suscepti- 


’ bility of dilute metal-ammonia solutions are at least 


qualitatively explained by the above proposed equilibrium 
between trapped pairs and trapped single electrons. How- 
ever, in attempting to calculate equilibrium constants 
from the magnetic data, it should be noted that the cavi- 
ties occupied by both are very large‘ (of the order of 
7X 10-8 cm in radius) and that consequently the diamag- 
netic susceptibilities are far from negligible. The absorp- 
tion spectrum studies offer hope of greater accuracy. 

1G. E. Gibson and W. L. Argo, J. Am. Chem. Sdc. 40, 1327 (1918). 

2 E. Vogt, Zeits f. Electochemie 45, 597 (1939). 

3 For literature references, see N. F. Mott and R. W. Gurney, Elec- 
tronic Processes in Ionic Crystals (Oxford Press, 1940). 

4 Richard A. Ogg, Jr., J. Am. Chem. Soc (January, 1946), in press. 


5 E. Huster, Ann. d. Physik 33, 477 (1938). 
8S. Freed and N. Sugarman, J. Chem. Phys. 11, 354 (1943). 





Copolymerization of Systems Containing 
Three Components 
TURNER ALFREY 
Polytechnic Institute of Brooklyn, Brooklyn, New York 
AND 
GEORGE GOLDFINGER . 
University of Buffalo, Buffalo, New York 
January 16, 1946 


N a previous letter! on this subject, we presented equa- 

tions for the composition of the initial copolymer 

formed froma given mixture of 3 monomers, (A, B, and C).? 
_ [8% AB+y"8°A*+78°AC] [A +e°B+0¢C] 


a 
b~ [oyAB+ ya B+PaBC] [B+RA+BC) 


a_[BPAB+y-6A?+y"8AC] [AteB+eC] 
c [eatp*AC+BaeBC+BaC?] [C+yA+yYB] ~~ 








Here a, 6, and ¢ are molar concentrations of the three com- 
ponents in the initial polymer; A, B, and C are the molar 
concentrations in the monomer mixture; a, a*, —, etc. 
are ratios of propagation rate constants, as follows: 


o? = Rap/Raa, B* = Real kes, ¥* = Real kee; 
of = Rac/Raa, B° = Roc/ Rw, oy =RalRee. 


It was stated that “equations 1 and 2 indicate the possi- 
bility of predicting the composition of a multipolymer 
from information gathered from a separate study of each 
pair of monomers involved.” This is, in general, true. 
Dr. G. L. Barrett? has pointed out, however, that there 
is an important class of monomers for which such a cal- 
culation is impossible. Members of this class (e.g., maleic 
anhydride, maleic esters) possess vanishingly small rate 
constants for self-propagation.‘ If monomer C in a 3- 
component polymerization is of this type, then k,. is equal 
to zero; y* and 7’ are infinite. The structure of the tri- 
polymer depends upon all 9 of the propagation constants. 
In the usual case, where no constant vanishes, knowledge 
of the six ratios used in Eqs. (1) and (2) means implicit 
knowledge of all other necessary ratios. Thus (ea/k.) is 
fixed if the ratios (Rea/kec) and (ka/kee) are known. If Ree 
vanishes, however, experiments on the separate 2-com- 
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ponent systems yields only the fact that both (k./Rec) 
and (kca/Rec) are essentially infinite. Knowledge of (Rea/k.) 


must be obtained somehow, even though this ratio does - 


not appear explicitly in Eqs. (1) and (2). Letting Ree ap- 
proach zero, Eqs. (1) and (2) take the form: 


a_[@°AB+R6°A2+R6°AC] [A+e°B+a°C] 
6 [Re®&AB+o°B?+0°BC] [B+R°A+e°C]’ 


a_ [B°AB+R6°A?+RBAC] [A+a°B+a°C] 
c LapAC+arpBC+or6°C?] [RA+B] ’ 


R = Rial Res. 


(3) 








(4) 


where 


If a’, a’, B*, B°, y*, and y® are known from the separate 
2-component studies, a single tripolymerization experi- 
ment will yield more than enough information to deter- 
mine R. (R can be determined by Eq. (3), which allows 
Eq. (4) to be used as a check.) 

There are also cases where ‘wo monomers (say B and C) 
have vanishing rates of self-propagation. If ki. and kp as 
well as ky and k,-- vanish, then all necessary ratios can be 
determined from separate copolymerization reactions using 
the monomers two at a time, and the equations relating 
copolymer composition to monomer composition take on 


the simple form: 


a B' @ B’ (5) 
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(7) 


If monomers B and C can add to one another, then a tri- 
polymerization experiment is necessary. The equations for 
this case also follow simply from (1) and (2): 


a_[RiAB+RiR2A?+R2AC] [A+e°B+o0°C] 
b [R2bAB+o°B?+a°BC] [RiA+C] ’ 





(8) 


a_[RiAB+RiR:A*+RAC] [At+@B+e°C] 
¢ [RiAC+eBC+ eC] ' 





[R:A+B] 


where Ri=hoal koe; R2=Realke. 


If separate 2-component studies have yielded the values 
of o and a%, plus the information that B and C do not add 
to themselves (ki=k-c=0) but do add to one another 
(B and C form a copolymer), then a tripolymerization ex- 
periment will yield enough information to determine the 
unknown ratios R: and Ro». 

1T. Alfrey and G. Goldfinger, J. Chem. Phys. 12, 322 (1944). 

2 These equations have been extended to the general case of 
TTT C. Walling and E. R. Briggs, J. Am. Chem. Soc. 67, 


3 Private communication. 
4T. Alfrey and E. Lavin, J. Am. Chem. Soc. 67, 2044 (1945). 





